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Executive Summary 



Purpose of this White Paper 

There has been superb progress in understanding the neutrino sector of elemen- 
tary particle physics in the past few years. It is now widely recognized that 
the possibility exists for a rich program of measuring CP violation and matter 
effects in future accelerator v experiments, which has led to intense efforts to 
consider new programs at neutrino superbeams, off-axis detectors, neutrino fac- 
tories and beta beams. However, the possibility of measuring CP violation can 
be fulfilled only if the value of the neutrino mixing parameter ^13 is such that 
sin^(20i3) greater than or equal to on the order of 0.01. The authors of this 
white paper arc an International Working Group of physicists who believe that 
a timely new experiment at a nuclear reactor sensitive to the neutrino mixing 
parameter ^13 in this range has a great opportunity for an exciting discovery, a 
non-zero value to ^13. This would be a compelling next step of this program. 
We are studying possible new reactor experiments at a variety of sites around 
the world, and we have collaborated to prepare this document to advocate this 
idea and describe some of the issues that are involved. 

Purpose of the Experiment 

In the presently accepted paradigm to describe the neutrino sector, there are 
three mixing angles. One is measured by solar neutrinos and the KamLAND 
experiment, one by atmospheric neutrinos and the long-baseline accelerator 
projects. Both angles are large, unlike mixing angles among quarks. The third 
angle, ^13, has not yet been measured to be nonzero but has been constrained 
to be small in comparison by the CHOOZ reactor neutrino experiment. 

The basic feature of a new reactor experiment is to search for energy de- 
pendent disappearance using two (or more) detectors, to see Vf. j>e dis- 
appearance. The detectors need to be located underground in order to reduce 
backgrounds from cosmic rays and cosmic ray induced spallation products. The 
detectors need to be designed identically in order to reduce systematic errors to 
1% or less. Control of the relative detector efficiency, fiducial volume, and good 
energy calibration are needed. 

A measurement of or stringent limit on ^13 would be crucial as part of a 
long term program to measure CP parameters at accelerators, even though a 
reactor — > Ve disappearance experiment docs not measure any CP violating 
parameter. A sufficient value of 6'i3 measured in a reactor experiment would 
strongly motivate the investment required for a new round of accelerator v ex- 
periments. A reactor (experiment's unambiguous measurement of ^13 would also 
strongly support accelerator measurements by helping to resolve degeneracies 
and ambiguities. The combination of measurements from reactors and neutrino 
results from accelerators will allow early probes for CP violation without the 
necessity of long running at accelerators with anti-neutrino beams. 
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Anticipated Sensitivity 

The best current limit on 613 comes from the CHOOZ experiment and is a 
function of Am^^^, which has been measured using atmospheric neutrinos by 
Supcr-Kamiokande and others. The latest reported value of Am^j^ from Super- 
Kamiokande is 1.2 < Am^j^ < 3.0 x 10~'^eV^ with a best fit reported at 
2.0. The CHOOZ limits for Am^^^ of 2.6 and 2.0 x lO'^eV^ are sin^(26'i3) < 
0.14 and 0.20. Global fits using the solar data limit the value for small Am^,„j 
to less than 0.12. In order to improve on the CHOOZ experiment, a new reactor 
experiment needs more statistics and better control of systematic errors. The 
relative sensitivity at low Am^^^ can be improved by locating the far detector 
further than 1 km. Increased statistics can be achieved by running longer, 
using a larger detector, and judicious choice of a nuclear reactor. The dominant 
systematic errors in an absolute measurement of the reactor neutrino flux, such 
as cross-sections, flux uncertainties, and the absolute target volume, will be 
largely eliminated in a relative measurement with two or multiple detectors. 
Good understanding of the relative detector response and the backgrounds is 
required for a precise relative measurement of the reactor neutrino flux and 
spectrum. Experiments are being considered which increase the luminosity from 
the CHOOZ value of 12 t GW y (ton-Gigawatt-years) to 400 t GW y or more. 
This will allow a mixing angle sensitivity of sin'^(20i3) > 0.01. For example, 
400 t GW y would be obtained with a 10 (40) ton far detector, and a 14 (3.5) 
GW reactor in 3 years. One design consideration of the new experiment is 
the possibility for upgrades to achieve even greater luminosity and sensitivity. 
The ability to phase upgrades to achieve a luminosity of 8000 t GW y is being 
considered. 

Major Challenges 

A new reactor experiment will build on the experience of several previous re- 
actor experiments, such as CHOOZ, Palo Verde and KamLAND (described in 
Section 4 of this white paper). These experiments had different goals, mostly 
being designed for signals due to large mixing. Important experience on cali- 
bration, control of systematic errors and the reduction of background has also 
been obtained by the Super-Kamiokande, SNO and Borexino collaborations. 

A next-generation reactor experiment will be designed to make a precision 
measurement of the reactor electron anti-neutrino survival probability at dif- 
ferent distances from the reactor and search for subdominant oscillation effects 
associated with the mass splitting of the mi and m3 mass cigcnstates. A mea- 
surement at the 0(1%) level will require careful control of possible systematic 
errors. Most of the technical requirements of this experiment are well under- 
stood but the details of the detector design still need to be optimized. Some 
of the open questions under consideration are the following: liquid scintillator 
loaded with 0.1% of gadolinium has been used in the past, but there are concerns 
regarding its stability in solution and possible attenuation length degradation 
which need to be fully understood. If movable detectors are chosen, there must 
be confidence that moving the detector does not introduce additional time- 
dependent effects. The use of a second detector will certainly help to control 
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many systematic errors, but also will present a challenge in maintaining a known 
relative calibration over time. Another challenge is reduction of cosmic ray as- 
sociated backgrounds such as neutrons and ^Li spallation and their accurate 
estimation. The reduction of gamma ray background is also important because 
it will affect the ability to reduce the threshold to below 1 MeV. These and 
other design issues are discussed in Sections 5-8 of this white paper. 

Experimental Prospects 

The International Working Group has held two workshops (April 30-May 1, 
2003 at the University of Alabama and October 9-11, 2003 at Technical Univer- 
sity of Munich) and we are planning a third one (March 20-22, 2004 at Niigata 
University.) During the past year, the International Working Group has iden- 
tified a large number of reactors as possible sites for a new experiment. Many 
of these sites are discussed in Section 9, and a few are described in more detail 
in seven Appendices. These include the Angra reactor in Brazil; the possibility 
of a new experiment at CHOOZ, called Double-CHOOZ (or CH6li36'i3Z); Daya 
Bay near Hong Kong in China, Diablo Canyon in California; a reactor in Illi- 
nois; the reactor complex at Kashiwazaki in Japan, and the Krasnoyarsk reactor 
underground at Zhclczhnogorsk in Russia. 

It is not the role of this document to provide a cost estimate or schedule for 
any of the experiments which will be proposed. But it is appropriate to try to 
set the scale of the endeavor in order to compare to other kinds of initiatives in 
neutrino physics. A two-detector system as described in this document seems 
to cost in the range $5M to $15M. The civil construction costs to place these 
detectors underground will be very site dependent and require a detailed en- 
gineering cost estimate as described in Section 11. Estimates are in the range 
of several tens of millions of dollars, depending on site condition and tunnel 
length. Since reactors with an underground site already exist, such as those at 
CHOOZ and Krasnoyarsk, there is a strong incentive to consider those sites for 
the earliest experiment, though there may be physics trade-offs which must be 
considered. Some of the envisioned reactor experiments might start taking data 
in 2007-2008. First results could be achieved as early as 2009. 

None of these efforts has yet resulted in a proposal to a funding agency, but 
site specific proposals and R&D proposals will be submitted during 2004. This 
white paper is a step in that direction. Given the importance of the measurement 
of ^13 and the enthusiasm of the proponents, we are hopeful that two or more 
of these experiments will move forward on a favorable time scale. 
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1 Introduction 



The discovery of neutrino oscillations is a direct indication of physics beyond 
the Standard Model and it provides a unique new window to explore physics 
at high mass scale including unification, fiavor dynamics, and extra dimensions. 
The smallness of neutrino masses and the large lepton flavor violation associated 
with neutrino mixing are both fundamental properties that give insights into 
modifications of current theories. Other possibilities that may reveal themselves 
in the neutrino sector include extra "sterile" neutrinos, CP violation in the 
neutrino mixing matrix, and CPT violation associated with the neutrino mass 
hierarchy. Since neutrino oscillations have now been established, the next step 
is to map out the parameters associated with neutrino masses and mixings. 
The experimental program to accomplish this goal will require a wide range of 
experiments using neutrinos from solar, atmospheric, reactor, and accelerator 
sources. Due to the relations between these various measurements, it will be 
important for the world-wide community to set up a structured program to work 
through the experimental measurements in a coherent and logical manner. 

The existing experimental results fit rather nicely into a picture with three 
massive neutrinos, which corresponds to the simplest scenario for three gen- 
erations (for recent global analyses see, e.^.. References mEI). Neutrino os- 
cillations then involve two mass-squared differences {Am^i and Ato^jj where 
Amfj = m{uiY — m{iyj)'^), three mixing angles (0i2, 023, and Oia), and a CP- 
violating phase ((5). The present status of these parameters is summarized in 
Figure ^ Atmospheric neutrino data [S] and the first results from the K2K 
long-baseline accelerator experiment ^ determine |Am|2| = (Slj g) x 10""^ eV^ 
(errors at 3a) and 6*23 ~ 45° Ej, whereas most solar data EIj com- 
bined with the results from the KamLAND reactor experiment [7], lead to 
Amii = (6.9+1;^) X 10-5 eV2 and sin2 6'i2 = O-ito ^or at 3cr p. 

The neutrino sector may contain more than three neutrinos by including 
mixing to sterile neutrinos (for example to account for the LSND |S] anomaly), 
but in these cases the mixing matrix most likely factors to a good approximation 
into a (3 X 3) submatrix with the parameters given above. The investigation 
of oscillations involving sterile neutrinos will demand measurements such as 
MiniBooNE as well as improved disappearance measurements at high Am^. 

The current experimental situation can thus be summarized by two more or 
less decoupled oscillations governed by the "atmospheric" and "solar/reactor" 
quadratic mass splittings Am^^^ — Am^g and Am^^^ = Am^i, respectively, 
and the corresponding mixing angles 612 = Ogoi and 623 = 9atm, which turned 
out to be surprisingly large. This leads in the future to two equally impor- 
tant experimental directions: The first task is to improve the knowledge of the 
above (leading) oscillation parameters and to make precision measurements. 
Conceptually at least equally important is the fact that three flavors imply also 
three flavor oscillations and thus one further mixing angle, 613 as well as a CP 
violating phase 6 ^. The CP phase (5 is a very interesting, but so far a com- 

^Note that Majorana neutrinos imply also two further CP violating phases, but these do 
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Figure 1: Status of neutrino oseillation parameters from a combined analysis of 
current global data 1 . Left panel: allowed regions of solar (Am^i, sin^ 6*12 ) and 
atmospheric {Am^i,sm^ 623) parameters at 90%, 95%, 99% and 3cr CL. Right 
panel: upper bound on sin^ 26*13 from the CHOOZ experiment at 90% (dashed) 
and 3cr (solid) CL for 1 DOF as a function of Am|]^. The light (dark) shaded 
region is excluded from CHOOZ + solar + KamLAND data at 90% (3cr) CL 
for 1 DOF. The horizontal lines indicate the current best fit value and the 3a 
allowed regions for Am|]^. 



pletely unknown, parameter. The fact the LMA solution has been confirmed 
means that S is in principle accessible in future experiments if ^13 is not too 
small. In many models of neutrino masses the see-saw mechanism leads to con- 
nections of the leptonic CP phase 6 to the CP phases in the heavy Majorana 
sector and thus to leptogenesis, one of the best known mechanism to explain 
the baryon asymmetry of the universe (see e.g. [5]). Neutrino masses may 
therefore explain a second indication for physics beyond the Standard Model, 
since the observed baryon asymmetry cannot be generated from CP violation in 
the Standard Model with massless neutrinos. Future neutrino experiments aim 
therefore indirectly at another key question in physics, namely what causes the 
baryon asymmetry in the Universe. 

The mixing angle ^13, the parameter relevant for three flavor effects in neu- 
trino oscillations, is known to be small from the CHOOZ 10, 11 and also from 
the Palo Verde experiment 12 . The current bound from global data is summa- 
rized in the right hand panel of Figure^ It depends somewhat on the true value 
of the atmospheric mass squared difference, since the bound from the CHOOZ 
experiment gets rather weak for Amgj^ ^ 2 x 10^'^ eV^. However, in that re- 
gion an additional constraint on ^13 from global solar neutrino data becomes 
important J^. At the current best fit value of Arri^i = 2 x 10"'^ eV^ we have 

not enter into neutrino oscillations. 
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the bounds at 90% (3cr) CL for 1 DOF 

sin^ 26*13 < 0.16 (0.25) , sin^ 613 < 0.053 (0.066) , 613 < 10.8° (14.9°) . (1) 

Genuine three flavor oscillation effects occur only for a finite value of 613 and 
establishing a finite value of 613 is therefore one of the next milestones in neutrino 
physics. Leptonic CP violation is also a three flavor effect, but it can only be 
tested if 6*13 is finite. There is thus a very strong motivation to establish a finite 
value of 613 in order to aim in the long run at a measurement of leptonic CP 
violation (see e.g. 

Future measurements of disappearance using a two detector reactor ex- 
periment and long-baseline i'^ — > ly^ and j/^ v^, experiments will be crucial 
in determining ^13, the sign of Am^j, and the CP phase 8. If 6*13 > 0.01, the 
design of experiments to measure the sign of Am|2 and the CP phase b become 
straight forward extensions of current experiments. For this reason, there is 
general agreement that a B\3 measurement should be the prime goal of the next 
round of experiments. On the theoretical side these experiments could test if 
the small value of ^13 could be a numerical coincidence or if e.g. some symmetry 
argument is required to explain a tiny value. 

Future measurements of Q\3 are possible using reactor neutrinos and ac- 
celerator neutrino beams. As will be shown in subsequent sections, reactor 
measurements have the property of determining ^13 without the ambiguities as- 
sociated with matter effects and CP violation. In addition, the needed detector 
for an initial reactor measurement is small (< 50 tons) and the construction of a 
neutrino beam is not necessary. For this reason, a precision reactor experiment 
could lead the way in establishing the future oscillation program by setting the 
scale of the d\3 mixing angle. The previous most accurate measurements were 
by the CHOOZ and Palo Verde experiments where a single detector was placed 
about 1 km from the reactor. Future reactor experiments using two detectors 
(~ 50 tons) at near (100 - 200 m) and far (1-2 km) locations will have signifi- 
cantly improved sensitivity for Q\3 down to the 0.01 level. With Q\3 determined, 
measurements of —f and i/^ oscillations using accelerator neutrino 

beams impinging on large detectors at long baselines will improve the knowledge 
of 6*13 and also allow access to matter or CP violation effects. For the field to ex- 
ploit the physics opportunities available for neutrino oscillation measurements, 
it is clear that a suite of experiments including both reactor and long-baseline 
accelerator measurements will be necessary. 

In addition to the general physics arguments, there are two factors that lend 
urgency to this initiative. Our studies indicate that a reactor experiment to 
measure sin^ 2^13 to the level of 0.01 could be done at significantly less cost and 
on a more rapid time scale than an accelerator long-baseline neutrino experiment 
with comparable sensitivity. This conclusion is infiuenced by several recent de- 
velopments, including the High Energy Physics Roadmap for the future, Febru- 
ary 2002 , the prioritizations made in "High Energy Physics Facilities on the 
DOE Office of Science Twenty- Year Roadmap" issued by the U.S. Department 
of Energy in March, 2003 17 and the "Facilities for the Future of Science: a 
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Twenty- year Outlook" issued by the Office of Science of the DOE in November, 
2003 ^S]. In particular, the latter document envisions that a high-intensity 
neutrino beam is more than 15 years in the future. For comparison with an 
off-axis long baseline experiment, we use cost and time estimates based on cur- 
rent work for the Fermilab proposal P929, the NuMI Off-Axis Experiment. We 
emphasize that a new reactor experiment does not reduce the motivation for 
the latter experiment; it obtains information complementary to that obtained 
by the reactor 613 experiment, such as the mass hierarchy between m^vg) and 
to(i/i). Instead, we would envision that, since the reactor experiment can be 
performed more quickly, its findings concerning ^13 will provide very important 
guidance for the long baseline program. 

In this White Paper, we outline the capabilities of next generation reactor 
experiments and summarize the design considerations that groups are consid- 
ering in developing this program. The International Working Group on ^13 is 
sharing ideas on how to best design a new reactor experiment, and one goal 
of this White Paper is to document the present status of our understanding of 
these issues. 

In the next section, we discuss in more detail the physics opportunities and 
the motivation for a new reactor experiment. The following Section |21 deals 
with the optimal baseline, luminosity scaling and the impact of systematic er- 
rors. Previous reactor experiments are described in Sectional and in Section |S1 
we present some thoughts about the general layout of the detector, a multi- 
layered volume of scintillator designed to define the fiducial volume well, and 
also carefully control other potential systematic errors. In Sectional the cal- 
ibration requirements for the detector are reviewed. Section 13 considers the 
issues of backgrounds and how they affect the required overburden. Depths 
that provide an overburden of 400 mwe to 1100 mwe are desirable. The goal of 
carefully minimizing systematic errors is qualitatively different than has been 
required of neutrino experiments at reactors in the past. We are confident that 
the two detector concept will provide lower systematic errors than have been 
previously achieved, but the ultimate limit on achievable systematic error has 
yet to be identified. A discussion of a variety of systematic errors is presented 
in Section|Sl Characteristics of a large number of sites are reviewed in Section|51 
and some more detailed experimental site plans for seven of the possible loca- 
tions are included in the Appendices to this document. Next we discuss other 
physics that can be done in Sectional Depending on the site, the costs of a new 
reactor experiment will potentially be dominated by the civil construction of a 
shaft or tunnel. Those civil engineering issues are reviewed in Sect ion [TTl Safety 
issues are discussed in Sectional Section is finally devoted to outreach and 
educational issues. The appendices contain further details of potential sites. 
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2 Physics Opportunities and Motivation 



2.1 Road Map for Future Neutrino Oscillation Measure- 
ments 

There is now a world-wide experimental program underway to measure the pa- 
rameters associated with neutrino oscillations. The current experiments include 
K2K that measures disappearance over a 250 km baseline from KEK to SK. 
Another experiment is MiniBooNE that is searching for appearance 
signal in the LSND Am^ region from 0.2 to 1 eV^. Upcoming longer-baseline 
(~ 700 km) experiments are NuMI/MINOS at Fermilab and CNGS at CERN 
that will study oscillations in the atmospheric Am? region. Groups in all the 
world-wide regions are also pursuing sites and experiments for a precision reactor 
experiment using detectors with fiducial volumes of 5 to 50 ton. Several near- 
term new long-baseline experiments are planned which will use off-axis beams 
including the approved J-PARC (previously called JHF) to Super-K (22.5 kton) 
experiment and the developing NuMI off-axis experiment (50 kton detector). 
Following these experiments, the next stage might be neutrino superbeam ex- 
periments with even longer baselines that could possibly be combined with large 
proton decay detectors. Four such projects under consideration are: (i) BNL 
with an AGS upgrade, (ii) Fermilab with a proton driver upgrade, (iii) J-PARC 
(phase II), and (iv) a CERN Superconducting Proton LIN AC experiment. Fu- 
ture neutrino factories, using a muon storage ring, will provide the ultimate in 
sensitivity and precision in oscillation measurements. 

It is clear that developments in the field will dictate how the community 
should proceed through these studies. As stated previously, the size of ^13 is 
the small parameter that sets the scale for fmtlicr studies in a three neutrino 
scenario. It is also clear that the final resolution of the LSND anomaly by 
MiniBooNE could significantly affect the direction for new investigations. To 
bring this information together in a coherent way, we present a roadmap for 
neutrino oscillations which tries to point out the relations between the various 
measurements: 

• Stage 0: The Current Program 

— There are improved measurements of Ami2 (5-10%) by solar neutrino 
and the KamLAND experiments. 

— NuMI, CNGS, and K2K experiments check the atmospheric oscilla- 
tion phenomenology and measure Am^sto ^ 10%. 

— MiniBooNE makes a definitive check of the LSND effect and measures 
the associated Am^ if the effect is confirmed. 

• Stage 1: Measurement or tight constraint's on the ^13 angle^ 



^The combination of all these experiments may give the first indications of matter and CP 
violation effects. 
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— The NuMI/MINOS on-axis experiment probes sin^ 2^13 > 0.06 at 
90% CL. 

— Two-detector, long-baseline reactor experiments probe sin^ 2^13 > 
0.01 at 90% CL. 

— The NuMI and J-PARC off-axis experiments with 20-50 kton de- 
tectors investigate —^ Ve transitions for oscillation probabilities 
greater than 1%. 

• Stage 2: Measurements of the sign of Am23 and CP violation using su- 
perbeams and very large detectors (500 to 1000 kton) 

(This is feasible if sin^ 26'i3 > 0.01 and if 6 is large enough.) 

— Measurements of i/^ —t v^. at several baselines need to be combined 
with either precision reactor measurements of v^, — > or with i?^ — > 

— Increased neutrino beam rates are needed, especially for the run- 
ning, which make high intensity proton sources necessary. 

• Stage 3: Measurements with a Neutrino Factory 

— New facilities probe a mix of ^v^ ^j,, — s- ^i^'e/p transitions with sensi- 
tivities below the 0.001 level 

— They also map out CP violation with precision for sin'^ 26*13 > 0.001. 
A flow chart with these ideas is shown in Figure |21 

2.2 Where do reactor oscillation experiments fit in? 

Any oscillation effect in Ve survival is governed, assuming three flavor mixing, 
by the equation 

P{Ve ^ve) = l~ sin^ 2013 sin^( ^";°*"'^ ) - cos^ ^13 sin^ 26,2 sm'{^^). 

(2) 

This equation is plotted in Figure |31 as a function of L/E with the current best 
values for the Am^s and mixing angles (sin^(20i3) is set to the maximum value 
allowed by current limits). One can clearly see the two oscillations governed 
by the two Am^s. Experimentally, a judicious choice of L/E should be able to 
distinguish one effect from the other. The KamLAND experiment is the first 
reactor experiment to see oscillation effects, by measuring a 40% disappearance 
of Pe- Given that the average baseline for KamLAND is 180 km, the detected 
deficit is presumably associated with the third (Am^2) term in Equation (jj)). 

The current best limit on 6*13 comes from a lack of observed oscillations at 
CHOOZ and Palo Verde (sin^ 26*13 < 0.20 for Amlt^ = 2.0 x lO-^eV^). These 
experiments were at a baseline distance of about 1 km and thus more sensitive 
to the second (Am^^^) term of Equation Those experiments could not 
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have had greatly improved sensitivity to ^13 because of uncertainties related 
to knowledge of the flux of neutrinos from the reactors. They were designed to 
test whether the atmospheric neutrino anomaly might have been due to — > Ue 
oscillations, and hence were searching for large oscillation effects. 

Since the effective disappearance will be very small (see Figure any new 
experiment which is designed to look for non-zero values of ^13 would need to 
move beyond the previous systematic limitations. This could be achieved by 
utilizing the following properties: 

• two or more detectors to reduce uncertainties to the reactor flux 

• identical detectors to reduce systematic errors related to detector accep- 
tance 

• carefully controlled energy calibration 

• low backgrounds and/or reactor-off data 

Note that CP violation does not affect a disappearance experiment, and that 
the short baseline distances involved in a reactor measurement of oscillations at 
the atmospheric Am^ allow us to safely ignore matter effects. 

A next generation reactor oscillation experiment would use at least two de- 
tectors placed at various distances from a high power reactor (Figure 01). The 
reactor provides a high intensity, isotropic source of neutrinos with a well-known 
spectrum as shown in Figure El The neutrino cross section for this process is 
well known as described in Reference |19j . Antineutrinos are detected through 
the inverse-/? decay process followed by neutron capture. 

i/g + P ^ -|- n. 

The detector would most likely be composed of a vat of scintillator oil viewed 
from its surface by an array of photomultipliers. In order to reduce the back- 
ground from cosmic-ray spallation, the detectors will need to be underground 
with at least 300 mwe of shielding. A detected event would correspond to a 
coincidence signal of an electron and capture neutron. The incident neutrino 
energy is directly related to the measured energy of the outgoing electron. The 
search for oscillations would then involve comparing the neutrino rate in the 
two detectors and looking for a non-l/r^ dependence. 

As stated above, ^13 is a key parameter in developing the future neutrino os- 
cillation program. Reactor experiments offer a straightforward and cost effective 
method to measure or constrain the value of this parameter. The sensitivity of 
a two detector experiment is comparable to that of the proposed initial off-axis 
long-baseline experiment. Since a reactor experiment would be much smaller 
and use an existing reactor neutrino source with a well understood neutrino 
rate, the experiment should be able to be done fairly quickly and at reduced 
costs. It is likely that an early measurement of 013 will be necessary before 
the community invests a large amount of resources for a full off-axis measure- 
ment. For the longer term, a reactor experiment would be complementary to the 
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off-axis experiments in separating the measurement of from other physics 
parameters associated with matter effects and CP violation. A follow-up reactor 
experiment with much larger detectors at various baselines will continue to be 
an important component of the neutrino oscillation program. 
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Figure 2: Flow chart for a possible future neutrino program 
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Figure 3: Probability of disappearance versus L/E for ^13 at its current upper 
limit 
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Figure 4: Schematic layout of a two detector reactor neutrino oscillation exper- 
iment. 
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Figure 5: Antineutrino flux, cross section, and relative event rate using a reactor 
source pi] . 
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2.3 Reactor experiment as a clean laboratory for the ^13- 
measurement 



In this section, we demonstrate that a reactor measurement of ^13 is a clean 
measurement which is free from any contamination, such as from effects of the 
other mixing parameters or from the Earth matter effect [Mj- This key feature 
is one of the most important advantages of the reactor experiments. We use the 
standard notation of the lepton flavor mixing matrix: 



U 



C12C13 S12C13 sise'"-^ 

-S12C23 - Ci2S23Sl3e"^ C12C23 - 5125235136"^ S23C13 
S12S23 - £126235136"^ -C12S23 - Sl2C23Sl3e"^ C23C13 



(3) 



Due to the low neutrino energy of a few MeV, the reactor experiments are 
inherently disappearance experiments, i.e., they can only measure the survival 
probability P{Pe — > i^g). Unlike the case of the Ve appearance probability, it is 
well-known that the survival probability does not depend on the CP phase S in 
arbitrary matter densities (for more details, see |H2). For reactor experiments, 
the matter effect is very small because the energy is quite low and the effect can 
be ignored to a good approximation. This can be seen by the comparison of the 
matter and the vacuum effects 

-^^34.10-^r '^"^3il VV^U p UM 

V2.5-10-3eVV V4MevA2.8g-cm-3A0.5y' ■ ^ ^ 

Here E is the neutrino energy. In addition, a — \/2Gi?A^e denotes the index 
of refraction in matter with the Fermi constant Gf and the electron number 
density Ne in the Earth (which is related to the Earth matter density p by 
-^e = Yep/mp with the proton fraction Ye). 

Since we know that the matter effect is negligible, we immediately under- 
stand that the survival probability is independent of the sign of Am§]^. There- 
fore, one can use the vacuum probability formula for the analysis of a reactor 
measurement of ^13. The expression for P{i>e — > i^e) in vacuum is given by j21| 

1 - P{De ^De) = sin^ 26li3 sin^ A31 

+ 1/2 c^2sin^ 26*13 sin2A3i sin2A2i 

-I- Ci3 sin^ 26^12 sin^ A21 

-I- sin^ 26*13 cos 2 A31 sin^ A21 , (5) 

where Aij = AmfjL/(AE) and Cy = cos 9ij. Defining the mass hierarchy param- 
eter a as a = Amli/Arrigi, where |a| ~ 0.03, the second term in Equation ^ 
is suppressed relative to the main depletion term (first term) by a factor of 
\a\ sin^ 2013 < 6 x 10"^^ the fourth term even by a factor of sin^ 2^13. Thus, 
we can re-write Equation ((Sj for A21 <C 1 (for the baselines considered) as 

l-P(Pe-^Pe) ^ sin^ 26*13 sin^Aai +Q!^A3iCi3 sin^ 26*12. (6) 
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Though the second term on the right-hand side of this equation could be of 
the order of the first term for very large \a\, it can be neglected for the first 
atmospheric oscillation maximum (where the first term is large) and sin^ 26*13 
larger than about 10"'^. Therefore, the disappearance probability can be well 
approximated by the two-flavor depletion term in vacuum, which is the first 
term in Equation Assuming that |Am§i| is accurately determined by a long- 
baseline disappearance measurement, the reactor experiments thus serve for 
a clean measurement of ^13 independent of other mixing parameters. 

2.4 Comparison to superbeams 

We have demonstrated in the last section that reactor measurements allow a 
degeneracy- free measurement of sin^ 2^13 . In order to qualitatively discuss the 
difference between reactor experiments and superbeams, we can compare the 
oscillation probabilities of the dominant oscillation channels. For the super- 
beams, one can expand the appearance probability P^e (or Pps) in terms of the 
small mass hierarchy parameter a = Am^i/Am'^^ and the small mixing angle 
sin 2^13 using the standard parameterization of the leptonic mixing matrix U 
in Equation (O. As a first approximation for a qualitative discussion, one can 
use the vacuum formula from References |14[ 1221 123, with the terms up to the 
second order (i.e., proportional to sin'^ 2^13, sin 2^13 • a, and a^): 

P^e |sin26li3 sin6l23 sin A3ie'('^^^±''<=p) -h cos 6*13 cosfe sin26li2 sin Asij^ 

« sin^ 2013 sin^ 6*23 sin^ A31 

=Fa sin 2^13 sin(5cp cos ^13 sin 26*12 sin 26*23 A31 sin^ A31 

+a sin 2^13 cos Sqp cos 6*13 sin 26i2 sin 2^23 A31 cos A31 sin A31 

cos^ 6I23 sin^ 26*12 . (7) 

Here Ay = J\mf-L/{AE) = {mf — m'^)L/{AE) and the sign of the second term 
refers to neutrinos (minus) or antineutrinos (plus). We have used the approxi- 
mations that sinA2i ~ q;A3i ^ 1 and that A32 ~ A31. 

For the reactor experiments, we have, up to the same order in sin 2^13 and 
a. Equation ©. Comparing Equation ((HJ to Equation jT)) clearly demonstrates 
that the superbeams are quite rich in physics and much more complex to an- 
alyze. Depending on the true values of a and sin 26*13, each of the individual 
terms in Equation Q obtains a relative weight. The result is then determined 
by the mutual interaction of the four terms in Equation Q leading to multi- 
parameter correlations and degeneracies. Correlations and degeneracies are de- 
generate solutions in parameter space, where the correlations are connected 
solutions and the degeneracies are disconnected solutions in parameter space 
(at the chosen confidence level). For example, many of the degeneracy prob- 
lems originate in the summation of the four terms in Equation Q especially 
for large a and sin 2^13, since changes of one parameter value can be often 
compensated by adjusting another one in a different term. This leads to the 
well-known {5,9i^) [S^j sign(Am|;^) and (6*23, 7r/2 — ^23) 1^ degeneracies. 
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Figure 6: The sin^ 26*13 sensitivity limit as function of the atmospheric (left) 
and solar (right) mass squared differences for the JPARC-SK (black curves) 
and Reactor-I (gray curves) experiments from References |29[ I3U| (five years of 
neutrino running for JPARC-SK and an integrated luminosity for Reactor-I of 
400tGWy). In addition, JPARC-SK@443 km is shown for a modified basehne 
of 1.5 X 295 km ~ 443 km in the left-hand plot (dashed curve). For the oscillation 
parameters, Am|i = 2.0 • 10"^ gy^^ sin^ 26*23 = 1-0, Am|i = 7.5 • 10"^ eV^ and 
sin^ 2(?i2 = 0.8 are used (if not varied) j31[ IT]. Furthermore, the Sa excluded 
regions are gray shaded and the current best-fit values are marked. The anal- 
ysis includes systematics, multi-parameter correlations, and degeneracies as in 

Refs. EaEni. 



i.e., an overall "eight-fold" degeneracy which can severely affect the poten- 
tial of many experiment types [2H|- On the other hand, the reactor Equation ^ 
contains the product sin^ 29i3 • sin^ A31 as the main contribution, which leads 
to a simple two-parameter correlation between sin^ 2^13 and sin^ A31. In this 
correlation, sin^ 26*13 acts as the (energy independent) amplitude of the mod- 
ulation and sin^ A31 contains the spectral information. Thus, with sufficiently 
good spectral information and the current knowledge about Arri^i, it is easy to 
disentangle these two parameters. In addition, the reactor measurement hardly 
depends on the true value of Am^i. 

The dependence on the true values of the atmospheric and solar mass squared 
differences is, for the current best-fit values and ranges, illustrated in Figure El 
The figure compares a reactor experiment with an integrated luminosity of 
400 1 GW y to the JPARC to Super-Kamiokande first-generation superbeam ex- 
periment with a running time of five years (neutrino running only). The two 
plots illustrate that the considered reactor experiment would be better than 
the JPARC-SK superbeam at the current best-fit values of Am§i and Am|i, 
as well as in most of the still allowed parameter ranges. Since the energy spec- 
trum of a reactor experiment is broader than the one of an off-axis (narrow 
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band) superbeam, the reactor experiment is less affected by the smaller value 
of Am^i after the Super-Kamiokande re-analysis |33| . In addition, it is hardly 
affected by the true value of Am^i as we discussed above. In the left plot of 
Figure ini we also show the JPARC-SK experiment for the (artificial) baseline 
of 1.5 X 295 km ~ 443 km, which means that for this baseline the oscillation 
peak is shifted to Amj^ = 2.0 • lO^^gy^. The fi gure clearly demonstrates that 
this shifting would not solve the problem due to the 1/L^ luminosity scaling. 
Because of the over-proportional loss of events for a lower neutrino energy due 
to the production mechanism and the cross section energy dependence, a lower 
energy instead of the longer baseline would also not help. 

We have now demonstrated that the reactor measurement could provide a 
more robust limit on sin^ 29i3 with respect to the (within certain ranges) true 
parameter values of Ato§]^ and Amfj^. However, it is obvious from Equation JB)) 
that reactor experiments at a baseline of a few kilometers are not sensitive to 
the mass hierarchy or 6cp , which means that superbeams will still be needed to 
test these parameters. On the other hand, a large reactor experiment could help 
to resolve the degeneracies in Equation by measuring sin^ 2^13. In this case, 
one could talk about synergies between a reactor experiment and a superbeam. 
For example, it has been demonstrated in Reference [SOj that there are several 
advantages from a large reactor experiment (e.g., with an integrated luminosity 
of 8000 1 GWy). First of all, a reactor experiment could help to determine the 
mass hierarchy very well independently of the true value of Am|^. Secondly 
it could improve the CP sensitivity by allowing to operate the superbeam with 
neutrinos only (instead of using a fraction of antineutrino running). A reactor 
experiment performed on a shorter timescale than a superbeam would change 
the main goal of a superbeam from finding a non-zero value for sin^ 2^13 to 
measuring 5 and the sign of Arn^i- 
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2.5 Theoretical Motivation for non-zero ^13 

One may ask if there exist theoretical reasons why ^13 should be within the reach 
of a new experiment, with a sensitivity down to sin^ 2^13 ~ 0.01. This question 
is of course connected to the origin of neutrino masses. For example, there exist 
apparent regularities in the fermionic field content which make it very tempting 
to introduce right-handed neutrino fields leading to both Dirac and Majorana 
mass terms for neutrinos. The diagonalization of the resulting mass matrices 
yields Majorana mass eigenstates and generically very small neutrino masses. 
This is the well known see-saw mechanism . It can be nicely accommodated 
in embeddings of the SM into a larger gauge symmetry, such as SO(IO). 

A reason for expecting a particular value of 613 does clearly not exist as long 
as one extends the SM only minimally to accommodate neutrino masses. 613 
is then simply some unknown parameter which could take an arbitrarily small 
value, including zero. The situation changes in models of neutrino masses. 
Even then one should acknowledge that in principle any value of ^13 can be 
accommodated. Indeed, before the discovery of large leptonic mixing, many 
theorists who did consider lepton mixing expected it to be similar to quark 
mixing, characterized by small mixing angles. Experiment led theory in showing 
the striking results that sin^ 2^23 — 1 and tan^ 6*12 ~ 0.44, while 613 is small. 
Indeed, the most remarkable property of leptonic mixing is that two angles 
are large. Therefore, today there is no particular reason to expect the third 
angle, ^13, to be extremely small or even zero. This can be seen in neutrino 
mass models which are able to predict a large 612 and 023 . They often have 
a tendency to predict also a sizable value of 613. This is both the case for 
models in the framework of Grand Unified Theories and for models using flavor 
symmetries. There exist also many different texture models of neutrino masses 
and mixings, which accommodate existing data and try to predict the missing 
information by assuming certain elements of the mass matrix to be either zero or 
equal. Again one finds typically a value for 613 which is not too far from current 
experimental bounds. A similar behavior is found in so-called "anarchic mass 
matrices" . Starting essentially with random neutrino mass matrix elements one 
finds that large mixings are actually quite natural. 

An overview of various predictions is given in Tabled For more extensive 
reviews, see for example |SSI ESI EH EH) • The conclusion from all these consid- 
erations about neutrino mass models is that a value of 6*13 close to the CHOOZ 
bound would be quite natural, while smaller values become harder and harder 
to understand as the limit on 9i3 is improved. 

Besides, neutrino masses and mixing parameters are subject to quantum 
corrections between low scales, where measurements are performed, and high 
scales where some theory predicts ^13. Even in the "worst case" scenario, where 
9i3 is predicted to be exactly zero, they cause ^13 to run to a finite value at 
low energy. Strictly speaking, ^13 = cannot be excluded completely by this 
argument, as the high-energy value could be just as large as the change due 
to running and of opposite sign. However, a severe cancellation of this kind 
would be unnatural, since the physics generating the value at high energy are 
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not related to those responsible for the quantum corrections. The strength of 
the running of ^13 depends on the neutrino mass spectrum and whether or not 
supersymmetry is realized. For the Minimal Supersymmetric Standard Model 
one finds a shift Asin^ 2^13 > 0.01 for a considerable parameter range, i.e. one 
would expect to measure a finite value of ^13 [3^] . Conversely, limits on model 
parameters would be obtained if an experiment were to set an upper bound on 
sin^ 2013 in the range of 0.01. In any case, it should be clear that a precision 
of the order of quantum corrections to neutrino masses and mixings is very 
interesting in a number of ways. 

Altogether there exist very good reasons to push the sensitivity limit from 
the current CHOOZ value by an order of magnitude and to hope that a finite 
value of 9i3 will be found. But as already mentioned, at this precision even a 
negative result would be very interesting, since it would test or rule out many 
neutrino mass models and restrict parameters relevant for quantum corrections 
to masses and mixings. From a larger point of view the experiments discussed 
in this white paper might probe if a small value of ^13 is a numerical coincidence 
or the result of some underlying symmetry. 



17 



Keierence 


sm b'ls 


sm 20i3 


d(J(1uJ 






Lron, Monapatra, iNg |4Uf 


A 1 O 


O.lo 


Urbijola bU(lO) 






Asaka, Buclimiiller, Covi |41| 


U.l 


U.U4 


SO(IO) + flavor symmetry 
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0.18 


Albright, Barr gS] 
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Chen, Mahanthappa IIHI 
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0.1 
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Buchmiiller, Wyler [KT)] 


0.1 


0.04 


Bando, Obara EH 


0.01 .. 0.06 


4- 10-4 .. 0.01 



Flavor symmetries 



Grimus, Lavoura [S21ISS! 








Grimus, Lavoura [S^I 


0.3 


0.3 


Babu, Ma, Valle ^ 


0.14 


0.08 


Kuchimanchi, Mohapatra [ES] 


0.08 .. 0.4 


0.03 .. 0.5 


Ohlsson, Seidl ^ 


0.07 .. 0.14 


0.02 .. 0.08 


King, Ross [S7j 


0.2 


0.15 


Textures 






Honda, Kaneko, Tanimoto (SH] 


0.08 .. 0.20 


0.03 .. 0.15 


Lebed, Martin P?) 


0.1 


0.04 


Bando, Kaneko, Obara, Tanimoto 60 


0.01 .. 0.05 


4- 10-4 .. 0.01 


Ibarra, Ross 61 


0.2 


0.15 


3x2 see-saw 






Appelquist, Piai, Shrock 


0.05 


0.01 


Frampton, Glashow, Yanagida [Hil 


0.1 


0.04 


Mei, Xing |65j (normal hierarchy) 


0.07 


0.02 


(inverted hierarchy) 


> 0.006 


> 1.6 • 10-4 


Anarchy 

dc Gouvea, Murayama 66 


> 0.1 


> 0.04 


Renormalization group enhancement 
Mohapatra, Parida, Rajasekaran 


0.08 .. 0.1 


0.03 .. 0.04 



Table 1: Incomplete selection of predictions for 6*13. The numbers should be 
considered as order of magnitude statements. 
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3 Optimal Baseline Distances, Luminosity Scal- 
ing, and the Impact of Systematics 

3.1 Total Flux vs. Baseline 

The equation for the survival probabiUty of reactor neutrinos under full three 
flavor mixing was previously described in Equation ||2Jl. It was pointed out by 
reference to Figure |31 that a judicious choice of baseline distance could restrict 
one to oscillations dominated by one or the other of the Am^ oscillations. For 
this experiment, we are choosing to focus on the shorter baseline, which cor- 
responds to Ato^(„j. Thus, neglecting the other oscillation term. Equation |2J) 
reduces to 

P{v, ^D,) = l- sin2(20i3) sin2(Am2,„L/4£;) (8) 

where the current best estimate from Super-Kamiokande has Am^^^ = 0.002 
|33|. To measure this oscillation effect, the optimal choice in baseline distance 
will depend on the energy of the neutrinos. As shown previously in Figure |31 
the detected spectrum for reactor neutrinos is between 1-10 MeV with a peak 
at about 3.8 MeV. In addition, recall that the flux of neutrinos will fall as 
the square of the baseline distance. A comparison of the expected flux with 
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Figure 7: Illustrative expected number of detected neutrino events as a function 
of baseline distance from the reactor core. The two curves show the expectations 
for the case of no oscillations and for an oscillation with Am^ = 0.002 and an 
amplitude of sin^(20i3) = 0.5 which is 2.5 times the current limit. The curves 
are calculated for a luminosity of 600 1 GW y and a mono-energetic neutrino flux 
at 3.8 MeV. 

and without oscillations is shown in Figure |7| for a mono-energetic neutrino 
beam of 3.8 MeV. Note that the amplitude of the oscillation shown is set to 
sin^(20i3) — 0.5 which is 2.5 times the current limit from CHOOZ in order to 
amplify the effect. As would be expected, one sees the disappearance effect at 
regular intervals. However, when the full energy spectrum, shown in FigureEl is 
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folded in, the regular disappearance effect is washed out even with the magnified 
amplitude (see Figure [S)). 




Oscillation: Am^ = 0.002; sin ^(28,3) = 0.5 
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Figure 8: Expected number of detected neutrino events as a function of baseline 
distance from the reactor core. The two curves show the expectations for the 
case of no oscillations and for an oscillation with Am^ = 0.002 and an amplitude 
of sin^(26'i3) = 0.5 which is 2.5 times the current limit. The curves are calculated 
using the complete energy spectrum from a nuclear reactor and a luminosity of 
eOOtGWy. 



In order to make the oscillation effect in Figure |H1 visible, the ratio of the two 
curves is shown in FigureEl Notice that the largest deviation occurs at a baseline 
distance of just over 2 kilometers. This corresponds to the first oscillation for 




Figure 9: The ratio of the expected number of neutrino events with and without 
oscillations as a function of distance from the reactor core. This calculation 
was made for a luminosity of 600 1 GW y and includes the true neutrino energy 
spectrum. The oscillation is assumed to have Ato^ — 0.002 and an amplitude 
of sin^(26'i3) = 0.5 which is 2.5 times the current allowed limit. 
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El, = 3.8 MeV as shown in Figure [7| This makes sense since this is the peak 
of the neutrino energy spectrum and therefore has the most statistical power. 
However, it is clear that as basehne distance increases, the effect of other parts of 
the energy spectrum being at their respective maxima and minima of osciUation 
effectively neutrahzes any ability to detect a specific oscillation signature. 

3.2 Spectral Shape Information 

From Equation JH)), it is clear that neutrinos of differing energies oscillate with 
different frequencies. Figure El shows that observable oscillation effects in the 
total number of neutrinos detected wash out with increasing baseline distance. 
But by looking at the specific energy distribution of the detected neutrinos, 
more information is available. In Figure [TUl two comparisons of the normalized 




Figure 10: Comparisons of the expected measured neutrino energy spectra at 
various baseline distances in the case of oscillations. The oscillations are as- 
sumed to have Am? = 0.002 and an amplitude of sin^ (26*13) = 0.5 which is 
2.5 times the current limit in order to magnify the effect. These plots show 
the expected statistics for a luminosity of 600 1 GW y at the specified baseline 
distance. 

energy spectra are shown. These plots show statistical errors only for 0.2 MeV 
bins and a luminosity of 600tGWy at each location. As with the previous 
plots, the amplitude of the oscillation has been magnified by a factor of 2.5 
(sin^ (26*13) = 0.5) and a mass difference of Ato^ = 0.002 has been used. In 
addition, an energy resolution of 10% /Ve has been assumed. The plot on 
the left compares the spectrum at 100 meters, which is effectively unoscillated, 
with the spectrum at 1.3 kilometers. One can notice that at this distance, the 
low energy part of the spectrum is showing a deficit with respect to the near 
detector spectrum. However, one could also confuse this with an overall shift in 
the absolute energy scale between the two detector locations. 

The right plot of Figure ^| however, compares two oscillated spectra from 
1.5 and 3 kilometers. Notice that the shapes of the spectra are vastly different. 
This arises from the fact that at 3 kilometers it is the high energy part of the 
spectrum which is fully oscillated away while the low energy part has returned 
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to full presence. It is interesting to note that the statistical significance of the 
spectral distortions in the two plots is nearly identical. While the total number 
of events is significantly higher in the plot with baselines of 100 m and 1.3 km, 
the spectral distortion is much more pronounced in the plot with baselines of 
1.5 and 3 km. These two effects appear to compensate one another. 

3.3 Combining Shape and Rate Information 

It turns out that the most statistically significant spectral shape distortion, 
given the assumed oscillation parameters above, is achieved for a near detector 
at the closest possible location and a far detector at about 1 kilometer. Thus the 
spectral shape information has a different optimal baseline than the depletion 
of the total fiux, which was previously observed to be maximal at just over 
2 kilometers from the source. 

Since the spectral shape measurement requires the use of a normalized energy 
spectrum at each location, the statistical significance of each measurement (each 
bin) is weighted by the total flux at that location. This gives a 1 /L^ reduction in 
the statistical sensitivity. Therefore, from a strictly statistical perspective, the 
total flux measurement will have slightly more than a factor of two more power. 
However, the total flux measurement is susceptible to systematic differences 
between the two detectors. Since the normalized energy spectra are normalized 
to the total number of events at that location, all systematic effects, except those 
which will be uncorrelated bin-to-bin within a detector, will be removed. This 
additional freedom from systematic effects implies that in the limit of infinite 
statistics, a more precise measurement can be made with the information from 
the energy spectrum. 

This interplay between the two methods can be seen graphically by refer- 
ring to the plot in Figure ITTI (see also the discussion in Sec. I3.4|l . There, the 
sensitivity to sin^(20i3) is shown as a function of luminosity. At low luminos- 
ity (£ < 400t GW y) the sensitivity is directly proportional to statistics and is 
dominated by the total flux measurement. Then, for luminosities between 400 
and 6000 tGWy, additional statistics do not make significant gains in the sen- 
sitivity. In this region, the systematic uncertainties between the two detectors 
(called (Jnorm in this plot) become dominant. However, beyond 6000 tGWy, 
notice that the sensitivity again becomes proportional to the statistics. This is 
caused by the fact that enough statistics have been gained to allow the spectral 
measurements to dominate over the systematically limited flux measurement. 

Realization of this interplay between the methods implies that the optimal 
choice of baseline distance depends on the expected luminosity of the experi- 
ment. Since most of the discussion in this paper does not expect a luminosity of 
greater than 6000 tGWy (which would require kiloton sized detectors), we will 
focus on measurements in which the total fiux measurement is not systematically 
limited. To make optimal use of the available statistics, we therefore wish to 
combine the information from the total fiux and energy spectral measurements 
at both detectors. One can rather simply create a chi-squared comparison of the 
two statistical distributions which takes into account both sets of information 
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Figure 11: The sensitivity to sin^ 26*13 at 90% CL as a function of the integrated 
luminosity for different values of the normahzation error cTnorm and the energy 
cahbration error (icai- The vertical lines mark the luminosities of the CHOOZ 
experiment and our standard setups Reactor-1 and Reactor-II. We take Arn^i = 
2 X 10"'^ eV^ and Lpu = 3 km. The horizontal line shows a typical sensitivity 
limit obtainable by the JHF-SK superbeam experiment (also known as JPARC- 
u) for the same parameter values. 



with the following definition: 




where L refers to the baseline distance to the near or far detectors and refers 
to the number of events in the i-th bin of the measured energy spectrum at that 
detector. By using the definition of Equation @, one can determine the 90% 
confidence level limit on a measurement of sin^(20i3) for a given luminosity. This 
is shown as a function of baseline distance in Figure 1121 assuming a luminosity 
of 600tGWy. For this estimation, the near detector is assumed to be fixed at 
300 meters and a 1% systematic limit has been used. 

One can see that for the current best fit value of Am^ — 0.002, the optimal 
baseline distance for the far detector is at 1.6 kilometers. However, given that 
the optimal distance depends quite strongly on the value of Am^ , we also show 
the sensitivity plots for values of Am^ which match the upper and lower bounds 
of the 90% allowed region from Super-Kamiokande. Notice that for each curve, 
the sensitivity to sin^(20i3) is relatively flat around the optimum. Therefore, a 
reasonable sensitivity can be reached, with flexibility to various values of Am^ , 
for a far detector baseline distance between 1.2 and 2.4 kilometers. 
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Figure 12: Statistical sensitivity to sin^ (26*13) as a function of baseline distance 
to the far detector. The statistical power is calculated using Equation © with 
a luminosity of GOOtGWy and a 1% systematic limit bin-to-bin. Curves are 
shown for three values of Am^ representing the best fit and the upper and lower 
limits of the 90% allowed region from the Super-Kamiokande experiment. 



3.4 Statistical Analysis and Luminosity Scaling 

In this section some general analysis methods are proposed to investigate the 
sensitivity to 9i3 in a reactor experiment with two detectors and one single 
reactor. The impact of the integrated luminosity, positions of the near detector, 
various systematic errors, and a possible background on the sensitivity limit 
are discussed. As a measure for the "size" of the experiment the integrated 
luminosity L is useful, which is defined as £ = fiducial detector mass [tons] x 
thermal reactor power [GW] x running time [years] (assuming 100% detection 
efficiency and no deadtimes). We define two benchmark setups Reactor-I (£ = 
400 t GW y) and Reactor-II (£ = 8000 t GW y) corresponding roughly to 31 500 
and 630 000 reactor neutrino events for no oscillations, respectively, assuming 
a PXE-based scintillator. 

We take into account that spectral information is available in the near, as well 
as in the far detector in form of iVbins bins in positron energy. For the theoretical 
prediction for the number of reactor neutrino events in the ith energy bin of the 
near (A — N) and far {A — F) detector, respectively, we write 

T^ = {l + a + b^ + c,)N,^+g^M,^, (10) 

and consider a x^"function including the full spectral information from both 
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detectors: 




Here, is the expected number of events in the «th energy bin of the cor- 
responding detector, which depends on the oscillation parameters. Nf' is cal- 
culated by folding the reactor neutrino spectrum, the detection cross section 
for inverse beta decay, the Ve survival probability, and the energy resolution 
function. In the numerical calculations we assume an energy resolution of 
5%/^Svis[MeV], and we use 62 bins in the range between Ep^ = 1.8 MeV and 
Eu^ = 8.0 MeV, corresponding to a bin width of 0.1 MeV.^ is the observed 
number of events. In the absence of real data, we take for Of the expected 
number of events for some fixed "true values" of the oscillation parameters, 
e.g., to calculate a sensitivity limit the expected number of events for 6'i3 = 
is used. If the near detector is so close to the reactor that no oscillations will 
occur, the T/^ will not depend on the oscillation parameters, and in that case 
one can set = Nf. The quantities Tf^ and Of correspond only to reactor 
neutrino events. If a certain background has to be subtracted from the actual 
number of events it will contribute to the statistical and systematic errors. In 
Equation (|ll|l Bf is the number of background events in the ith bin of detector 
A, and we assume that it is known with an (uncorrelated) error ctbg- 

For each point in the space of oscillation parameters, the x^-function has to 
be minimized with respect to the parameters a, , , , , and q modeling 
various systematic errors. 

1. The parameter a refers to the error on the overall normalization of 

the number of events common to near and far detectors, and (Tq is typically 
of the order of a few percent. The main source for such an error is the 
uncertainty of the neutrino flux normalization. 

2. The parameters b^ and parameterize the uncorrelated normaliza- 
tion uncertainties of the two detectors. Here contributes, for instance, 
the error on the fiducial mass of each detector. We assume that in this 
case an error below 1% can be reached. 

3. The energy scale uncertainty in the two detectors is taken into ac- 
count by the parameters g^ and g^ . To this aim we replace in the 
visible energy i?vis by (1 + g'^)i?vis, and to first order in g^ we have 
Nf{g^) w Nf{g^ = 0) + g^ M^. A typical value for this error on the 
energy calibration is deal ~ 0.5%. 

•^The results depend very weakly on our choices for the energy resolution and the number 
of bins. 
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4. In order to model an uncertainty on the shape of the expected en- 
ergy spectrum, we introduce a parameter q for each energy bin, known 
with an error (Jshapc- This corresponds to a completely uncorrelated error 
between different energy bins, which is the most pessimistic assumption 
of no knowledge of possible shape distortions. However, we choose this 
error fully correlated between the corresponding bins in the near and far 
detector (the same coefficient is used for the corresponding bins in the 
two detectors), since shape distortions should affect the signals in both 
detectors of equal technology in the same way. 

5. We include the possibility of an uncorrelated experimental system- 
atic error (Jcxp- In this way we assume that the observed number of 
events in each bin and each detector Of has in addition to the statistical 
error zLy/Of the (uncorrelated) systematic error zLac^pOf. We call this 
uncertainty "bin-to-bin error" . Taking it completely uncorrelated between 
energy bins as well as between the near and far detectors corresponds again 
to the worst case scenario. Values of CTexp at the per mil level should be 
realistic. 

Note that all the parameters describing the systematic errors are at most at 
the percent level, which means that the linear approximation in Equation H1U|I 
is justified. The following discussion of general features of such an analysis is 
based on the results reported in Reference |30| . 

Let us first assume that the near detector is close enough to the reactor, 
such that no oscillations develop {Lnd < 200 m). Furthermore, we first assume 
that the background is negligible, and Coxp ^ 0.1% and iTshape ^ 2%, which 
means that these errors can be neglected. Then the x^-analysis can be signifi- 
cantly simplified. In particular, it is not necessary to explicitly include the near 
detector in the analysis and Equation becomes |^ 

^ Cj V '^cal J \ CTnorm / 

with 111 

CTa 

For example, assuming typical values of CTq = 2% for the flux uncertainty and 
CTf, = 0.6% for the detector-specific uncertainty, we obtain with Equation H13|l 
an effective normalization error of (T,iorm — 0.8%. This is the value which is used 
for the numerical calculations. 

In Figure El we show the sensitivity to sin^ 26'i3 as a function of the inte- 
grated luminosity C In this figure, the lower diagonal curve corresponds to the 
idealized case of statistical errors only, and shows just the expected scal- 
ing, whereas the values (Xnorm = 0.8% and tJcai = 0.5% lead to the thick curve. 
At a luminosity around lOOtGWy, we detect a departure from the statistics 
dominated regime into a flatter systematics dominated region. This effect is 



26 



dominated by the error on the normaHzation (Tnormj whereas the energy cahbra- 
tion error acai only plays a minor role. In fact, we find that for all considered 
cases (including Reactor-I and Reactor-II) the impact of the energy scale un- 
certainty is very small, as long as the oscillation minimum is well inside the 
observable energy spectrum. Hence, we will neglect this error in the following. 

At large luminosities > lO^tGWy, the slope of the curve changes, and we 
are entering again a statistics dominated region with a \/\/L scaling. This 
interesting behavior can be understood as follows. For low luminosities the 
sensitivity to sin^ 2^13 comes mainly for the total number of events, hence the 
absolute normalization is important. The turnover of the sensitivity line into 
the second statistics dominated region occurs at the point where the spectral 
shape distortion becomes more important than the total rate, which implies 
that the overall normalization, and consequently also the actual value of tJnorm, 
becomes irrelevant. We illustrate this by the upper thin black line, which shows 
the luminosity scaling for the case of larger systematic errors. As an example, 
we choose values of CTnorm = 1.7% for the normalization and Ucai = 1.1% for the 
energy calibration. We find that, in this case, the transition to the systematics 
dominated regime occurs at much smaller luminosities. However, for large lu- 
minosities, the same limit is approached as for the more optimistic case. The 
diagonal gray curve shows the limit for no constraint at all on the normalization 
and energy calibration.^ Even in this extreme case, we obtain the same limit 
for high luminosities. 

3.5 Systematics, Background, and the Position of the Near 
Detector 

In the following we investigate the impact of various systematic effects beyond 
the simple overall normalization. To this end we apply the full function as 
given in Equation Hll|l . In FigurelT^the luminosity scaling of the sin^ 26'13-limit 
is shown for various choices for the experimental bin-to-bin error dexp (see itemEl 
above) and background levels in the far detector. For the sake of concreteness 
we assume a flat background in each detector Bf — with an error of ctbg = 
10%. The size of the backround is measured by /bg, which is defined as the 
fraction of the total number of background events relative to the total number 
of reactor neutrino events for no oscillations, i.e., B^ = /bg X^i -^z^nooscZ-^bins- 
From the figure we find that Reactor-I is not afi^ected by a bin-to-bin error up 
to 0.5%, nor by backgrounds in the far detector up to 5%. In contrast, such 
errors and backgrounds are important to some extent for big experiments like 
Reactor-II. In that case values of cTexp > 0.1% start to deteriorate the sensitivity 
limit, and the background in the far detector should be smaller than 1% of the 
reactor neutrino signal. We note that backgrounds in the near detector up to 
a few percent do not affect the result. Regarding the huge number of reactor 
neutrino events in the near detector it should be possible to obtain backgrounds 
below 1%. 

^Although we leave the normalization free in the fit, we assume that the shape is known. 
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Figure 13: Sensitivity to sin^ 29i3 at 90% CL as a function of the luminosity for 
an uncorrelated experimental systematic error ("bin-to-bin error") dcxp = 0.1% 
and 0.5%, and background levels in the far detector relative to the total number 
of events for no oscillations of 1% and 5%. Here L^d = 0.2 km, Lpu = 1.7 km, 
Am|]^ = 2 X 10~^ eV'^, and dshapc = 2%. Identical detector masses are assumed 
for near and far detectors. 
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Figure 14: The sensitivity to sin^ 26*13 at the 90% CL for Reactor-I and Reactor- 
II as a function of the near detector position. The far detector is situated at 
1.7 km and we assume identical detector masses and Amgj^ = 3 x 10"'^ eV^. Fur- 
thermore, the impact of an uncorrelated theoretical shape uncertainty cTshape = 
2% is shown. 



For practical reasons it might be hard to find a reactor station where a near 
detector can be situated very close (< 200m) to the core with sufhcient rock 
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Reactor-I Reactor-II 


Effective normalization cTnorm = 0.8% 
Energy calibration deal — 0.5% 
Exp. bin-to-bin uncorr. error doxp ^ 0.1% 
background in far detector /bg ^1% 
near detector baseline 

theor. shape uncertainty cTshapc = 2% 


important not important 
not important not important 
not important important 
not important important 
Lnd ^ 400 m Lisiu ^ 1 km 

important for ?J 1 km 



Table 2: Relevance of various systematic errors for the two reactor benchmark 
setups Reactor-I and Reactor-II. 

overburden. Therefore, it is interesting to investigate the impact of larger near 
detector baselines on the sin^ 2(?i3 limit. In this case the information provided 
by the near detector on the initial flux normalization and energy shape is al- 
ready mixed with some oscillation signature. The results of such an analysis are 
presented in Figure We find that for the case of Reactor-I the limit starts 
deteriorating around a near detector distance of 400 m, whereas for Rcactor-II 
the limit even improves slightly up to near detector baselines of ^ 1 km. Due 
to the high statistics in the case of Reactor-II, flux normalization and shape 
are very well determined by the near detector even in the presence of some ef- 
fect of ^13, and the additional information on oscillations improves the limit a 
bit. Furthermore, we find from Figure [T^ that the shape uncertainty becomes 
important for near detector baselines > 1.1km, especially for Reactor-II. A re- 
duction of this theoretical error would be helpful in such a situation. We note 
that assuming the shape error to be completely uncorrelated corresponds to the 
worst case. A more realistic implementation of the shape uncertainty including 
correct correlations will lead to results somewhere in between the curves for no 
and 2% shape error in Figure We have verified that for Reactor-I the wors- 
ening of the limit comes mainly from the fact that with increasing near detector 
baselines the number of events decreases rapidly, i.e., it is statistics dominated, 
whereas for Reactor-II the loss in sensitivity is driven by the systematic shape 
uncertainty and cannot be compensated by larger near detectors. 

To summarize, for the case of the Reactor-I setup (400 t GW y) the main 
information comes from the total number of events and the systematic normal- 
ization error dominates. In order to obtain a reliable limit, it should therefore 
be well under control. For large luminosities, such as for the Reactor-II setup 
(8000 t GW y), the sensitivity limit comes mainly from spectral shape infor- 
mation and is independent of normalization errors. In that case a bin-to-bin 
uncorrelated experimental systematic error should be below 0.1% and the back- 
ground should be at the 1% level. Furthermore, for the case of Reactor-I-like 
experiments one should look for a site where the near detector can be placed 
at a distance of at most 400 m from the reactor. For large detectors, such as 
Reactor-II, near detector baselines of up to 1 km will perform well. For near 
detector baselines longer than about 1 km the correct treatment of the theo- 
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retical shape uncertainty becomes important. These results are summarized in 
Table 12 
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4 Previous Reactor Experiments 



Here we review the three raost recent previous neutrino experiments at reac- 
tors. The best current Umit on 613 conies from CHOOZ. An experiment with 
similar distance and running time, but smaUcr overburden, was conducted at 
Palo Verde. Finally the KamLAND experiment was conducted with much larger 
overburden, and larger size, but much longer baseline, which resulted in a re- 
duced sensitivity to ^13 but unprecedented sensitivity to 6*12. These previous 
experiments are reviewed in order to present several lessons which are needed 
to show that a future experiment can control systematics to the level below one 
percent. 

4.1 CHOOZ 

The CHOOZ experiment |1(J[ IllL 1^ was located close to the Chooz nuclear 
power plant, in the North of France, 10 km from the Belgian border. The power 
plant consists of two twin pressurized water reactors (PWR), the first of a series 
of the newly developed PWR generation in France. The thermal power of each 
reactor is 4.25 GW (1.3 GW electrical). These reactors started respectively in 
May and August 1997, just after the start of the data taking of the CHOOZ 
detector (April 1997). This opportunity allowed a measurement of the reactor- 
OFF background and a separation of individual reactor's contributions. 

The detector was located in an underground laboratory about 1 km from the 
neutrino source. The 300 mwe rock overburden reduced the external cosmic ray 
muon flux, by a factor about 300, to a value of 0.4 m~^ s~^. This was the main 
criterion to choose this site. Indeed, the previous experiment at the Bugey reac- 
tor power plant showed the necessity of reducing by 2 orders of magnitude the 
flux of fast neutrons produced by muon-induced nuclear spallations in the ma- 
terial surrounding the detector. The neutron flux in CHOOZ was measured at 
energies larger than 8 MeV (endpoint of the neutrino flux from nuclear reactors) 
and found to be 1/day, in good agreement with expectation. 

The detector envelope consisted of a cylindrical steel vessel, 5.5 m in diameter 
and 5.5m in height. The vessel was placed in a pit (7m diameter and 7m deep), 
and was surrounded by 75 cm of low activity sand. It was composed of three 
concentric regions: 

• a central 5 ton target in a transparent plexiglass container filled with a 
0.09 % Gd-loaded scintillator 

• an intermediate 70 cm thick region, filled with non loaded scintillator and 
used to protect the target from PMT radioactivity and to contain the 
gamma from neutron capture on Gd. These 2 first regions were viewed by 
192 PMTs 

• an outer veto, filled with the same scintillator. 

The scintillator showed a degradation of the transparency over time, which 
resulted in a decrease of the light yield (live time around 250 days). The event 
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position was reconstructed by fitting the charge balance, with a typical precision 
of 10 cm for the positron and 20 cm for the neutron. The time reconstruction 
was found to be less precise on source and laser tests, due to the small size of the 
detector. The reconstruction became more difficult when the event was located 
near the PMTs, due to the divergence of the light collected (see Figure 31 

of cni). 

The final event selection used the following cuts : 

• positron energy smaller than 8 MeV (only 0.05% of the positrons have a 
bigger energy) 

• neutron energy between 6 and 12 MeV 

• distance from the PMT surface bigger than 30 cm for both positron and 
neutron 

• distance between positron and neutron smaller than 100 cm 

• only one neutron 

• time window between positron and neutron signals is from 2 to 100 fis. 

The 6 MeV cut on the gamma ray's total energy from a neutron capture 
on Gd cannot be computed from a simulation, because only the global released 
energy is known. The number of gammas and their individual energies were 
very poorly known. The scintillating buffer around the target was important to 
reduce the gamma escape. This cut was calibrated with a neutron source (0.4% 
systematic error). The 3 cuts on the distances were rather difficult to calibrate, 
due to the difficulty of the reconstruction described above. This created a 
tail of mis-reconstructed events, which was very difficult to simulate (0.4% 
systematic error on the positron- neutron distance cut). The positron threshold 
was carefully calibrated, as shown in Figure 39 of TU'. The value of the threshold 
depends upon the position of the event, due to the variation of solid angles and 
to the shadow of some mechanical pieces such as the neck of the detector (0.8 % 
systematic error). The time cut relied on MC simulation. The time spectrum 
happened to be exponential to > 20 ^s, but there was no reason for this (the 
Gamow law, which allows to demonstrate an exponential behavior is wrong 
for Gd, whose capture cross section is only epithermal). The corresponding 
systematic error was estimated to be 0.4%. 

The final result was given as the ratio of the number of measured events 
versus the number of expected events, averaged on the energy spectrum. It 
was: 

R = 1.01 ± 2.8% (stat) ± 2.7% (sys) 
Two components were identified in the background : 

• a correlated one, which has a flat distribution for energies bigger than 
8 MeV, and is due to the recoil protons from fast spallation neutrons. It 
was extrapolated to 1 event/day. 
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• an accidental one, which is obtained from the measure of the single rates. 

The background was measured while the reactor was off, and by extrapolating 
the signal versus power straight line (see Figure 49 of Q^ii). It is in good agree- 
ment with the sum of the correlated and accidental components measured as 
1.41 ± 0.24 events per day. These numbers have to be compared to a signal of 
26 events /day at full power. 

The systematic errors were due mostly to the reactor uncertainties (2%), 
to the detector efficiency (1.5%), and to the normalization of the detector, 
dominated by the error on the proton number from the H/C ratio in the liq- 
uid (0.8%). The resulting exclusion plot is shown in Figure 58 of ^01- The 
corresponding limit on sin^ 26* is 0.14 for Am^ — 2.6 x lO^'^eV^, and 0.2 for 
Am^ = 2.0 X IQ-^eV^. Due to specific source-detector distance of about 1 km, 
no limit on sin^(26') can be set for Am^ = 0.8 x 10~'^eV^, due to the limited 
distance from the between the cores and the CHOOZ detector. 
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4.2 Palo Verde 



The Palo Verde experiment was motivated by the discovery of the atmospheric 
neutrino anomaly |7D1 1^ [75] , which could be explained by neutrino oscilla- 
tions with large mixing angle and a mass-squared difference in the range of 
10-3 - 10-2 eV^. The Palo Verde experiment, together with the CHOOZ ex- 
periment ^1 nil IH^ with a similar baseline, were able to exclude the Vg 
oscillations as the dominant mechanism for the atmospheric neutrino anomaly. 
While Palo Verde pursued its goal of exploring the then unknown region of 
small Am^, results from Super-Kamiokande [73] were published which favored 
the — > Vr oscillation channel over —f Vg. In this section, we provide a 
brief description of the Palo Verde experiment and its final results. Details on 
the experiment can be found in its physics publications I12L I75L [75] and in 
technical publications cited below. 

The Palo Verde experiment was carried out at the Palo Verde Nuclear Gen- 
erating Station, located about 80 km west of Phoenix, Arizona. The largest 
nuclear power plant in the Americas, Palo Verde consists of three identical pres- 
surized water reactors with a total thermal power of 11.63 GW. The detector, 
containing 11.3 tons of liquid scintillator for the neutrino target, was located 
at a shallow underground site, 890 m from two of the reactors and 750 m from 
the third. The 32 meter-water-equivalent overburden entirely eliminated any 
hadronic component of cosmic radiation and reduced the cosmic muon flux. 
The collaborating institutions on the experiment were the California Institute 
of Technology, Stanford University, University of Alabama, and Arizona State 
University. Data were collected from 1998 to 2000. 

A schematic view of the detector is shown in FigureEl The central detector 
was an 11 x 6 matrix of cells. Each cell was 9 m long, subdivided into a 740-cm 
central section filled with Gd-loaded liquid scintillator (TTj and an 80-cm section 
of mineral oil at either end. The cell was viewed at each end by a 5-inch PMT. 
Surrounding the central detector along the long sides were tanks providing a 
layer of water shielding 1 m thick. The water and mineral oil shielding sections 
attenuate gammas and neutrons emitted from the laboratory walls and outer 
components of the detector, e.g. the glass of the PMTs. The detector was fully 
enclosed by liquid scintillator detectors used to veto cosmic muons. 

Electron antineutrinos were detected via inverse beta decay, manifested ex- 
perimentally as a prompt energy deposit due to the kinetic energy and annihi- 
lation energy of the positron followed an average of ~28 fis later by a gamma 
cascade of 8 MeV total energy due to capture of the neutron on Gd. The cen- 
tral detector was segmented in order to improve the discrimination between 
positrons from inverse beta decay and electrons, gammas and recoil protons. 
The experimental signature required for a positron was an energy deposit in one 
cell greater than '^l MeV (kinetic energy of the positron) and energy deposits 
in adjacent cells consistent with those expected from back-to-back 511 keV 
annihilation gammas. 

The event trigger was based on a so-called triple. For triggering, the anode 
output of the PMT was split and sent to two sets of discriminators, one set 
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Figure 15: Schematic view of the detector and the inverse /3-decay reaction 
producing a triple coincidence pattern inside the detector. 



having a threshold corresponding to an energy deposit of ~50 keV ( "LO" ) in the 
cell and the other set having a threshold corresponding to ^--^500 keV ("HI"). The 
discriminator outputs were fed into a fast trigger processor |ZHI which generated 
a triple if there was a coincidence between at least 2 LOs and 1 HI in any 5x3 
cell submatrix in the detector. The occurrence of a triple initiated digitization 
of the associated event. Readout was carried out if two triples occur within 450 
IIS of each other. Given the proximity in time of the "prompt" and "delayed" 
part of a candidate event, two banks of Fastbus ADCs and TDCs had to be 
used for digitization. The trigger rates for triples and correlated triples were 
approximately 50 Hz and 1 Hz, respectively. 

The muon veto hit rate was about 2 kHz. A hit in the veto generated 5 fis 
of deadtime for the triple trigger processor. Otherwise, muon hits were only 
clocked and latched for readout, and the main fj, veto cuts were applied off-line. 

Detector calibration for energy and position reconstruction was carried out 
using 7 point sources, blue LED's, and a fiber optic flasher system. The detector 
simulation program used to estimate the triple trigger efficiencies was tuned and 
checked against data taken with ^^a and ^^e sources for the case of positrons 
and with a ^^^f source and a tagged Am-Be source for the case of neutron 
capture. Detector stability between calibrations was monitored using the LED's 
and fiberoptic flasher system. 

The expected flux was calculated from the reactor power and fuel compo- 
sition. The expected Dc interaction rate in the whole target, both scintillator and 
the acrylic cells, is plotted in Figure [TBI for the case of no oscillation from July 
1998 to July 2000. Around 220 interactions per day are expected with all three 
units at full power. Four periods of sharply reduced rate occurred when one of 
the three reactors was off for refueling, the more distant reactors contributing 
each approximately 30% of the rate and the closer reactor the remaining 40%. 
The short spikes of decreased rate are due to accidental reactor outages, usually 
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less than a day. The gradual decline in rate between refuelings is caused by fuel 
burn-up, which changes the fuel composition in the core and the relative fission 
rates of the isotopes, thereby affecting slightly the yield and spectral shape of 
the emitted i/e flux. 
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Figure 16: The calculated interaction rate in the detector target for the case 
of no oscillations. The four long periods of reduced flux from reactor refuelings 
were used for background subtraction. The decreasing rate during the full power 
operation is a result of the changing core composition as the reactor fuel is 
burned. 

Inverse beta decay candidates were selected according to the following cri- 
teria. Each subevent (prompt and delayed) had at least one hit with energy 
greater than 1 MeV and at least two additional hits with energy greater than 
30 keV. The energy thresholds of this cut were chosen to select events in the 
energy ranges where the triggers were efficient. Any event with hits greater than 
8 MeV in either subevent was discarded. The magnitude and pattern of energy 
deposits in the prompt subevent were required to resemble what was expected 
from the kinetic energy of the positron and its annihilation. The prompt and 
delayed subevents of the event were required to be correlated in space and time. 
To further suppress backgrounds, an event was accepted if it started at least 
150 /is after the last veto hit and at least 3.5 MeV of energy was deposited 
in either the prompt or delayed subevent. For the case of no oscillations, the 
energy-dependent combined efficiency of the trigger and selection cuts on neu- 
trino interactions is about 18%. The deadtime induced by the veto-dependent 
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hardware and software cuts further reduced the efficiency to about 11%. The 
event rate of ~55 day"'^ after selection may be compared to an expected signal 
rate of about 20 day~^ for no oscillations. Below, "positron cuts" refer to cuts 
applied to the prompt subevent and "neutron capture cuts" to the cuts applied 
to the delayed subevent. 

Backgrounds surviving the event selection may be naturally classified as 
uncorrelated and correlated. Uncorrelated backgrounds are due to random co- 
incidences between triple triggers within the delayed coincidence window. The 
dominant source of uncorrelated events is natural radioactivity. Correlated back- 
ground events are events in which both subevents are due to the same process. 
The main source of this type of background are neutrons from muon spallation 
or capture. These events are mainly comprised of proton-neutron events-in 
which a single neutron deposits its kinetic energy by scattering from protons 
and is then captured-and double neutron events-in which two (typically ther- 
mal) neutrons from the same spallation event are captured in the detector. The 
interevent time distribution for uncorrelated background events followed an ex- 
ponential function with a time constant of 500 ^s, as would be expected given 
the muon veto rate of ^2 kHz and the veto-dependent event selection require- 
ments. This time dependence is slow compared to that of signal and correlated 
backgrounds, hence the contribution of the uncorrelated background was iso- 
lated and studied by looking at long interevent times. Based on these studies, 
the contribution of uncorrelated backgrounds to the event rate after selection 
was estimated to be about 7 day~-^. 

To estimate the contribution of the correlated backgrounds, two different 
approaches were used. In the first approach, the so-called "reactor power" 
method, the correlation between reactor power and observed event rate was an- 
alyzed: the signal rate would vary with reactor power while the background rate 
is independent of reactor power. In the second approach, the so-called "swap" 
method |7H], the cuts on the prompt and delayed subevent were interchanged 
and the resulting event rate was subtracted from the event rate obtained with 
the standard cuts. This approach efficiently removed backgrounds such as un- 
correlated events and double neutrons, which are symmetric in the prompt and 
delayed subevents while keeping most (~80%) of the neutrino signal. The re- 
maining important source of background, namely proton-neutron events, was 
estimated from simulation. Owing to the fact that only one reactor was refu- 
eled at a time and the refueling time was short (approximately 30 days every 6 
months), the swap analysis had more statistical power than the reactor power 
analysis, but the swap analysis had an additional contribution to the systematic 
error from the uncertainty in the background. 

Table |3| summarizes the observed and corrected rates for the main data 
taking periods of the experiment. Data was accumulated at full reactor power 
for four periods and one reactor out of three was off for four periods. 

Figure IT7I shows the energy spectrum of the neutrino candidates after back- 
ground subtraction using the reactor power method. Also shown are the ex- 
pected spectra for no oscillations and for oscillations based on the Kamiokande 
best fit (assuming it is due to i/^ Ve transitions). The observed spectrum is 
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Table 3: Data taking periods, efficiencies (including livetime), measured event 
rates A''i and N2 [Ni - event rate after applying the neutrino selection cuts, N2 
- event rate obtained by applying the "swapped" selection cuts), ABpn - the 
residual contribution to {N1 — N2), mainly due to the proton-neutron component 
of the correlated background, and estimates of the background. N^, and Ncaic 
are corrected measured neutrino event rate and calculated expected rate for no 
oscillations. Uncertainties are statistical only. 



Period 


1998 


1999-1 


Reactor 


on 


890m off 


on 


750m off 


time (days) 


30.4 


29.4 


68.2 


21.8 


efficiency (%) 


8.0 


8.0 


11.5 


11.6 


Ni (day-^ 


39.6 ± 1.1 


34.8 ± 1.1 


54.9 ± 0.9 


45.1 ± 1.4 


.¥2 (dav^^ 


25.1 ± 0.9 


21.8 ± 0.9 


33.4 ± 0.7 


32.0 ± 1.2 


AZJp, (day 


0.88 


0.89 


1.11 


1.11 






efficiency con 


ected results 




Background 


292 ± 11 


255 ± 10 


265 ± 6 


266 ± 10 




202 ± 19 


182 ± 18 


212 ± 10 


124 ± 17 


^calc 


216 


154 


218 


129 


Period 


1999-11 


2000 


Reactor 


on 


890m off 


on 


890m off 


time (days) 


60.4 


29.6 


83.2 


27.5 


efficiency (%) 


11.6 


11.6 


10.9 


10.8 


A^i (day-^ 


54.2 ± 0.9 


49.4 ± 1.3 


52.9 ± 0.8 


43.1 ± 1.3 


N2 (day-i 


32.5 ± 0.7 


32.6 ± 1.0 


30.2 ± 0.6 


30.4 ± 1.1 


ABpn (day-i) 


1.11 


1.11 


1.07 


1.07 






efficiency corrected results 




Background 


256 ± 6 


265 ± 9 


249 ± 5 


272 ± 9 




214 ± 11 


161 ± 15 


237 ± 10 


129 ± 16 




220 


155 


218 


154 
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consistent in shape and normalization with the hypothesis of no osciUations. 
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Figure 17: The prompt energy spectrum after oti- Ojff subtraction averaged over 
the 4 pairs of on/ Ojff periods. The histograms show the corresponding expecta- 
tions for no oscillations (solid line) and the Kamiokande best fit (dashed line). 

Estimates of the systematic uncertainties in the reactor power and swap 
analyses are presented in Table 01 The systematic uncertainty received contri- 
butions from the detection efhciency and the flux calculation. In addition, the 
reactor power method suffered a systematic error from background variations, 
and the systematic uncertainty in the swap method had a contribution from the 
uncertainty in the estimate of ABp^. 

A analysis, using the Feldman-Cousins prescription |8()| for determining 
the 90% CL acceptance region and taking into account statistical and systematic 
uncertainties, was performed on the data to determine the regions in Am^ — 
sin^(26') space excluded at the 90% CL. The analysis was performed for the 
reactor power method and for the swap method. The results are shown in 
Figure [El 
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Table 4: Contributions to the systematic error of the reactor power and swap 
analysis. 



Error source 


on-ofr(%) 


swap(%) 


e+ efficiency 


2.0 


2.0 


n efficiency 


2.1 


2.1 


Pe flux prediction 


2.1 


2.1 


Pe selection cuts 


4.5 


2.1 


Bp„ estimate 


N/A 


3.3 


Background variation 


2.1 


N/A 


Total 


6.1 


5.3 




Figure 18: Regions of Am^ — sin^ 26 plane (two flavor oscillations) excluded 
at the 90% CL by the reactor power analysis (dashed curve) and swap analysis 
(solid curve). Also shown is the Kamiokande allowed region and best fit and 
the region excluded by the CHOOZ experiment. 
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4.3 KamLAND 



The KamLAND experiment, located in the underground Kamioka laboratory 
was conceived to test, in a wholly terrestrial experiment, the MSW LMA neu- 
trino oscillation solution to the solar neutrino problem by searching for disap- 
pearance of VeS emitted by Japanese nuclear power plants. The power-weighted 
average distance between Kamioka and Japanese nuclear power plants is approx- 
imately 175 km. The total reactor power is such that, for no disappearance, one 
would expect to see about 1.5 inverse beta decay events per day at Kamioka 
in a 1 kton liquid scintillator detector. Assuming two-flavor oscillations with 
mass-squared difference Am^ and mixing angle 0, the probability that a Ve of 
energy E survives over a distance L is given by 

, , o 1.27 Arri^L 

P = 1 - sin^ 29 sin^ 

E 

for Arn^ in eV^, L in m, and E in MeV. Given that the mean cross-section- 
weighted neutrino energy is about 5 MeV, KamLAND would thus be sensitive 
to values of Am^ below 10~^ eV^ for large mixing angle, reasonable running 
time, and well-controlled backgrounds. 

The KamLAND experiment is being carried out by a collaboration of uni- 
versities and laboratories from the United States and Japan. The international 
character of the collaboration has been essential for realizing an experiment 
which has the the advantages of an excellent site, sophisticated detector tech- 
nologies, and a team of physicists with experience and expertise in mounting 
reactor neutrino experiments. 

In this section we briefly describe the KamLAND experiment and summarize 
its first published results j7| . While the design of the experiment allows it to 
measure other quantities, e.g. the ^Be component of the solar neutrino flux, the 
focus of this section is exclusively on the reactor neutrino measurement. Further 
details on the detector and data analysis can be found in 7 . 

A schematic view of the KamLAND detector is shown in Figure ^1 The 
experiment target consists of 1 kton of ultra pure liquid scintillator [20% pseu- 
documene, 80% mineral oil (dodecane), and 1.52 g/1 PPO] contained in a 13-m- 
diameter transparent balloon. The target is viewed by 1,879 PMTs, of which 
1,325 are 17-inch and the balance 20-inch, providing a photocathode cover- 
age of about 34%. The PMTs are bolted to a 18-m-diameter stainless steel 
sphere (containment vessel). The volume between the balloon and stainless steel 
sphere is filled with mineral oil, buffering the target against natural radioactiv- 
ity in the PMTs, stainless steel sphere, and surrounding rock and against fast 
neutrons generated by muon spallation outside the containment vessel. A UV- 
transparent acrylic sheet 3 mm thick mounted just in front of the PMT faces 
acts as a barrier to radon emanating from the PMTs. The stainless steel sphere 
is enclosed by the outer detector, a water Cerenkov muon detector which con- 
sists of a cylindrical tank of pure water viewed by 225 PMTs. The detector is 
located in the Kamioka Underground Observatory, in the cavern formerly occu- 
pied by the Kamiokande Experiment. The rock overburden exceeds 2,700 mwe. 
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Figure 19: Schematic view of the KamLAND detector 

resulting in a muon rate through the experiment target of about 0.34 Hz. 

There are 16 commercial nuclear power plants in Japan, accounting for 97% 
of the neutrino flux at KamLAND. The bulk of the flux-about 80%-is due to 
reactors 138-214 km away. The mean Dp energy, weighted by the inverse beta 
decay cross section, is about 5 MeV. Records on thermal power, burn-up, and 
fuel exchange are furnished to the experiment on a continuing basis by the 
plant operators. The total reactor power varies by 20-30% throughout the year 
(down to 50% in 2003), which provides a means of subtracting backgrounds by 
investigating the correlation of the event rate with reactor power. 

The primary trigger for the experiment presently requires 200 PMT hits, 
corresponding to an energy threshold of about 0.7 MeV. Following each primary 
trigger, the threshold is lowered to 120 PMT hits for 1 ms to detect low energy 
delayed activity. In addition to being discriminated for the trigger, the PMT 
voltages arc sampled and digitized by waveform analyzers (ATWDs). There are 
two ATWDs per PMT, allowing two-step sequential events to be fully recorded. 
The ATWD sampling rate is about 630 MHz, and 128 samples are acquired per 
waveform. Signals in the outer detector are recorded as part of the standard 
data stream for offline analysis. The trigger rate is about 20 Hz. The amount 
of data recorded per day is about 150 GB. 

Energy estimation, vertex reconstruction, and detection efficiency are cal- 
ibrated using gamma sources ranging in energy from 0.279 MeV to 7.7 MeV, 
neutron sources (Am-Be), and light flashers (LEDs and lasers). Sources are 
deployed along the vertical axis of the experiment by winch operated from a 
sealed glove box at the top of the detector; in addition, there are blue LEDs 
permanently mounted on the stainless steel sphere. Besides detailed calibrations 
carried out from time to time, the detector is monitored on a weekly basis with 
gamma sources. Natural sources, namely spallation neutrons, cosmogenics, and 
natural radioactivity, are also used for detector calibration. 
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The signal process is inverse beta decay 

i>e + p — > + n. 

The experimental signature is a prompt energy deposit of 1-8 MeV, due to the 
positron kinetic energy and annihilation, followed an average of 200 fis later by 
emission of a 2.2 MeV gamma from neutron capture on hydrogen. Exploiting 
the delayed coincidence is key to controlling backgrounds. 

Following the first publication Ij, we now describe event selection, back- 
ground estimation, experimental uncertainties, and interpretation of data in 
the context of Df, disappearance and neutrino oscillations. 

The cuts applied to select inverse beta decay candidates are the following: 

1. Both the prompt and delayed subcvent vertices lie within 5 m of the 
detector center. 

2. The delayed subevent occurs within 0.5-660 fis after prompt subevent. 

3. The prompt and delayed subevent vertices are separated by less than 
1.6 m. 

4. The delayed subevent energy lies between 1.8 and 2.6 MeV. 

5. The delayed subevent vertex lies more than 1.2 m from the central vertical 
axis. 

6. The event occurs at least 2 s after a showering muon (energy deposit 
greater than ^^3 GeV) and at least 2 s after, or more than 3 m away from, 
any other muon track. 

7. The prompt subevent energy is greater than 2.6 MeV. 

Cut 1 is applied to control backgrounds due to natural radioactivity in the 
balloon system, the PMTs, and beyond and backgrounds due to muon spallation 
in the surrounding rock. Cut 5 controls backgrounds from natural radioactiv- 
ity in the thermometers suspended in the central detector. Cut 6 is made to 
suppress backgrounds from cosmogenics, e.g. ^He and ^Li, and spallation due 
to muons passing through the central detector. We have included Cut 7 to 
eliminate the low-energy region in which terrestrial radioactivity is expected to 
contribute. 

The efficiency of the Cuts 2-5 plus Cut 1 applied to the delayed subevent 
on inverse beta decay events is 78.3±1.6%. The effect of Cut 6 is taken into 
account in the livetime calculation. 

Applied to the data set acquired March 4 - October 6, 2002, which includes 
in all 145.1 d of live time, 54 events are selected. Figure EUl shows the energies of 
the delayed versus prompt subevents after all cuts have been applied but Cuts 
4 and 7. 

The contribution to the event sample from accidental delayed coincidences 
has been estimated by repeating the event selection with an off-time delayed 
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Figure 20: The energies of the delayed subevents in KaniLAND versus the 
energies of the respective prompt subevents for events passing all selection cuts 
but Cuts 4 and 7 (see text). 

coincidence window. The contribution is found to be small, namely 0.0086 ± 
0.0005 events. 

The residual cosmogenic background in the candidate sample has been esti- 
mated by analyzing the time and spatial correlations between muon tracks and 
event candidates. The result of this analysis is that the expected cosmogenic 
contribution to the sample is 0.94 ± 0.85 events. 

Fast neutrons from muon spallation can easily mimic inverse beta decay 
events and thus are a potentially dangerous background. Their contribution is, 
however, suppressed by the fiducial volume cut (r < 5 m). The contribution 
of this background has been estimated from a sample of muon events-where 
the muon misses the inner detector-which contain an inner detector event that 
passes the event selection cuts. The vertex distribution for this sample is ex- 
trapolated into the fiducial volume and then normalized to the muon track 
reconstruction inefficiency. This approach yields the result that less than 0.5 
events are due to fast neutrons. 

In summary, the total expected background in the event sample is 0.95±0.99 
events. 

For comparing the observed events with the expected number of events under 
different scenarios (no oscillations, two-flavor neutrino oscillations, etc.), a num- 
ber of quantities enter which carry systematic uncertainties. These quantities 
and the estimated magnitudes of the systematic uncertainties are: total liquid 
scintillator mass-2.1%; fraction of mass within fiducial volume-4.1%; energy 
scale at 2.6 MeV-2.1%; selection cut efficiency-2.1%; experiment live time- 
0.07%; reactor power-2.0%; fuel composition-1.0%; finite lifetime of fission 
products-0.3%; neutrino spectra~2.5%; and inverse beta decay cross section- 
0.2%. 
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Adding the individual contributions in quadrature, the total systematic un- 
certainty is estimated to be 6.4%. 

The integrated number of events expected for no disappearance is 86.8 ± 
5.6 events. Figure [771 shows, including accidental background and De's from 
terrestrial radioactivity, both (a) the expected prompt energy distribution for 
the case of no disappearance and (b) the corresponding distribution for the 
observed events. 
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Figure 21: Upper panel: Prompt KamLAND energy distributions for the ex- 
pected (no disappearance) case. Lower panel: same as upper panel except 
binned in energy and including energy spectrum of observed events. Also shown 
is the result of a fit to the observed events above 2.6 MeV in terms of two-flavor 
neutrino oscillations; the shaded area superimposed on this fit indicates the 
systematic uncertainty in the best-fit result. 

Accounting for the ^1 background event, the 54 events observed above 
2.6 MeV are inconsistent with the hypothesis of no disappearance at the 99.95% 
confidence level. 

One may instead interpret the observed events in terms of two-flavor neu- 
trino oscillations. Figure EH shows the result for analysis of the total rate and of 
the rate and shape of the energy spectrum. The best fit to the rate and shape 
gives Am^ = 6.9x10"^ eV^ and sin" 26* = 1.0. Assuming CPT invariance, it 
can be seen that the KamLAND observation is consistent with the MSW LMA 
solution to the solar neutrino problem and further restricts the allowed LMA 
region. 

In summary, since starting to take data in January 2002, KamLAND has ob- 
served reactor anti-neutrino disappearance with a very high level of confidence 
(99.95%). The MSW LMA solution is the only neutrino oscillation solution to 
the solar neutrino problem which is consistent with the KamLAND result and 
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Figure 22: Excluded regions of neutrino oscillation parameters for the rate 
analysis and allowed regions for the combined rate and shape analysis from 
KamLAND at 95% CL At the top are the 95% CL excluded region from CHOOZ 
and Palo Verde experiments, respectively. The 95% CL allowed region of the 
LMA solution of solar neutrino experiments is also shown. The solid circle shows 
the best fit to the KamLAND data in the physical region 
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CPT invariancc. 

KamLAND continues to take data and make improvements to the detector 
and data analysis. These improvements include inclusion of the 20-inch tubes in 
the data analysis, installation of a An deployment system for calibration sources, 
deployment of a /3 source and additional gamma sources for better understanding 
of the detector energy response, and development of a full detector simulation. 
These improvements arc expected to reduce significantly the uncertainties on 
the amount of mass contained within the fiducial volume, the energy scale at 
2.6 MeV, and the selection cut efficiency. Improvement in the accuracy to which 
the reactor related quantities, e.g. reactor power, arc known is also anticipated. 
Achieving a total systematic error in the range of 3-4% appears realistic, in 
which case the dominant source of error for the experiment after several years 
of data taking will be statistical. 
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5 Detector Design 



From the discovery of the neutrinos by Reines and Cowan [SI] at Savannah 
River to the evidence for disappearance at KamLAND '7', reactor neutrino 
experiments have used the same fundamental design. A single detector, placed 
at a given baseline distance from the reactor core, is used to measure the ab- 
solute flux and energy spectrum of through the inverse beta decay process. 
As detectors have moved further away from the reactors over the years, it has 
become more important to shield them from background sources (primarily de- 
rived from cosmic muons) since the flux of reactor neutrinos falls with the square 
of the baseline distance. Experiments at baseline distances of one kilometer or 
more (Chooz, Palo Verde and KamLAND) have accomplished this by placing 
their detectors at significant depths underground. 

5.1 Detector Target and Buffer 

Any detector to measure reactor antineutrinos takes advantage of the inverse 
beta decay reaction: 

v^p e'^n (14) 

followed by neutron capture and measures a coincidence between the signal from 
the 6+ and the neutron. The prompt positron will exhibit 1-8 MeV of visible 
energy when it annihilates, with a minimum energy from the e"'"e~ masses. The 
inverse beta decay reaction takes place on hydrogen, an element which occurs 
in water and all forms of liquid scintillator. Liquid scintillator has been used in 
CHOOZ, KamLAND and earlier experiments. Since scintillator consists of long 
organic molecules, one issue is to accurately determine the hydrogen fraction. 

A neutron will capture on hydrogen and form deuterium giving gamma rays 
with an energy 2.2 MeV. This is the process measured in the large KamLAND 
experiment. Smaller experiments find it advantageous to load the scintillator 
with about 0.1% Gadolinium (Gd) which has a very large capture cross sec- 
tion for neutrons and also leads to a higher energy gamma shower, 8 MeV. 
In CHOOZ, 86.6 ± 1.0 % of the neutrons were captured in Gadolinium |10| . 
However the Gadolinium lowers the optical attenuation length and, more impor- 
tantly, time dependent effects have been noticed regarding the optical properties. 
Time dependent effects will need to be minimized or thoroughly understood, at 
least in the way that they might manifest themselves differently in multiple de- 
tectors. Gd loaded targets are presently favored for a new reactor experiment, 
but if the time dependent effects cannot be adequately controlled, then a larger 
detector for a future 8000 t GW y experiment might consider not using it. 

There are some issues to be considered in the choice of scintillator. In some 
forms. Gadolinium loaded scintillator is not compatible with an acrylic vessel. 
One example is pure PXE with 0.1-0.15% Gd. The compound Gd-ACAC or 
Gd-Carbolylate is being investigated at Max-Plank-Institut fiir Kernphysik in 
Heidelberg. The aromatic component is CigHis which has a low concentration 
of hydrogen. The alternative is 40% PXE/PC with 60% mineral oil and 0.1% 
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Gd, the cocktail already used by Palo Verde. It is compatible with acrylic, 
though there have been problems of degradation of attenuation length versus 
time seen in certain solutions of as much as 2 cm/day. Aromatics and alkenes 
provide more hydrogen per unit volume than pure PXE. To compare, 10 cubic 
meters of PXE weigh about 10 tons and contain 5.15 x 10^^ atoms of H, com- 
pared with 10 cubic meters in CHOOZ, which weighed 8.54 tons and contained 
6.87 X 10^9, an increase of 33%. 

The glass in the phototubes themselves emits gamma rays from which the 
fiducial volume must be shielded. This will be accomplished with a scintillating 
buffer, a non-scintillating buffer, or both. 

The advantage of a scintillating buffer is that the positron energy is fully 
contained. No neutron will capture on the Gd in the buffer (because there isn't 
any), but a positron which crosses the target/buffer boundary will not lose any 
less visible energy if the buffer is scintillating, so the minimum energy of 1 MeV 
is seen, i.e. Ethreshoid < Eminie'^)- This reduces a systematic error which was 
0.8% in CHOOZ to 0% 

The disadvantage of a scintillating buffer is that it contributes to a greater 
accidental background from '"^K radioactivity present in the PMTs. This would 
argue for the presence of both scintillating and non-scintillating buffers. One 
conceptual design of a detector with two buffers is shown in Figure!^ 

The inner volume of the conceptual 10 ton detector design shown in Figure 

is 1.4 m (radius) Gadolinium loaded scintillator to serve as the neutrino 
target. The size of this volume is determined by the desired target mass, but 
may be constrained by the overall size of the target hall. This volume will 
serve to well define the fiducial volume in which neutron capture takes place. 
The 2nd volume is 0.35 m scintillator without Gadolinium, the gamma catcher. 
This distance is defined by the absorption length of gamma rays and assures 
that all gammas from neutron capture in the fiducial volume are visible. With 
a scintillating buffer, there is no need for a fiducial volume cut. The 3rd volume 
is 1.0 m mineral oil without scintillator, the PMT buffer. This boundary, a few 
photon interaction lengths, buffers the central region against photons coming 
from radioactivity associated with the phototubes. The outside volume from 0.6 
to 1.0 m thick is filled with water has an independent set of phototubes is used 
to veto cosmic ray muons and other entering tracks. Alternatively, there may 
be layers of active and passive material to serve as a cosmic veto. There may 
also be an additional meter of water in some directions to serve as a neutron 
catcher. 



The disadvantage of the second buffer is that it makes a detector larger for a 
fixed fiducial volume. The size of the detector will affect the civil construction 
costs of the detector hall and tunnel, which are expected to be the dominant 
costs for the experiment. This is especially true if movable detectors are con- 
sidered, because the path on which the detector moves will need to be that 
much larger. Therefore the advantage of the buffer will have to be evaluated 
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Figure 23: Conceptual detector design. Inner volume is Gadolinium loaded 
scintillator. The second volume is non-doped scintillator. The third volume is 
mineral oil. Outside is veto region. 

quantitatively together with the construction costs for any site. For example, 
in a design with two buffers, 2.7 m of diameter is needed for the buffers. This 
means that there is a maximum of 15 ton target volume if the detector needs 
to be moved through a tunnel with a diameter of 6 meters. 

5.2 Mechanical Structure 

The structure will be a significant part of the mechanical design of the detector. 
There will need to be an optical structure, to separate the target from any buffer 
regions and another one to separate the buffers from each other. The optical 
structures will perhaps be made from acrylic. KamLAND and Borexino have 
had success using nylon balloons, but a flexible container is probably not a good 
way to maintain fiducial volumes to better than the 1% accuracy required for 
this experiment. There will also need to be an opaque structure to separate 
the buffers from the veto region and to support a frame which holds the PMTs. 
This will perhaps be made from steel. The acrylic vessels will need to have 
access ports through which the volumes can be filled and also in which sources 
can be put in and moved around. All structures will need to be designed so that 
proper buoyancy can be maintained during fill. 
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5.3 Muon veto 



A sufficient volume of water is needed to shield the detector from neutrons 
created by muon cascades outside the veto region. The veto region will need 
to have enough phototubes to efficiently tag all muons which enter the muon 
region. If movable detectors are chosen, it is not foreseen that the veto detector 
would be movable. But it would then need to be designed so that the fiducial 
target could be removed. 

Another design for a muon veto involves layers of passive shielding, such as 
concrete, and a layer or layers of high efficiency planar detectors to measure the 
incoming muon. 

There are reasons to also consider passive shielding outside the veto region, 
depending on the conditions of local radioactivity and the orientation of any 
shafts or tunnels. 

5.4 Liquid Handling System 

There will need to be a safe and careful system for handling and testing the dif- 
ferent liquids. There will need to be very accurate measurements of the volume 
of each liquid during fill, particularly in the target region. Spills of liquid scin- 
tillator will need to be avoided, particularly with regards to any contamination 
of ground water. 

The conceptual design for the detector shown in Figure ESI is accompanied 
with the need to handle four different liquids: Gd loaded scintillator, scintillator, 
mineral oil, and water. All will need to be maintained with well-known optical 
properties. One possible solution for the issue of optical changes in Gd loaded 
scintillator is mixing through recirculation, with or without filtering. A recir- 
culation system will have to be designed to evenly affect each liquid throughout 
its volume. Filtering might be designed to deal with optical degradation due to 
organic compound production, but it may also remove Gd, and a system which 
removes and precisely reloads solutions in the liquid may be difficult to accom- 
plish. Another concern about recirculation is the introduction of microbubbles, 
which could both change the optical properties as a function of time and also 
lead to small time-dependent density changes. The advantages and disadvan- 
tages of recirculating each volume will need to be carefully reviewed, probably 
by extensive testing. 

5.5 Detector Shape 

The ideal detector has a spherical design. This design offers the lowest ratio of 
surface to volume, which implies the least number of photomultipliers per ton, 
the lowest surface radioactivity, the most buffer material to absorb external 
neutrons, the lowest cross section to cosmic rays, and hence the lowest total 
mass per ton of target. With SNO, KamLAND, MiniBooNE and Borexino, 
there is a great deal of experience with the design and construction of spherical 
liquid scintillator detectors. They also have the simplest parameterization of 
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optical effects, wliicli will be even more important due to the low attenuation 
length of Gd loaded scintillator. 

Nevertheless, if movable detectors are chosen, cylindrical detectors offer some 
advantages. The target mass must be moved in an area less than 6 meters 
squared (the maximum size of most economic tunnels). To achieve a 50 ton or 
greater mass, the detector must be cylindrical. As a comparison, in Figure EH 
there is a 50 ton cylindrical target and a 25 ton spherical one. Both detectors 
have a 70 cm pure oil buffer, and 20 cm region to shield the PMTs, together 
with 20% PMT coverage. Assuming that only one near detector is needed in 
each case, then for a 50 ton active volume for the far detector, 2372 PMTs are 
needed for 2 cylindrical detectors compared to 2382 for 3 spherical ones. For 
75 tons, the result is 3156 for 2 cylindrical detectors and 3176 for 4 spherical 
ones. The costs of the vessels cannot be ignored, but in the limit where the 
detectors costs are dominated by PMT and channel counts, this would argue 
toward multiple spherical detectors over larger cylindrical ones. 



5.6 Movable Detectors 

The purpose of designing movable detectors is to be able to move the identical 
far detector to the near detector location for part of the running period and 
have a head-to-head calibration of the relative number of events. Assuming 
that there is no additional error introduced by the act of moving the detectors, 
the uncertainty of this efficiency (e) is 



where N is the number of events measured during the calibration. 

Designing the detectors to be movable will certainly have a great effect on 
the detector design. However, it will have an even larger effect on the tunneling 
requirements. Let's consider the options for an experiment with two detectors 
and a flat overburden. Three scenarios are shown in Figures HOI In 
Scenario One, there are two shafts and a tunnel connecting them in which the 
detector can be moved. The shaft for the near detector may be placed away 
from existing infrastructure such as substations and parking lots, outside the 
security fence. This is probably the easiest solution for detector design, but the 
most expensive for civil construction. The second solution takes advantage of 
sloped access and the fact that the near detector does not need to be as deep as 
the far detector. The angle of approach for the access ramp could be from any 
convenient direction. In scenario three there would be two shafts and a surface 
detector transport procedure. The detector would be raised to the surface by 
a crane or elevator system and transported on surface rails. This may be the 
cheapest option, but it may also be a more complex and risky procedure. 




(15) 
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Figure 24: A comparison of spherical and cylindrical detectors 
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Figure 25: A scenario with two shafts and a connecting tunnel 
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Figure 26: A scenario with sloped access 
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Figure 27: A scenario with two shafts 
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5.7 Multiple Detectors 

Assuming that backgrounds and deadtime issues are under control, the statis- 
tical power of the experiment is driven by the size of the far detector. The near 
detector will have a larger event rate by a factor of so for the same 

statistical power, it might be reasonable to make the near detector smaller. Any 
difference in size between the two detectors, however, compromises the ability to 
cancel systematic errors, and this experiment's goal is to measure the neutrino 
flux with a much better systematic error than has previously been achieved. 
This dilemma has led to consideration of several other detector configurations: 
Multiple Far Detectors It might be desirable, for example, to make a 20 ton 
near detector, and multiple identical 20 ton far detectors, perhaps five. This 
would be a way to get more fiducial mass at the far detector, and also a possible 
way to stage some of the experiment. The multiple far detectors would be a 
smaller fiducial mass than a single larger detector with the same channel count. 
However, they would provide additional checks on some of the systematic errors 
of multiple detectors. 

Multiple Locations Another variation of this idea is shown in FigureESlwhere 
multiple identical detectors are used. In this example, a small 5 kiloton detector 
is placed near the reactor, and another one a few hundred meters away, along 
with a larger 30 ton detector, which is repeated at 1.7 km. Again, there are 
multiple opportunities to study possible systematic errors and staging possibili- 
ties. This scheme provides additional opportunities for cross calibration of each 
detector. 

One possible attraction of the multiple location idea is that the optimization 
of distance for the 400 t GW y experiment could be at a different location than 
for the more sensitive 8000 t GW y experiment. This is because the optimum 
location, discussed in Section |2| is different for the rate test and the energy 
shape test. The rate test is more powerful for the lower statistics and the shape 
test for the higher statistics experiment. Unfortunately, the optima discussed 
in Section 01 does not seem to lend itself to a solution that involves staging 
detectors in a natural way. 
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Figure 28: Possible layout for multiple detectors. 
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6 Calibration 



6.1 Introduction 

Calibration provides the information required to: (a) reconstruct the energies 

and positions of the prompt and delayed sub-cvcnts, (b) time the experiment's 
simulation software to match as closely as possible the data without introduc- 
ing biases and so compute the detection probability as accurately as possible, 
and (c) reliably estimate the magnitude of uncertainties in detector response. 
Although it will not be further discussed in what follows, calibration should 
also enhance the capability to discriminate between signal and background pro- 
cesses, e.g., through accurate characterization of detector response to neutrons 
with energies typical of spallation products and through development of a trigger 
to record muon- produced neutron capture in the experiment. 

What calibration accuracy is required? Based on the record of previous 
experiments similar in design to those being considered for a next generation 
reactor neutrino experiment, the contribution from calibration to the overall 
uncertainty for a single detector can be limited to a few percent. The calibration 
program outlined should aim toward controlling the single detector uncertainty 
at this level. To achieve an overall systematic uncertainty on the relative far- 
to-near detector event rate at the percent level, it is further required that the 
calibration of the near and far detectors-as well as detector design, operation, 
and data analysis-be as identical as possible. 

6.2 Calibration System Design Considerations 

A calibration program capable of achieving the accuracy specified above will 
have the following characteristics: 

• Calibration sources that are precisely deployable throughout the entire 
active target region of each detector. At representative locations in each 
detector, it should be possible for some calibration sources to pass through 
and beyond the outer boundary of the active target. The requirement that 
sources be deployable throughout the entire target region and somewhat 
beyond its boundary at representative locations can be relaxed if events in 
the boundary region are excluded from analysis. Whatever region is used 
for analysis, it should be considered important that calibration sources 
can be placed at any point in this region. The uncertainty in the position 
of the calibration source should be about 5 cm or less, assuming the lineal 
dimensions of the detectors are on the level of several meters; and the 
deployment system and calibration sources should be designed to minimize 
shadowing and absorption. These requirements require novel mechanical 
designs for source insertion and removal that may be quite challenging. 
Indeed, considerable R&D may be required to establish that calibration 
methods successfully demonstrated on the ~12 ton CHOOZ detector can 
be scaled up to 50-|- ton devices being contemplated. 
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• Multiple calibration sources that include the following: 

1. Point gamma sources spanning the energy range from inverse beta de- 
cay threshold to the highest achievable energy for calibrating energy 
estimation, vertex finding, and detection efficiency. Readily available 
radio-isotopes that would be suitable include Cs-137 (0.662 MeV), 
Ge-68 (2 X 0.511 = 1.022 MeV), Zn-65 (1.116 MeV), and Co-60 
(1.333 + 1.173 = 2.506 = MeV). Sources with energies in the range 
of 5-8 MeV would also be desirable, but are typically difficult to fab- 
ricate and deploy due to short half lives and low branching ratios for 
the gammas of interest. A development effort to design and deploy 
gamma sources in this energy range may well be worthwhile. At least 
one of the gamma sources should have its visible activity known to 
~1% for the purpose of measuring detection efficiency, and a range 
of gamma energies is needed in order to help measure the effects of 
quenching and Cerenkov radiation on light output. The Ge-68 source 
is particularly attractive because it probes the detector response at 
inverse beta decay threshold. 

2. j3 source. Using gammas alone, it may prove difficult to determine 
the absolute energy scale for positrons to better than 1%. Some 
studies of future reactor neutrino experiments have assumed that 
the absolute energy scale is known to 0.5%. To achieve this level 
of precision, use of a /3 source, preferably a (3'^ source, should be 
planned. Controlling uncertainties in the [3 energy spectrum from 
capsule shadowing and absorption in the encapsulation materials will 
be a challenge, however. A deployment of a Ge-68 source in solution 
with the liquid scintillator has been suggested |HH1 to solve to this 
problem and it would also provide a truly homogeneous and in situ 
calibration for gamma rays. 

3. Tagged neutron source for measurement of the neutron detection 
efficiency. Am-Be is a good candidate for such a source. The 4.4 MeV 
gamma can be used to tag single neutron emission. Corrections for 
the differences between the Am-Be neutron energy spectrum and the 
inverse beta decay neutron energy spectrum would be made using 
simulation. Important also is that such a source provides additional 
gamma calibration points: (i) 2.225 MeV from capture of neutrons 
on protons and (ii) 4.443 MeV provided that a significant fraction 
of the neutrons can be moderated before they enter the scintillator. 
Moderated Am-Be sources have been used successfully for energy 
calibration in other experiments. The CHOOZ experiment also made 
good use of Cf neutron sources. These produce a continuous neutron 
and gamma spectrum but offer the advantages of higher rate for 
studies of homogeneity. 

4. Variable intensity light flasher. Possible sources include UV lasers 
and UV LEDs. As practically point-like sources of fast pulses of con- 
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troUablc amplitude, such sources are valuable for calibrating vertex 
reconstruction (which is usually based on PMT timing with charge- 
dependent corrections), measuring PMT gain and interpolating the 
detector response between gamma calibration points. 

• Integration of the calibration system as an integral part of the data ac- 
quisition, event reconstruction, and detector simulation software: A large 
fraction of events read out over the course of the experiment will be vari- 
ous calibration triggers. The acquisition of this data must not adversely 
affect dead time for the experiment. Careful design of triggering and 
run control is essential, and processing of calibration data must take place 
quickly to help identify shifts in detector performance. 

• Simulations. Because it is not possible to deploy a source that mimics in- 
verse beta decay in the reactor neutrino energy range, detector simulation 
ultimately plays an important role in estimating the detector response to 
this process. As such it requires an accurate and detailed description of 
the detector (regarding materials, geometry, and optics) and of the cal- 
ibration sources and the devices used to deploy them. Generation and 
transport of photons to the PMT photocathodes must be modeled ac- 
curately, as should also be the digitization of PMT pulses and effects of 
discriminator thresholds. Simulated data must be reconstructed with the 
same software as real data. The event generators for inverse beta decay 
and the calibration processes must be accurate and detailed. 

• Construction of identical calibration programs for the near and far detec- 
tors in both hardware and software. So far as is possible, the calibration 
sources should be the same, or at least fabricated by identical means; and 
the deployment systems should likewise be identical in design. If the near 
and far detectors have different sizes, designing the deployment systems 
to be the same to within overall scale should be a consideration. The 
scope and magnitude of the calibration programs should be planned so 
that any calibration carried out for one detector can also be carried out 
for the other at nearly the same time. 

6.3 Concluding Remcirks 

What is ultimately required for adequate calibration will not be fully known 
until experience has been gained on the detectors running under normal condi- 
tions. However, given past experience and the experimental design and goals, 
it is clear that the baseline calibration program should have the capability of 
precisely deploying well-understood calibration sources at any desired position 
throughout-and preferably beyond, in the case of gamma and neutron sources- 
the active volume of the detector. Control of systematics will be aided by the 
ability to deploy the same sources in both the near and far detectors. The cali- 
bration program effort should be matched by an effort to develop a full detector 
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simulation tuned and checked against the calibration data; this simulation pro- 
gram can then be used to extend understanding of the detector into areas for 
which calibration data is not available. 
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7 Detector Overburden and Backgrounds 



Backgrounds will be an important consideration in site selection and detector 
design. They can be classified into two distinct sources: (1) internal back- 
grounds associated with radioactivity in the materials of the detector and lab 
site, and (2) external activity associated with cosmic rays and their interaction 
daughters. The former is a well-understood technology that has been solved 
many times over several decades by a combination of materials selection, self- 
shielding, and signal enhancement via addition of neutron absorbers that boost 
the energy of the neutron capture reaction. It is not expected that this will be 
a major problem. 

Cosmic ray activity is typically more serious and in most cases requires an 
underground location of at least at least a few tens of meters of water equivalent 
(m.w.e.) for even rather small detectors of less than 1 ton. The limiting factor 
is typically the flux of neutrons associated with muon nuclear interactions in 
the detector and in the surrounding rock. Table [S] is a brief summary of some 
previous reactor experiments done over the last 25 years. As can be clearly seen, 
longer distances have necessitated larger detectors, which have in turn required 
underground sites to reduce the cosmic ray muon flux and hence the flux of 
external neutrons. 



Table 5: Summary of some previous reactor experiments 



Name 


powert/i 
(GW) 


mass 
(tons) 


distance 
(m) 


depth 
(mwe) 


comments 




ILL IHSI 


0.057 


0.32 


8.76 


7 


■^He + scint. 


. PSD 


Gosgen [HDl 


2.8 


0.32 


38/46/65 


9 


same as ILL 




Rovno [SH 


1.4 


0.43 


18 




Gd scint. 




Krasnoyarsk 


1.6 


0.46 


57/231 


600 


'^He only 




Bugey 3 HSHm 


2.8 


1.67 


15/40/95 


23/15/23 


^Li + scint. 


PSD 


Savannah R. 


2.2 


0.25 


18.2/23.8 


- 10 


Gd -1- scint. 


PSD 


CHOOZ 10] 


4.4 


5 


1000 


300 


Gd -I- scint. 




Palo Verde jTS] 


11.6 


11.3 


800 


32 


Gd + scint. 




KamLAND 13 


- 80 


408 


- 180,000 


2,700 


scint. 





The final design of this experiment will be a trade-off between several con- 
flicting goals: (1) the desire to have the near and far detectors be as identical 
as possible to allow cancellation of systematic errors due to detector geometry, 
(2) the desire to have the near detector be small to avoid muon-related back- 
grounds and expensive overburden requirements, and (3) the desire for a far 
detector of size large enough to give an acceptable event rate. Backgrounds will 
be a major consideration in the final design and so we have developed some 
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general constraints based on the measurements of previous experiments. 



7.1 Types of Backgrounds 

The double coincidence afforded by + p ^ n + e+ followed hy n + p-^d + j 
requires that two events appear within a given time and space window, typi- 
cally 25-200 ^s and < 1 m. The space requirement is usually enforced either by 
vertex fitting or detector segmentation. Events that are inherently "singles" are 
therefore strongly suppressed as they must appear by accidental coincidence in 
the correct window. Examples are internal radioactive decay, stopping cosmic 
ray muons, and external gamma events. 

Correlated backgrounds are more serious. Such events can mimic the reactor 
anti-neutrino signal and there can be significant uncertainties associated with 
the rate. The majority of correlated backgrounds come from events associated 
with cosmic ray activity. Examples are fast neutrons from muon nuclear inter- 
actions, muon capture reactions from stopped muons, and excited nuclei from 
muon spallation reactions that may de-excite by neutron emission. 

In many cases experiments have been placed deep enough underground to 
sufficiently suppress the cosmic ray muon flux (e.g. KamLAND). In other cases 
a segmented detector was used to allow more precise cuts on vertex location 
and partial detector vetoing (e.g. Palo Verde and Bugey). In cases where the 
detector is small and close to the core, only minimal shielding was required (e.g. 
Savannah River). Table shows the effective detection rate for some previous 
experiments. This is the daily rate of neutrinos detected (corrected to a distance 
of 100 m) while the detector is operating per ton of fiducial target per GW of 
reactor power ^ . Variations of a factor of two are seen, mostly due to the overall 
efRciency constraints imposed by the necessity to reject backgrounds. Detectors 
which are segmented and near the surface (ILL, Gosgen, Bugey, Palo Verde) 
tend to have less efficiency than monolithic ones (Savannah River, CHOOZ, 
Rovno, Krasnoyarsk, KamLAND) with sufficient shielding relative to the detec- 
tor size. 

In the discussions below we use the experience of previous reactor experi- 
ments and estimate depth and shielding requirements for both a "near" (100 m) 
and a "far" (1.5 km) detector. It is assumed that the basic detector module is 
roughly 10 tons, with the far site requiring several such modules and the near 
site one. 



^There are differences expected on the order of 5% due to fuel composition 
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Name 



effective rate at 100 m 



ILL 



15.7 
14.4/14.5/13.6 



Gosgen 
Rovno 
Krasnoyarsk 
Bugey 3 
Savannah River 
CHOOZ 
Palo Verde 
KamLAND* 



22/19/19 



46 
35/41 



50/54 
51 
12 
61 



Table 6: Effective detection rate. (*) The KamLAND rate is based on the 
expected rate of 86.8 events in 141 days. 



7.1.1 Accidental Coincidences 

Single, uncorrelated events may result from radioactive decay from internal 
materials such as steel and glass, external gammas from the shielding or sur- 
rounding laboratory, or from cosmic ray muons if the depth is shallow. The 
coincidence rate R for two uncorrelated events with rates Ri and i?2 to occur 
in a time window w is given by ^: 



This rate must be compared with expected neutrino interaction rates for 
a typical 10-ton scintillator detector. For the near (far) module this would 
typically be about 1000 (4.5) events/day/GW with 100% efficiency. A typical 
reactor is assumed to be 3 GW thermal, giving a signal rate of 3000 (13.5) 
events/day. Using eauationll6l and setting a criterion that the accidental rate 
should be less than 1% of the actual event rate, restrictions on various sources 
of "singles" can be estimated. 

Cosmic Ray Muons: Cosmic ray muons can produce accidental coincidences 
in several ways: (1) straight-through muons mimicking positrons or neutron 
capture (unlikely), (2) decay of stopped muons giving an followed by an ac- 
cidental coincidence of a muon capture that gives a neutron, and (3) a stopped 
muon decaying to an followed by the accidental coincidence of a neutron 
from spallation. In (2) and (3), the rate must be less than the rate of acci- 
dentals from (1) since they are initiated by somewhat rare reactions from (1). 
Correlated background from these types of events is treated in a later section. 



This assumes that Rw <C 1 



R = R1R2W 



(16) 
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Figure 29: The rate of muon accidental coincidences as a fraction of the expected 
signal as a function of depth. Left (right) is for the near (far) detector. 



Taking Ri — R2 = i?^, where is the muon rate through the detector, and 
a typical neutron capture time of 25/l(S which would imply a window of roughly 
100 lis for 98% efficiency, then allowed muon rate through the detector is 1.9 
Hz (0.12 Hz) for the near (far) detector. This is to be adjusted for muons that 
can be tagged. Most all previous experiments had muon vetoes (see Table CJ 
that were better than 95% efficient in tagging muons. Therefore the acceptable 
muon rates can be conservatively adjusted up by at least a factor of 20. This 
would correspond to a depth of about 85 m.w.e. (350 m.w.e. for far) for a 
10 ton module. Figure EHl shows the percentage of the signal this background 
would be as a function of depth. 

This estimate is conservative in that it does not take into account the low 
probability that a muon is misidentified as a positron and/or capture gamma. In 
addition, vertex cuts may further reduce the possibility of a chance coincidence. 
What can be said is that at these depths single muon accidentals will not be 
a problem as long as the neutron capture time is kept short by the addition of 
gadolinium or other absorber. 

Internal Radioactivity: Backgrounds from internal radioactivity come from many 
sources: (1) U/Th in detector materials (such as scintillator), (2) U/Th and 
'^^K activity from PMT glass, and (3) U/Th in stainless steel, to mention a few. 
Radon incursions may also take place from penetration through piping flanges, 
calibration device insertion, or any place where detector scintillator contacts the 
radon-laden air common in underground sites. 
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Table 7: cosmic ray rejection 



Name muon rate veto efficiency dead time 



^l!^ TO TO 

ILL 250 99.8 8 

Gosgen 260/300/340 99.8 8 

Rovno 350 99 7 

Bugey 3 99.5 2 

CHOOZ 1 2 

Palo Verde 2000 96 2 

KamLAND 0.3 92 11 



The singles rate from these sorts of backgrounds is a strong function of the 
detector threshold. This is because most long-lived beta activity is below 2 
MeV and energetic alphas are strongly quenched, typically with factors of 0.1 
or less such that they have visible energies of less than 1 MeV. The effect of 
these events on the coincidence rate is therefore a strong function of threshold. 
Table IHl shows the threshold for recent reactor experiments. 1 MeV is clearly 
achievable if sufficient care is taken in materials selection, radon exclusion, and 
the allowance of sufficient shielding from PMTs via a non-scintillating buffer 
region and self-shielded fiducial volume. As an example, assuming a 238[/ igvel 
of 10"^** g/g in the scintillator (about 100 times less stringent than the lev- 
els required by Borexino) and only taking into account the 3 beta decays with 
Q > \ MeV results in a singles rate above 1 MeV of about 0.004 Hz in a 10 ton 
detector. This is conservative in that it assumes all the betas are above thresh- 
old. This level is consistent with a more detailed study done for the KamLAND 
detector in the U.S. proposal. In this study similar rates are obtained for ^^^Th 
and ^^K at levels that can be achieved by reasonably careful control of detector 
materials and possibly scintillator washing and/or filtration. Care is also taken 
to move the PMTs well back from the fiducial volume. 

Radon- laden air in enclosed spaces can sometimes reach as high as 1200 
Bq/m^ |96,. Since the ^'^'^Rn chain includes ^'^^Bi{Q=S.272 MeV) allowing 
such levels inside the scintillator would result in unacceptable background lev- 
els. In addition, ambient radon levels are often seasonal and/or dependent on 
local weather conditions. It is therefore necessary to control radon levels in the 
laboratory areas and to take precautions to ensure radon cannot diffuse into the 
scintillator via use of radon- free nitrogen blankets, radon impermeable o-rings 
and gaskets in all piping, and the use of detector materials low in radon emis- 
sivity. For example, the radon level in the Super-Kamiokande water is reduced 
to around 3 x 10~^ Bq/rrt' by taking such precautions. Similar low levels in a 
10-ton detector would result in a coincidence rate of around 0.003 c?~^, which 
is completely negligible. 
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Table 8: e"*" energy measurement performance of past reactor experiments. 
Name threshold precision resolution variation 



(MeV) (%) (%FWHM) (%) 

ILL 0.9 2 18 5 

Gosgen 0.7 2 18 5 

Rovno 0.7 1.5 30 1 

BugeyS 118 1 
Savannah River 2 

CHOOZ 1 15 
Palo Verde 1 

KamLAND 0.9 2 8 1 



External Gamm,a Activity: Gammas from U/Th decay chains in the rock sur- 
rounding the detector can be a source of background. Typical concentrations in 
rock can be as high as 10~^g/g. This necessitates placing a buffer region around 
the detector which efficiently shields these gammas. In addition, steel shielding 
may be necessary if the gamma rates are very high. In the KamLAND detector 
the calculated external gamma rate is about 0.2 Hz in a fiducial volume roughly 
5 gamma attenuation lengths from the rock wall. Assuming similar shielding for 
a 10-ton detector and scaling by the surface area results in an external gamma 
rate of about 0.05 Hz, or a coincidence rate of about 0.005 d~^, which is neg- 
ligible. Obviously, most of these events occur on the edge of the detector and 
the rate changes dramatically with the amount of shielding. Three attenuation 
lengths instead of five would increase the coincidence rate to about 14 d~^. Care 
must be taken to measure the rate of external gammas at a potential site and 
design the shielding accordingly. 



7.1.2 Correlated Backgrounds 

Correlated events that could be mistaken for a reactor neutrino event can come 
from many different sources: (1) direct cosmic ray hadrons, (2) muon decay from 
stopped muons, (3) gammas from muon capture on carbon followed by neutron 
emission, (4) spallation products from muons, most notably those which may 
de-excite via neutron emission, and (5) fast neutrons which elastically scat- 
ter off protons in the scintillator and are subsequently moderated and captured. 
All these are dependent upon the cosmic ray rate and energy spectrum at depth. 

Direct Nucleons: Cosmic ray protons and neutrons at sea level are not from 
the primary flux but are all essentially from secondary interactions in the at- 
mosphere. The flux is highly peaked in the vertical direction (~ cos™ 9, where 
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m ~ 8) and can be estimated by the expression : 

IniE) = (3 X 10'5)(^/ioO)-^ 

where I is the differential vertical neutron intensity in cm"^ sr~^ M eV~^ , 
E is the kinetic energy in MeV, and 7 is a spectral index which is a weak func- 
tion of energy, about 2.4 in the range around 1 GeV and softening to 2.9 at 
higher energies. 

The effective neutron attenuation length (A) has been determined for many 
types of rock in association with geological dating studies that rely on produc- 
tion of rare elements via neutron interactions 98 . Values range from about 120 
- 170 g cm~^ . Conservatively assuming all the flux to be vertical and taking the 
longest attenuation length one can then estimate the flux of these neutrons un- 
derground above 10 MeV to be less than 10~^m~^c?~^ at 30 mwe. Even taking 
a factor of ten for the uncertainty due to secondary production, these nucleons 
are negligible for a 10 ton detector below this depth. 

At depths less than 30 mwe it will be necessary to perform a more detailed 
calculation, taking into account secondary production, attenuation length of the 
actual shielding material, and the flux expected at the particular geomagnetic 
cutoff. 

Stopped Muon Decay and Capture: Stopped muons can contribute to correlated 
backgrounds in two ways: (1) prompt muon ionization signal followed by muon 
decay, and (2) muon capture on ^^C to produce ^^B. 

The stopping muon rate can be estimated as a function of depth by looking 
at the change in the vertical muon rate. Figure ISUl shows the estimated muon 
stopping rate in a 10-ton module as a function of depth. It is assumed that the 
stopping power difference between shielding and oil is a factor of two. 

The rate of untagged muons that decay in the detector is given by: 

R={1- e)RJdkfEe''-^^ 

where is the muon rate, e is the efficiency of the muon veto, fdk is the fraction 
of those muons that stop and decay, /b is the fraction of the decay electrons in 
the range of reactor neutrinos (~ 10%), is the veto time after any event, and 
T is the effective lifetime of muons in oil (~ 2.126 /is). For ty = 5 fj,s and a 95% 
veto efficiency there would be about 330 untagged stopped muons/day. Clearly 
one must improve the muon recognition and also accept signal loss due to dead 
time. 30 mwe was roughly the depth of the 11-ton Palo Verde detector, which as 
TableElshows did have a poor efficiency compared to other experiments. Going 
to 90 mwe would reduce this to an acceptable 30 events/day, or 1% of signal 
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Figure 30: The rate of stopping muons in a 10-ton detector as a function of 
depth. Left (right) is relevant for the near (far) detector. 

rate for a near detector. To achieve a 1% level for the far detector would require 
a depth of roughly 600 mwe. 

In oil about 7.9% of the /i~ will capture on ^'^C. Of these, only about 20% 
will go to the ^^B ground state, the rest going to particle unstable states that 
typically result in neutron emission. Such neutron-producing reactions are a 
very dangerous background since they are correlated with the initial muon and 
can look like a reactor neutrino event. We correct the rate formula for untagged 
stopped muons (ignoring Je) for the fraction of /i~ in cosmic rays (= 0.44), 
and the capture fraction to particle-unstable states, (0.44) (0.079) (0.8) = 0.028. 
Thus at 30 mwe there will be roughly 190 untagged captures/day that result 
in a correlated neutron. The rate at 50 mwe would be roughly 33/day, which 
is acceptable for the near detector. For the far detector a depth of 400 mwe is 
needed to reach the 1% level. 

Muon-Induced Spallation: Most spallation products that can be made by muons 
passing through scintillator are short-lived and/or do not make double-coincidence 
events that can be misconstrued as reactor neutrino events. In that case the 
coincidence rate must be less than the coincidence rate of single, untagged 
muons, or they simply appear as a low level increase in the overall singles 
rate. In cither case they do not present a problem. There are two exceptions: 
(1) ^He{ti/2 = 119ms) ^ n + Li with a branching ratio conservatively 

estimated at 16%; and (2) ^Li{ti/2 = 178ms) e~ +^ Be + n followed by 
^Be 2a. In these cases the electron or two quenched alpha's constitute the 
prompt signal and in both cases a ~ 1 MeV neutron is produced to make a 
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delayed signal (and also add to the prompt signal via n-p elastic scattering). 

The production rate of such events is very difficult to determine, as muon- 
induced spallation is not well-understood theoretically, though there are ex- 
perimental measurements from a muon beam at CERN j99j . However, most 
previous reactor experiments '' have measured their correlated backgrounds via 
reactor on/off data and such events would be lumped together with the back- 
ground from fast neutrons, described below. 

Fast Neutrons: Fast neutrons come from muon-induced nuclear reactions on 
either on the scintillator itself or on the rock. These neutrons then elastically 
scatter off protons in the scintillator and produce an energetic proton "prompt" 
signal. They subsequently capture in the usual way, mimicking a reactor neu- 
trino event. The spectrum of proton-recoil prompt signals is fairly flat out to at 
least 50 MeV and has been measured by the CHOOZ experiment at 300 mwe 
to be about 0.031 d~^ton~^MeV~^. Extrapolation of this rate in the reactor 
prompt signal region from 1-8 MeV would give a rate of 2.2 events/day for a 10 
ton detector at 300 mwe of similar geometry. 

Comparison of different experimental results on neutron production per 
muon per g ■ cm^^ of target (i?) at depths (x) from 20-5000 mwe can be roughly 
fitted by the empirical function TITO : R = lO.O.x""!'^ x 10"^ n-g ^cm?/^. Thus 
one may extrapolate the CHOOZ result to different depths. For the near detec- 
tor a depth of 50 mwe is sufhcient to reduce the rate to roughly 30 events/day 
(assuming the same resolution and efficiency of CHOOZ). For the far detector 
a depth of 1100 mwe would be required to reach 0.125 events/day. This is a 
very stringent requirement which may be difhcult to achieve in practice. Thus 
the background due to fast neutrons will have to be well-understood at depths 
shallower than this. 



7.2 Summary of Overburden Requirements 

The final result of the experiment will depend on the measurement of the reac- 
tor flux to high precision. To achieve a similar relative systematic error in the 
flux measurements of the near and far detector they should have roughly the 
same background as a percentage of the expected signal. Thus the requirement 
of similar S/N for the near and far detector results in different overburden re- 
quirements. To keep the background from any single source less than 1% such 
that it is on the order of the required precision in the flux measurement the near 
detector should be at a depth of 90 mwe or greater. 

For the far detector, a depth of 1100 mwe would be desirable, but may not be 
available in practice. At shallower depths the understanding of the background 

^The exception is the KamLAND detector, which estimated roughly 2 untagged events/year 
in 400 tons at a depth of 2700 mwe. 
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from fast neutrons will likely be a limiting factor. At depths less than 600 mwe 
accidental coincidences from untagged stopped muons may become important. 
More restrictive cuts than those described here or a more efficient veto will likely 
have to be devised. At 400 mwc the correlated background due to /i~ capture 
will begin to become important. This is potentially a serious background which, 
as in the case of stopped muons, can be reduced by better muon tagging. 
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8 Systematics 



8.1 Overview 

A general discussion of systematic errors in a two-detector experiment has been 
presented in Sections 13.41 and l3. 51 For the sake of definiteness, we consider here 
a restricted class of experiments, as follows: 

• There are at least two identical detectors (one near, and one or more far) . 

• At least one detector might be movable between near and far locations. 

• There may be more than one reactor. 

• The total nmubcr of antineutrino events in the far detector (s), although 
larger than in previous experiments, is not large enough that spectral- 
shape information materially improves the result from total rate informa- 
tion. 

This is essentially a counting experiment, and we will ignore errors in spectral 
shape that do not contribute to overall normalization. 

As discussed in earlier sections, the great advantage of a near-far detector 
configuration is that systematic errors in knowledge of antineutrino flux and 
antineutrino interaction cross-section are absent or greatly reduced, and that 
detector-related systematics are reduced to those due to differences in detectors. 
It may be possible to further reduce detector-related systematics by moving one 
far detector to the near position so that rates can be compared directly, as 
described in Section 15.61 However, the act of moving may in itself be a source 
of systematic error due to changes in a detector during the move. These effects 
are best addressed after the sources of error in a non-moving detector are fully 
understood. 

In this section, we will look at a particular detector configuration and at- 
tempt to list and estimate the sources of systematic error. We consider the 
generic detector geometry shown in Figure 1231 The detector fiducial volume 
is monolithic (like CHOOZ and KamLAND) rather than segmented (like Palo 
Verde) , and liquid scintillator is contained in hard- walled spherical or cylindrical 
acrylic vessels rather than in soft- walled balloons. There is a muon veto system 
surrounding each detector, which would remain in place if detectors were moved, 
implying more than one veto setup for any movable detector. The volume viewed 
by the inward-looking photomultipliers (PMTs) that detect the positron from 
antineutrino interactions and the subsequent neutron capture has three concen- 
tric (or coaxial) regions: an innermost volume (1) containing Gadolinium-doped 
liquid scintillator, whose light transmission will probably get worse with time; 
a "gamma-catcher" region (II) containing undoped liquid scintillator; and a 
buffer region (111) containing non-scintillating oil. Region 111 both isolates the 
scintillating regions from radioactivity in the PMTs and also does not permit 
scintillation light to be generated very close to any PMT. 
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selection 


e{/o} 


rel. error (%) 


positron energy 


97.8 


0.8 


positron-geode distance 


99.9 


0.1 


neutron capture 


84.6 


1.0 


capture energy containment 


94.6 


0.4 


neutron-geode distance 


99.5 


0.1 


neutron delay 


93.7 


0.4 


positron-neutron distance 


98.4 


0.3 


neutron multiplicity 


97.4 


0.5 


combined 


69.8 


1.5 



Table 9: Summary (from Reference ^01) of the detection efficiencies in the 
CHOOZ experiment and their systematic uncertainties. 

For detector-related errors, we will be guided by the one-detector experience 
of the CHOOZ experiment 10 to identify the sources of systematic error that 
are potentially most serious, and will briefly discuss how much they might be 
reduced in a near- far experiment with identical detectors. We also postulate 
an event reconstruction method similar to that used by CHOOZ [TOI: both the 
position and the energy of the prompt (positron) and delayed (neutron-capture) 
subevents of each event are obtained from a fit to the distribution of charge 
signals from the large number of PMTs in the detector. This method relies on 
good knowledge of the calibration of each PMT and of the time-dependent light 
attenuation in each detector region. 

The detection efficiencies and their systematic uncertainties for the CHOOZ 
experiment are listed in Table |5| the combined systematic error is 1.5%, unac- 
ceptably large for future near-far experiments. 

8.2 Potential sources of systematic error 

For convenience, we list potential sources of systematic error in four general 
categories: "physical" errors, errors from triggering and data recording, errors 
from event-selection cuts, and errors from background subtraction. 

8.2.1 Potential sources of "physical" errors. 

This category includes errors in source and target parameters and in losses of 
primary reaction products (e^, n). 

Potential errors related to reactor sources (s). 

1. Distance from each reactor core to each detector. CHOOZ |69j quoted 
a precision of ±10 cm. If the near detector in a near-far experiment is 
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at 250 m, a similar precision would introduce a negligible uncertainty of 
0.08% in l/r2. 

2. Relative flux from each reactor core. This error cancels for the case of two 
cores and symmetric placement of detectors. The general case of multiple 
reactors and multiple detectors is discussed in the following subsection. 
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Potential errors in number of free protons. 

1. Mass of fiducial region. It should be possible to weigh each kind of liquid 
added to each detector to ~ 0.1%. Each acrylic- walled region of each 
detector would then have a thin standpipe to monitor the excess volume 
and its fluctuations due to temperature change. 

2. Fraction by weight of free protons. CHOOZ found that determining this 
quantity from combustion measurements was not trivial, and quotes an 
uncertainty of ±0.8%. With multiple detectors the issue is whether the 
free proton content is the same in all detectors. One could, for example, fill 
all detectors equally from each batch of mixed scintillator. This relative 
uncertainty should be reducible to the ~ 0.1% level. 

3. Effective fiducial volume. The question here is whether the boundary of 
the innermost (Gd-doped) detector region (I) actually delimits the fiducial 
volume. This would happen, for example, if all detected e"*" subevents from 
region I were accompanied by detected n subevents, and no e"*" subevents 
from region II were accompanied by detected n subevents. This will cer- 
tainly not be true. However, to first order, "spill-in equals spill-out": the 
loss of e+n events near the outside edge of region I is compensated by 
a gain of e+n events near the inside edge of region II. The difference in 
non-compensation among detectors should be small. 

An alternative is to make an explicit cut on the fitted positions of the e+ 
and n subevents, as was done for a single detector by CHOOZ (see Sec- 
tion below). To bound this error, we will will use their uncertainties 
resulting from this method. Whether using no position cuts at all is a 
superior procedure is clearly a subject for detailed simulation studies. 

Loss or absorption of reaction products. This category is coupled to 
the event-selection cuts because of the multiple gammas present in an event: 
two annihilation gammas from the positron, and typically three from neutron 
capture on Gd. These deposit visible energy by successive Compton collisions. 
The most likely symptom of gamma loss is thus a reduced positron or neutron 
signal rather than none at all; this generates event inefficiency via an energy 
cut. We consider here the rarer situations in which positrons or neutrons are 
lost entirely. 

1. Loss of positrons or neutrons by absorption in inert material (acrylic ves- 
sels) or through openings such as chimneys or source pipes. One can rely 
on simulations to indicate that this effect is small, and if so, that it essen- 
tially cancels when detectors are compared. 

2. Loss of neutrons escaping beyond the active region before they are captured. 
This effect should not be sensitive to small differences in detectors, and 
should nearly cancel when detectors are compared. 
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3. Capture of neutrons on the wrong material. Neutrons can be captured 
on hydrogen instead of on Gadolinium. In this case, the capture gammas 
have a total energy of only 2.2 MeV instead of « 8 MeV, and would be 
rejected by typical energy cuts on the delayed event (6 to 12 MeV). 

This is a serious source of systematic error in a single-detector experiment. 
CHOOZ used a tagged ^^^Cf neutron source in several positions in 
their detector to calibrate this effect. They quote an efficiency for n- 
capture on Gd of (84.6 ± 0.85)%, a 1% systematic. 

In a multiple-detector experiment, in addition to measurements with a 
tagged neutron source, one can compare neutron counting rates when 
the same untagged neutron source is placed in similar positions in each 
detector. The relative error among detectors should be reducible to a few 
tenths of a percent. 

8.2.2 Potential errors from triggering and data recording 

Potential errors in the time domain. Wc assume that all detectors will 
be interlocked so that down-time due to malfunctions such as high- voltage trips 
will be the same for all. 

1. Deadtime losses from the muon veto and from event recording. These 
losses should be small and easily measurable. Note, however, that if the 
muon-shielding overburdens of the detectors are very different, the dead- 
time losses could be substantially different. 

2. Trigger time window for the delayed event. Care should be taken that the 
hardware cut is looser than any anticipated software cut. 

Potential errors in the pulse-height domain. CHOOZ's lowest-level trig- 
ger required a minimum analog sum of pulse heights and a minimum number 
of PMTs firing at the « 0.5 photoelectron level. Care must be taken that such 
requirements do not compromise clean cuts on fitted energy, especially since 
periodic changes in trigger thresholds must be made due to increasing light at- 
tenuation in the Gd-doped scintillator. Most of CHOOZ's quoted systematic 
error due to positron energy (0.8%) came from this effect. 

Potential errors from differences in detector electronics. Just as mul- 
tiple detectors should be physically identical, so must their electronics be identi- 
cal in performance. Identical pulser-driven "pseudo-events" should be delivered 
frequently to all detectors and analyzed as data. 

8.2.3 Potential errors from event-selection cuts. 

The CHOOZ analysis ^01 imposed six criteria for selecting events. Four of these 
depend on the results of the position-energy fit to each event described above. It 
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is clear that in a near- far experiment, it is essential to ensure identical behavior 
of this fit for all detectors. 



1 . Fitted positron energy Ee < S Me V. The presence of the region II gamma- 
catcher guarantees that even when the positron has zero kinetic energy, 
the 1.1 MeV of energy from its annihilation gammas will be seen. No 
lower cut on positron energy is therefore needed, removing a potentially 
scrioiis source of systematic error. The high-side energy cut Ee < S MeV 
eliminates a negligible number of real events, and generates an uncertainty 
of only 0.05% even in a single-detector experiment. 

2. Fitted neutron energy 6 < En < 12 MeV. This cut is made to reduce 
background. Partial escape of capture gammas gave it an efficiency of 
(96.4 ± 0.4)%, as measured from the energy yield of neutrons from ra- 
dioactive sources. The relative error among detectors due to the loss of 
the small low-energy tail should be much smaller than this. 

3. Fitted positions: minimum distance from photomultipliers of subevents. 
(a) Fitted position of positron (prompt subevent) at least 30 cm from the 
"geode" of PMT faces {de+ > 30 cm; (b) fitted position of neutron (de- 
layed subevent) also at least 30 cm from the geode (d„ > 30 cm). These 
requirements were motivated for CHOOZ in large part by the ungraceful 
performance of the fitting procedure when light was deposited very near 
a PMT. It should not be necessary at all in a detector geometry with a 
non-scintillating oil buffer layer between the scintillating regions and the 
PMTs. CHOOZ assigned systematic uncertainties of 0.1% each to these 
cuts. One might guess that the the uncertainty in the difference in perfor- 
mance between detectors should be at least a factor of two smaller than 
this. However, the detector design should preserve a position resolution 
comparable to that of CHOOZ (ct w 4 cm at 8 MeV). 

4. Fitted positions: maximum distance between subevents. Difference in fitted 
positions of positron and neutron less than 100 cm (dg+n < 100 cm). This 
cut was made to suppress background, and it may well be desirable in 
a future experiment. CHOOZ used Monte Carlo simulation to establish 
its efficiency at (93.7 ± 0.4)%. The relative error among detectors in this 
quantity should be several times smaller. 

5. Time window for neutron subevent. Time delay of neutron subevent rela- 
tive to prompt subevent 2 < t^+n < 100 /is. CHOOZ imposed the cut on 
the low end of this window because of overshoot problems in AC-coupled 
photomultipliers (which can be avoided in future experiments); it caused a 
loss of (1.6 ±0.2)%. The high-end cut covers about three capture lifetimes 
in Gd-doped scintillator, and had an efficiency of (95.3 ± 0.3)%. Since 
there is no reason for this cut to depend strongly on details of detector 
performance, the relative error among detectors should be considerably 
smaller. 
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6. Neutron multiplicity: Nn = 1. CHOOZ imposed this cut to reduce the 
correlated background from muon spallation, which typically generates 
more than one neutron. It rejects some real events when the accidental 
presence of a gamma ray from natural background fakes the second neu- 
tron. The efficiency of this cut was (97.4 ± 0.5)%. Although the relative 
effects among identical detectors should be reduced, care must be taken if 
the background environments of the detectors are substantially different. 

8.2.4 Potential errors from background subtraction. 

The fraction of background in the near and far detectors will be different be- 
cause of the much higher real event rate in the near detector. Subtraction of 
background in the far detector(s) must therefore be well understood. 

A general discussion of backgrounds is given in Section |7| Briefly, there are 
two types: accidental backgrounds, in which the prompt and capture subevents 
are randomly associated; and correlated backgrounds, in which the subevents 
are not random (but may not be in the proper order). At 300 mwe, CHOOZ 
observed an accidental background rate of (0.42 ± 0.05) d^^, and a correlated 
background rate of (1.04±0.11) d~^. Their real event rate at full reactor power 
was 24.7 d~^, so the background fraction was 6% and contributed a systematic 
uncertainty of 0.5%. 

It should be noted that due to reactor commissioning, CHOOZ had almost 
as much reactor-off time to study backgrounds as they had reactor-on time, a 
circumstance unlikely for future experiments. It is also important to note that 
their background rates changed significantly whenever thresholds were changed 
to accommodate increased light attenuation in Gd-doped scintillator. 

Let us assume for a moment that in a future experiment each far detector 
will have about 5 times the fiducial tonnage of CHOOZ, but that they are likely 
to be located about 1.7 times farther away from a source of comparable power, 
and that background rate scales with detector mass. Under these assumptions, 
the background-to-signal ratio will be 1.7^ = 3 times worse than for CHOOZ, 
~ 18%. To keep the systematic error from background subtraction below 0.5% 
under these assumptions, it will be necessary to know the background to better 
than 3% of itself, more than three times as well as CHOOZ did with ample 
reactor down-time. Serious thought must be given to reducing background and 
to inventing ways to measure its spectrum and magnitude. 

8.3 Summary - potential systematic errors 

For a near-far experiment, the systematic errors in rate comparison between 
detectors associated with reactor sources and number of target protons can be 
kept below a few tenths of a percent with reasonable attention. Extrapolating 
from CHOOZ experience, the errors arising from differences in detectors in per- 
formance and for event-selection criteria can probably be reduced about a factor 
of three from the CHOOZ figure of 1.5%, provided that sufficient care is taken 
in experimental design and execution. The presence of gamma-catcher and inert 
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buffer regions appear to be essential elements of the design, and understanding 
the behavior of the position-energy fit to each event under conditions of decreas- 
ing scintillator transparency is crucial. Finally, the systematic error associated 
with background subtraction could well be the hardest to keep under control 
because of the limited beam-off time available to measure it. 
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8.4 Multiple-reactor scenarios 

In the following parts of this section, the effect of multiple reactors and multiple 
detectors is considered. The Kashiwazaki site, for example, has seven reactors 
in two groups. It is shown that the uncertainty due to the nature of multiple 
reactors can be made small, if the near detector is placed in such a way that 
the contribution from each reactor of the cluster is approximately equal. 

8.4.1 Reactor induced systematics - multiple reactors, one detector 

To develop notation, we consider first the simplest case involving one reactor 
and one detector. The effective systematic error is given by 

'^eff = <^u+ (^c + [<^c ) + ) 

where (T„ is the uncorrelated error of the detector, Cc is the correlated error of 

(r) (r) 

the detector, cjc is the correlated error of the flux from the reactor and ct„ is 
the uncorrelated error of the flux from the reactor. 

8.4.2 The case with one detector + multiple reactors 

When the number of reactors is larger than one, (^^ff becomes 

a=l ^ ^ 

where Tq (a = 1, • • • , N) stands for the yield from the a-th reactor and T = 
Sa=i is the total yield. If the yield from each reactor is equal, i.e., = 
T/N, then the contribution of the uncorrelated systematic error (cri'^-')^ from 

the reactor in Equation ^ becomes (cri^^)^ Y.a=i (Ta/Tf = {1/N){ai''^)^ , so 
that we have 

a,^,,=a,^+a,^ + (aM)2 + l(aMf. (18) 
Therefore, we see that the more reactors there are, the smaller the contribution 

(r) 

of the uncorrelated systematic error tri becomes. This is because the average 
of independent N fluctuations is smaller than a single fluctuation. 

8.4.3 The case with multiple reactors and detectors | 101 | 

It has been known that the correlated error is canceled in the case of a single 
reactor experiment with near and far detectors. Now a question arises: what 
happens to this cancellation in the case of an experiment with multiple reactors 
and detectors? To answer this question, let us consider the ideal case with N 
reactors and {N + 1) detectors, where each reactor has a near detector in its 
neighborhood and it produces the same yield at a far detector (See Figure I^T)) . 
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Figure 31: The ideal case with N reactors (R) and (A^ + 1) detectors (D). 
In this case the effective systematic error becomes 



where the conditions = /N and T^'' = S^^T'^'' have been assumed; 
and T^'^ are the yield at the far and 6-th near detector which is close to the 6-th 
reactor (6 = 1, - ' ' i ^) from the a-th reactor, and and are the total 
yield at the far and 6-th near detector. As in the case with one reactor, the 
dominant contribution to (Jeff comes from the uncorrelated error CTu, and the 
contribution of the uncorrelated error au'' is reduced in H19|) by a factor of N 
due to the averaging over the independent N fluctuations, but this reduction is 
irrelevant because such an effect comes in the correlated error which is almost 
canceled in the multi detector system. To conclude, the answer to the question 
at the beginning of this section is that the cancellation of the correlated error 
occurs also in the ideal case with N reactors and {N + 1) detectors. It should 
be noted that the number N of the near detectors in this case is sufficient but 
not necessary to guarantee this cancellation, as we will see below in the case of 
the Kashiwazaki plan. 

8.4.4 Kashiwazaki case study 

The next question which arises is: what happens to the cancellation of the 
correlated error in the Kashiwazaki plan. The Kashiwazaki-Kariwa nuclear 
plant consists of two clusters of reactors, and one cluster consists of four reactors 
while the other consists of three (See Figure|23). Before we discuss the effective 
systematic error for the actual Kashiwazaki plan, let us consider the ideal limit. 




il + l/N)al 



(19) 



1 + (1/iV) {l + (iV + 1) [a? + (ar^)2 + (ai''))Viv] /al} 



-1 ' 
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Figure 32: The Kashiwazaki plan with 7 reactors (R) and 3 detectors (D). 



Q® 



Figure 33: The ideal limit of the Kashiwazaki plan, where one reactor has the 
power four times as much as one of the Kashiwazaki-Kariwa reactors, while the 
other has the power three times as much as that of the KK reactors. 
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where the cluster of the four reactors shrinks to one reactor and the other of 
three reactors shrinks to another single reactor (See Figure lHS]) . It can be shown 
analytically that the effective systematic error for this ideal limit is given by 



^eff 



7 



where the correlated error has been canceled due to the near-far detector com- 
plex, and the reason that we have the factor a/ 74/49 instead of(l-|-l/A^)^/^| n=2 = 
•y/3/2 is because the yield from the first reactor to that from the second is 4:3 
instead of 1:1. Thus we see that the cancellation of the correlated errors occurs 
also in the ideal limit of the Kashiwazaki plan. Furthermore, it can be shown 
from the numerical calculations that the difference between a^f f for the actual 
Kashiwazaki plan and that for the ideal limit is very small: 



'^eff\ actual KK 



= 1.04, 



'^'^ffUdeal KK 

where cr„ = 0.6%, CTc — 1.6%, — 2.3%, ai^^ — 2.5% have been used as 
reference values for the systematic errors. This shows that the cancellation of 
the correlated errors occurs also in the actual Kashiwazaki plan. 

The systematic limit on sin^ 26*13 at 90%CL can be obtained from the effec- 
tive systematic error by 



(sin^ 29) 
where 
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stands for the averaged expectation value of the factor sin^ (^Am^L/AE), L-^ 
(X=near-1, near-2, far) is the distance between the detector X and the a-th 
reactor (a = 1, • • • ,7), and cu^ (X=near-1, near-2, far) is the fraction of the 
yield from the a-th reactor at the detector X (=(power of a-th reactor)/ (i^)^). 
The expectation (• • • ) is defined by 
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where e{E), f{E), a{E) stand for the detection efficiency, the neutrino flux, and 
the cross section, respectively. If |Am|]^| = 2.5 x 10~'^eV^, then we have 

sin^ ("^^^M-^O.TO 



AE 

W (^^^^^Jl^^^' -0.17 

in the setup of the Kashiwazaki plan. So, if o'„=0.6%, then Equation l|2()|l 
becomes 

U\n^2e)TT'' = ^—0.006 = 0.022. 

Furthermore, if we include the statistical error Ostat = l/\/60, 000 in the case 
of 20 t-yr, then we get 



(sin^ 2^?);/^;"'^ = :^V(74/49)0.0062 + 1/60,000 = 0.025. 



In Figure r^ a) comparison is given between the actual Kashiwazaki plan and 
a hypothetical experiment, which is depicted in Figure 03 with a single reactor 
and two detectors, where the same uncorrelated systematic error (T„ and the 
same data size (=20t-yr) are assumed (Only in this figure the statistical errors 
are taken into account). We observe that there is little difference between the 
sensitivities of the actual Kashiwazaki plan and the hypothetical experiment 
with a single reactor at |Am§]^| = 2.5 x 10~^eV^. Also it is remarkable that the 
sensitivity of the actual Kashiwazaki plan for higher value of |Am|]^| is better 
than the single reactor experiment. This is exactly because of the reduction of 
the uncorrelated error due to the nature of multiple reactors (cf. Equation (|18|l ). 
where the near detectors play a role as far detectors. On the other hand, to see 
how effectively the correlated errors are canceled in the actual Kashiwazaki plan, 
comparison is given in Figure IMT b') between the locations of the near detectors 
in the actual Kashiwazaki plan and a hypothetical plan which is depicted in 
Figure with one near detector very close to reactor #1 in the first cluster 
while the other detector is very close to reactor #5 in the second. In the 
hypothetical case, the sensitivity is deteriorated because the correlated errors 
do not cancel completely. From these two figures we see that the setup of 
the Kashiwazaki plan is almost optimized and it does not suffer from an extra 
uncertainty because there are more than one reactors. 

8.4.5 Energy spectrum analysis jl02j 

One can also derive a semi-analytic formula for the effective systematic error in 
the analysis of the energy spectrum, assuming that the uncorrelated bin-to-bin 
error ct^ is independent of bin. The result is 
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Figure 34: (a) Comparison of the actual Kashiwazaki plan (the solid line) and a 
hypothetical experiment (the dashed line), which is depicted in Figure with 
a single 24:.3GWth reactor and two detectors assuming the same uncorrelated 
systematic error cr„. (b) Comparison of the actual Kashiwazaki plan (the solid 
line) and a hypothetical experiment (the dashed line), which is depicted in 
FigureOni where the near detectors are placed in wrong location. The statistical 



errors as well as all the systematic errors except (t„ 
2.3%, are ignored for simphcity. 



(r) 



which is assumed to be 



Figure 35: The case where a single reactor has the power seven times as much 
as one of the Kashiwazaki-Kariwa nuclear plant. 
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Figure 36: A hypothetical experiment with the same location of the reactors 
as the Kashiwazaki plan but with stupid location of the near detectors. 



where n stands for the number of bins, CTc is the total correlated error, (Jshape 
is an uncertainty of the theoretical prediction for each energy bin which is un- 
correlated between different energy bins, cr^'') is the uncertainty of the flux, and 
the error of the energy calibration has been ignored because it turns out to be 
small. a„ and 6„ are given by 




XI - {i + i)cj 



(i = 1, - • • ,"-), 

and their numerical values are given in Table IIUI The integration over the 
energy in the expectation value to define Cj is from Ej to -Ej+i. The systematic 
limit on sin^ 2^13 at 90%CL is obtained from cTg/ / by 
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n 


an 




2 


151 


44 


4 


76 


20 


6 


60 


15 


6 


51 


12 


16 


41 


10 


> 30 


- 40 


~ 9 



Table 10: The coefficients a„ and fe„ which were calculated numerically. 



n 


(a„/2n)i/2 


6„/(2n + 2) 


2 


6.2 


7.3 


4 


3.1 


2.0 


6 


2.2 


1.1 


6 


1.8 


0.7 


16 


1.1 


0.3 


62 


0.6 


0.1 



Table 11: The numerical values of the coefficients (a„/2n)^/^ and 6„/(2n + 2) 
which appear in cTe/ / in the case of ctc = cr„ . 



For simplicity let us assume aghape = and y^ct*^''' '3> \/nac, <Ju- Then we have 

\ / J, 2 \ -1/2 

V 2n V 2 (T^ + ncr^y 
If we further assume (t„ = then (|21|1 is reduced to 

and the numerical values of hn/{2n + 2) are given in Table [TTl which shows that 
the contribution of the correlated error is not negligible for lower value of n. On 
the other hand, if we assume ^/nac ^ ctu, then we have 




and it indicates that the uncorrelated error gives the dominant contribution to 
the effective systematic error. The same conclusion is reached if ^/ncJc <^ cr^ is 
assumed. Therefore the realistic value of the uncorrelated bin-to-bin error fT„ 
has to be estimated carefully in any case. 
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9 Possible Sites 



In this section we examine the power performance of commercial reactor sites 
all over the world and summarize the assets of sites that have been identified as 
candidates to host a sin^ 2^13 experiment. 

Site selection will involve a comprehensive look at many criteria such as 
potential for earth shielding, locations for near and far detectors, construction 
costs and reactor power. Ultimately, choosing a site will require a a comparison 
of cost and sensitivity for site specific proposals. 

9.1 Top Performing Reactors Worldwide 

One characteristic on which various reactor sites can be directly compared is 
power. Tables IT^ 1131 and 1 1 41 list maximum and average power output for single, 
double and multi-reactor sites respectively. The reactor sites are sorted by their 
average power output as measured over the 7 year period from 1996 to 2002. 
The maximum rated power of each site is also listed. Sites with maximum power 
ratings of less then 3 GW per reactor core (or 9 GW for multi-reactor sites) are 
considered only if they have been identified as potential host sites. 



Reactor Site 


Country 


Avg MWth 


Max MWt 


Brokdorf 


Germany 


3900 


4214 


Emsland 


Germany 


3892 


4097 


Grohnde 


Germany 


3858 


4184 


Grand Gulf 


US 


3505 


3833 


Grafenrhcinfeld 


Germany 


3357 


3936 


Wolf Creek 


US 


3211 


3565 


Perry 


US 


3199 


3758 


Callaway 


US 


3176 


3565 


Leibstadt 


Switzerland 


3130 


3511 


Waterford 


US 


3152 


3390 


Watts Bar 


US 


3049 


3411 


Unterweser 


Germany 


3117 


4126 


Seabrook 


US 


2924 


3411 


Vandellos 


Spain 


2882 


3181 


Kruemmel 


Germany 


2868 


3851 


Confrentes 


Spain 


2858 


3160 


Hope Creek 


US 


2794 


3339 


Fermi 


US 


2750 


3430 


River Bend 


US 


2676 


3039 


TriUo 


Spain 


2672 


3119 


Columbia 


US 


2567 


3486 


Tokai 


Japan 


2086 


3219 


Krasnoyarsk 


Russia 


1600(?) 


2000(?) 



Table 12: Power performance for single reactor sites around the world jl()3lll()4] . 
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9.2 Reactors Sites Under Consideration 



Many possible host reactor sites have been identified by groups around the world. 
In this section we list some of the relevant features from each site mentioned for 
consideration. 

9.2.1 Angra, Brazil 

The Angra reactors are located on the coast of Brazil a few hours south of Rio 
de Janeiro. The site is ringed by sharply rising hill on three sides, with summits 
of over 200 meters within 1.5 km. 

The president of the reactor company is a former particle physicist and 
groups at several Brazilian institutions have shown an interest. 

9.2.2 Braidwood, Byron and La Salle, Illinois 

These reactor sites are all located within 100 km of both Fermilab and Argonne 
National Laboratory, making them easily accessible to the particle physics com- 
munity. The topology of northern Illinois is generally flat, so obtaining cosmic 
shielding at these sites will require extensive excavation below the surface level. 
Fortunately, the subterranean geology of this region is well understood jl()5| . 
At the Braidwood and La Salle sites there is a layer of dolomite, which is well 
suited for excavation, at a depth of 110 meters (approximately 300 mwe). At 
the Byron site, the western most of the three, the dolomite layer is much closer 
to surface. Below the dolomite, starting at a depth of 50 to 100 meters, is a layer 
of sandstone, which may be poorly suited for supporting excavated cavities. 

The Braidwood, Byron and La Salle reactors are owned and operated by the 
Exelon Corporation. Currently, negotiations are ongoing with Exclon for the 
use of these sites. 

9.2.3 Chooz, Prance 

This is the location of the previous 10 ton CHOOZ experiment, and a site for 
a new near detector has been identified near the reactors. The CHOOZ plant 
has a total 8.4 thermal GW and the CHOOZ-far lab is located 1.05 km from 
the two cores. 

9.2.4 Cruas, Prance 

The Cruas site is located in south eastern France. The most likely near detector 
baseline is about 1 km with a shielding potential of greater than 200 mwe. A 
far detector located at 1.8 km could get shielding greater than 400 mwe. 

9.2.5 Daya Bay, China 

There are two twin reactor cores located in Guangdong province near Hong 
Kong. The total reactor power is 11.6 GW. A third twin core is planned to 
be online in 2010. A near hall with an overburden of 400 mwe is potentially 
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Reactor Site 


Country 


Ave MWth 


Max MWth 




US 


6864 


7600 


1 1 1 V 


i 1 ciiUjC 


6799 


9135 




H yo n PP 


679'! 


8879 




v_.J V . 1 illdll y 


6734 


7865 






6491 


71 79 


V tie 


TIS 


64^6 


71 30 

I -LOU 


Byron 


us 


6442 


7172 


JJixJWliO -L l^ii Y 


TIS 


6377 


fiQl fi 


T liiTiPrir'K 


US 


6365 


6916 


Iscir 


( -rPi'in fi nv 

V. 1 V I iiicllly 


6313 


6985 




TIS 








u o 








T^rfi npp 


6197 


8088 


P ni liTnT^ niiTO' 

J. 11111 JJJJD U 111 


vjci Illdll y 


6187 


6976 


^iiQmiPnJiTi'nfi 

kJ LlD^J Ll'^llojlilld 




61 61 


6978 




H yo n PP 
± 1 dllljC 


61 36 

U -LOU 


7977 


1 P ^" CI Tir Pt CI 


ITS 


(SI 1 


(SS99 




JT 1 dlUjC 


(Sill 


7Q77 




TIS 


\J\J\J ± 


UO I \J 


x^ldUiU vycJjilvuil 


US 


6043 


6749 


(;r»Tn 51 TIP n A PppV 


TIS 


5986 


691 6 

Ui/ -LU 




jj 1 diiijt: 


^Ql 


OUOZi 


Neckar 


Germany 


5881 


6452 


McGuire 


US 


5880 


6822 


Flamanville 


France 


5879 


8088 


Biblis 


Germany 


5528 


7388 


Asco 


Spain 


5496 


6013 


Belleville 


France 


5377 


7977 


Kuo-Sheng 


Taiwan 


4749 


5764 


Angra 


Brazil 


4547 


5873 


Indian Point 


US 


4467 


6096 


La Salle 


US 


4323 


6978 


Salem 


US 


4281 


6918 


Ignalina 


Lithuania 


3985 


8778 


D.C. Cook 


US 


3281 


6661 


Millstone 


US 


3271 


6111 



Table 13: Power performance for double reactor sites around the world |1()3[ 
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u 








lia Tl a n a 

V_> CtilCtU-Clj 
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Chiiion 


France 


4 


8653 


11166 


Blayais 


France 


4 


8644 


11131 


Cruas 


France 


4 


8586 


11190 


Takahama 


Japan 


4 


8439 


9925 


Genkai 


Japan 


4 


8330 


10177 


Kori 


S. Korea 


4 


8314 


9203 


Ringhals 


Sweden 


4 


8307 


10841 


Tricastin 


France 


4 


8284 


11178 


Bruce 


Canada 


4 


8080 


10710 


Tihange 


Belgium 


3 


8075 


9127 


Hamaoka 


Japan 


4 


8031 


10584 


Forsmark 


Sweden 


3 


7773 


9408 


Dampierre 


France 


4 


7753 


10967 


Bugey 


France 


4 


7728 


10897 


Leningrad 


Russia 


4 


7642 


11705 


Balakovo 


Russia 


4 


7520 


11705 


Kozloduy 


Bulgaria 


6 


6618 


11002 


Kursk 


Russia 


4 


6577 


11705 



Table 14: Power performance for multi-reactor sites around the world |103[I104] . 
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located about 300 m from the core, and a far hall with 1200 mwe at a distance 
of 1500-2000 m to the core. 

9.2.6 Diablo Canyon, California 

The Diablo Canyon site is located on the coast of California. Within 1.5 km of 
the reactors are hills of over 300 meters. All the surrounding land is owned by 
the utility, PG&E, with which negotiations for the use of the site are ongoing. 

9.2.7 Kashiwazaki-Kariwa, Japan 

The Kashiwazaki-Kariwa reactor complex is the most powerful in the world. The 
site consists of 7 active reactors spanning approximately 700 meters. The site 
is relatively flat, so detectors would be located at the bottom of vertical shafts. 
Permission for the use of the site has been granted by the power company. A 
full description of a proposal at this site is given in Appendix A. 

9.2.8 Krasnoyarsk, Russia 

The Krasnoyarsk reactor is part of a fully underground facility at a depth of 600 
mwe. Halls suitable for detectors exist at baselines of 115 meters and 1000 me- 
ters. Like the Wolf Creek site this single reactor site will allow for a sensitive 
determination of the background rate. 

Russian physicists have a long-standing relationship with the facility and 
have mounted several neutrino experiments in the past |1()6| . 

9.2.9 Kuo-Sheng, Taiwan 

The Kuo-Sheng reactor complex is located in a hilly region of coastal Taiwan. 
The existence of a road tunnel within 2 km of the reactor indicates that tunneling 
is viable in the hills surrounding the site. Within 500 meters of the reactor there 
are hills 50 meters or higher suitable for shielding a near detector. Far detector 
shielding of 150 to 200 meters of rock is possible. 

Local physicists have an existing relationship with the reactor operators and 
are currently operating a neutrino magnetic moment experiment at a distance 
of 28 meters from one of the two reactor cores jl07) . 

9.2.10 Limerick and Peach Bottom, Pennsylvania 

These reactor sites are also owned and operated by Exelon. The topology of 
Pennsylvania is somewhat more hilly than in Illinois. Therefore, the potential 
for significant earth shielding with horizontal access may exist at these sites. 
In particular. Peach Bottom has elevation variations of greater that 60 meters 
within 1.5 km of the reactors. 

9.2.11 Penly, Prance 

The Penly reactor complex is embedded in 120 meter high chalk cliffs. 
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9.2.12 Wolf Creek, Kansas 

The topology of this single reactor site is similar to the flat sites in Illinois. At 
a depth of 110 meters there is relatively wide layer of shale. While the shale is 
not ideal for tunneling it is likely sufficient to support the structures required 
for this project. 

The main attraction of the Wolf Creek site is that as a single reactor facility 
it will get full reactor off running which can be used to measure the background 
rate. 
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10 Other physics 

The main goal of the experiment described in this white paper is to search for or 
further hmit the value of 6i^. It is reasonable to explore whether there is other 
physics that can be done with this experiment, or with modest enhancements. 
In this section we briefly consider three other physics topics: sterile neutrinos, 
9i2 and reactor physics. 

10.1 Sterile Neutrinos 

The discovery of sterile neutrinos would have a revolutionary impact on neutrino 
and particle physics. If there is non-negligible mixing of the electron neutrino 
with sterile neutrinos, they would contribute to the reactor antineutrino disap- 
pearance pattern. An idea how to look for sterile neutrinos at reactors along 
with sin6'i3 was proposed at the Kurchatov Institute in 1998 |l()8j . 

The notion of sterile neutrinos was originally introduced by B. Pontecorvo 
in 1967 ["109 and has been considered later by many authors, e.g., D. Caldwell 
and R. Mohapatra J.T. Peltoniemi, D. Tommasini and J.W.F. Valle [TTT] . 

S. Bilenky, C. Giunti and W. Grimus [TTl, K. Benakh and A. Smirnov [TT^ . 
B. Kayser |114) . Information on the theory of sterile (and mirror) neutrinos 
and references can be found in the recent paper by V. Berezinsky, M. Narayan, 
F. Vissani |TT5| . 

While solar, atmospheric, and laboratory (Super Kamiokande, SNO, Kam- 
LAND) studies are understood in the framework of only 3-active neutrino mix- 
ing (see, however, de Holanda and Smirnov |116| ) they do not exclude some 
admixture of sterile neutrinos [Q. An experimental hint in favor of sterile neu- 
trinos comes from the unconfirmed observation of the LSND collaboration |118| 
on z/^ — > Vf. transitions. In particular the so-called (3-1-1) and (3-1-2) neutrino 
schemes, which have been considered to explain the LSND signal, predict reac- 
tor neutrino disappearance with a Am^ ~ eV^ very close to the current upper 
bound from the Bugey experiment (see References |119| and |12(J) '). 

10.1.1 The effect of sterile neutrinos in (?i3 reactor experiments 

In the standard 3-active neutrino framework, antineutrino disappearance at 
distances L = 1000 — 2000 m from a reactor source is governed by Am^^j^ ~ 
2 X 10"'^ eV^ and by the mixing parameter sin^ 26*13 : 



Let us denote the three standard neutrino mass states by {i — 1,2,3), such 
that the atmospheric and solar mass-squared differences are given by Am^^^ = 
m| — m\ and Amg^^j = rn^ — m\ ^ 7 x 10~^ eV^, the solar mixing angle is 
determined by tan^ 6',oi = |C/e2|V|f/eiP - 0.4, and sin^ 26*13 = 4|C/e3p(l - 
|t/e3p)- If we now assume that there are additional neutrino states {i = 




(22) 
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4, 5, . . .) with masses such that^ |Amfj^| ^ 10"'^ eV^ we know from the Bugey, 
Palo Verde and CHOOZ experiments that the mixing of the electron neutrino 
with these new mass states has to be small (see shaded region in Figure 
i.e., 

\Ue^\^<S:l for i>3. (23) 

To first order in these small quantities one obtains for the survival probability 
at nuclear reactors 

P..^..-l-45:|C/..Psin^(^) . (24) 

i>3 ^ ^ 

In the case of only one additional neutrino one can write Equation (|24|l as 

P,,^,, « 1 - sin^ 2013 sin^ (^%^) - sin^ 20, sin^ (^%^) ' (25) 

where sin^ 29 , « 4|f7e4p and Am^^^ = — are the mixing parameters of 
the sterile neutrino, and sin^ 26113 ~ 4|C/e3p- 

From Equation l|25|) it is obvious that such a sterile neutrino would have 
some impact on an experiment of the Kr2Det type |12H or its modifications 
(as discussed in 2002-2003 at the meetings in Paris, Alabama and Munich), 
provided the associated mixing parameter sin^ 26s is not too small and the 
mass-splitting is in the relevant range. 

• If Am^p^ ~ Am^^jj^ it will be rather difficult to disentangle Pe — > I'fj./r 
oscillations with sin^ 29i3 from Vg Vs oscillations with sin^ 20s at a 
reactor experiment, especially if the main information comes from the 
total rate measurement. Only if enough spectral information is available 
and/or Ato^^^j and Am^^^ differ sufficiently one might be able to separate 
the two channels. If no effect is observed both, sin^ 20i3 and sin^ 26*5 can 
be constrained. 

• If Ato^(,^ Am^tm and the oscillations with Am^^^ are already averaged 
out at the near detector position at L ^ few x 100 m, no information about 
sterile neutrino mixing can be obtained from the comparison of the far and 
near detectors, and the transitions to the sterile state will not affect the 
013 measurement. In that case information on sterile neutrino mixing can 
be obtained from the near detector if relatively precise information on the 
initial reactor neutrino flux is available, or if a "very-near" detector at 
L 10 m could be installed (see below). 



*Note that any new mass-squared difference with Am^ <^ 10 eV^ will have no effect in 
reactor neutrino experiments with L ~ 2 km. 
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Figure 37: Example of a layout. Detector positions, target scintillator masses 
and overburden (ni.w.e.) are shown. 



Distance [m] 


Target mass [t] 


i^e rate/300day 


50 


5 


1 100 000 


300 


5 


30 000 


300 


30 


190 000 


1700 


4 X 30 


24 000 



Table 15: Detector positions, scintillator target masses and i/g detection rates 
per 300 days. 



10.1.2 Example of layout 

Imagine that a tunnel is built near one 3.2 GW thermal power reactor. We 
consider five identically designed 30 ton target scintillator (movable) detectors, 
four of them stationed in the far position at a distance of 1700 m from the 
reactor, and one detector stationed in the near position at 300 m from the 
reactor. To expand the explored mass parameter region toward larger values 
two small detectors are considered at 300 m and 50 m from the reactor. A 
possible layout of such an experiment is illustrated in Figure |23 The expected 
neutrino detection rates per 300 days are shown in Table IT^ 

We consider two types of data analysis: shape and rate. With only one 
reactor as source the shape analysis (as we already know) is independent 
of the exact knowledge of reactor power, energy spectrum of the Ue, flux and 
spectral time variations, target volumes and proton concentrations, and absolute 
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efficiencies of Ve detection. Backgrounds can be measured periodically during 
reactor-off periods. On the other hand, the analysis based on the comparison 
of the far/near Df. total rates requires good knowledge of the ratios of target 
volumes and detection efficiencies. In both cases no exact information from the 
reactor services on reactor power and fissile fuel composition is needed for the 
data analysis. 

10.1.3 Sensitivity 

In Figure 051 we show the sensitivity for sin^ 2^13/5 for the experimental config- 
uration described in the previous section. Within 3 years of data taking (300 
days/year) in a large part of the Am^ range from 0.001 to 0.5 eV^ a sensitivity 
of sin^ 26*13/^ ~ 0.01 — 0.02 can be obtained |122| . which is in general agreement 
with the analysis performed by P. Huber, M. Lindner, T. Schwetz and W. Win- 
ter |30| . The limits shown in Figure 1381 were obtained by assuming an energy 
resolution ue = and the systematics o-shape = 0.5%, tJiatc = 1%- As can 

be seen in the figure the CHOOZ hmit on sin'^ 2^13 at Amatm — 2 x 10^"^ eV^ 
can be improved by a factor 10. 

We conclude that a search for sterile neutrinos at reactors does not require 
much additional effort and can be done along with 613 . With one reactor and a 
number of detectors a high sensitivity to ^13 and sterile neutrino mixing can be 
reached. To increase the sensitivity in the range Am^ ^ Arn^^^ [in particular to 
reach the O(eV^) region indicated by the LSND evidence] a very- near detector 
at i ~ 10 m is necessary. 

10.2 Solar Neutrinos and 612 

10.2.1 The present and near-future for 9i2 

The first data from the KamLAND experiment 7 showed that the reactor i>e 
undergo oscillations on a distance scale of ^ 160 km. Under the plausible 
assumption of CPT-invariance, this result confirmed the Large Mixing Angle 
(LMA) solution of the solar neutrino {vq) problem. When combined with the 
solar neutrino data, the KamLAND data split the allowed LMA solution region 
of I'Q oscillation parameters Am^i and sin^ 612 into two sub-regions - low- 
LMA and high-LMA. The best fit values of Am^i and sin^ 9i2 in the two sub- 
regions read: Am|i = 7.2 x 10"^ eV^, sin^ 9i2 = 0.3 (low-LMA), and Ato^^ = 
1.5 X lO"'' eV^ sin^^ia = 0.3 (high-LMA) ^2M- Adding the salt phase data 
from the SNO experiment |6j the global analysis allows the high-LMA solution 
only at the 99.13% CL (nS]. Thus, the high-LMA solution is disfavored by the 
current data compared to the low-LMA one, although it is still not ruled out 
comprehensively. 

The S'e survival probability in the reactor experiments of interest, Pgg, de- 
pends on Am2i, sin^0i2, Arn^^ (driving the atmospheric neutrino oscillations), 
the angle 9i3 limited by the CHOOZ and Palo Verde experiments, and on the 
type of neutrino mass hierarchy |126[ I127| . The potential sensitivity of a reac- 
tor experiment to each of these parameters depends crucially on the baseline 
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Figure 38: Expected 90% CL limits from the comparison of detectors at dis- 
tances of 300/1700 m and 50/300 m. The shaded region is excluded by the 
CHOOZ, Palo Verde and Bugey experiments. 

of the experiment. Experiments with a baseline L ^ (1 — 2) km can be used 
to get information on sin^^is: on the indicated distances only oscillations in- 
duced by Am^j are operative and sin^ 2^13 determines their amplitude. For 
baselines L <; 50 km, the oscillations due to Am|j^ average out and we have: 
Pee ~ [1 — sin^ 2^12 sin^(Arn^L/4iJ)] cos^ ^13. The oscillations generated by 
Am2i have been seemingly detected in the KamLAND experiment. 

In the absence of oscillations, the maximal contribution to the signal in an 
experiment with reactor Ue detected via the inverse /?— decay reaction comes 
from Ue with energy E ~ 3.6 MeV. For a fixed Am^i, maximal sensitivity to 
sin^ 6*12 can be achieved if for E ~ 3.6 MeV, L is "tuned" to a z?e survival proba- 
bility minimum (SPMIN), i.e., if sin^(ATO^L/4iJ) « 1. The latter is reflected in 
the 6+— spectrum, measured in the experiment. The corresponding minimum in 
the spectrum is determined by P^e « 1 — sin^ 20i2, and thus is very sensitive to 
the value of sin'^ 26*12. If, in contrast, L is such that sin^ {Am^ L / 4E) = e <C 1, 
Pee would have a maximum (SPMAX): Pee ~ 1 — esin^ 2^12 w 1. In this case 
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the sensitivity to sin^ 20i2 is worse than in the preceding one. The positions of 
the extrema in both cases are highly sensitive to the value of Am2i. 

For values of Aml^ from the low-LMA region, Am^^ ^ 7.2 x 10~^ eV^, 
and E ~ 3.6 MeV, the SPMIN and SPMAX take place at L ~ 70 km and 
L ^ 160 km. For the KamLAND experiment, the most powerful reactors (the 
Kashiwazaki complex) are located at a distance of ^ 160 km, unfortunately 
close to the SPMAX for low-LMA. The values of Am^i and sin^ 612, allowed at 
99% CL by existing and prospective KamLAND data have been determined in 
References |128l 1129] . The uncertainty in Arn^i, determined using only the i/0 
data, reduces from 68% to 30% after the inclusion of the first KamLAND data 
in the analysis, while that in sin^ 9i2 does not change, remaining rather large 
(29%). The uncertainty in Amfi would further diminish to 9% (6%) after 1 kTy 
(3 kTy) data from KamLAND. However, there will be little improvement in the 
precision on the value of sin'^ 9i2 with the increase of KamLAND statistics. 

The L best suited for measuring 612 if Amfi lies in the low-LMA region is ~ 
70 km |128| . For a reactor complex with a power of 24.6 GW (e.g., Kashiwazaki) 
and data of 3 kTy from a KamLAND- like detector at L ~ 70 km, sin^ 612 can 
be determined at 99% CL with a ^^10% uncertainty p28' . 

A new reactor power plant, Shika-2, is expected to start operations in Japan 
in March 2006. It will be located at i ~ 88 km from KamLAND and will have 
a power of ~4 GW. This baseline is close to the "ideal" one of L 70 km. The 
implications of the new source of Pe on the KamLAND sensitivity to Arri^i and 
s\v? 9i2 were studied in |130| . It was concluded that due to averaging effects 
of the Ve fluxes from the Kashiwazaki and Shika-2 reactors, the sensitivity of 
KamLAND to sin^ 9x2 would not improve, while its sensitivity to Avn? would 
diminish. 

If contrary to the trend emerging in the solar neutrino experiments the next 
KamLAND results conform to a point in the high-LMA region, one would need 
an intermediate baseline reactor experiment with L ~ (20 — 30) km to get a 
SPMIN in the resultant e+ -spectrum [T^ [T77| . It was shown in [1771 that 
with an experimental set-up at intermediate L ~ (20 — 30) km from a reactor 
with power of 5 GW and 3 kTy of statistics, one could measure both Am^i 
and sin^ 9i2 with a ~ 3 — 7% error at the 99% CL. If in addition the detector 
has a sufficiently high energy resolution and sin^ 6'i3 ^ 0.03, one could observe 
the Am|]^— driven subdominant oscillations. This could be used to measure also 
Am|]^ with a high precision, and even to get information on the neutrino mass 
hierarchy [HHlinTj . 

10.2.2 Role of a new Experiment for 9\2 

The goals of a reactor experiment at 70 km and 2 km from a reactor are substan- 
tially different. However, it is worthwhile to consider the location of reactors 
when one is considering the site of an experiment. It might be useful to consider 
whether a particular reactor has reasonable sites both 2 km and 70 km from the 
core or multiple cores. It is also reasonable to consider whether multiple and 
movable neutrino detectors could have any role in a future 9i2 experiment. 
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10.3 Reactor Physics 



Any near detector for an experiment to measure will measure the flux and 
energy distribution of the reactor neutrinos with a greater accuracy than has 
been done before.^ This will allow comparison with both thermal power and 
reactor fuel loading measurements and calculations. We are not currently aware 
that there are any important checks that can be made on reactor design. How- 
ever, we will continue to work on this as a possibility and a possible benefit to 
the reactor companies. Another application could be the direct check of nuclear 
non-proliferation treaties. 

10.4 Supernovae Neutrinos 

A large vat of scintillator will be sensitive to antineutrinos from a galactic su- 
pernova in the 10 to 50 MeV range. This is higher energy than most neutrinos 
from a reactor. A 50 ton detector is not large, but could be sensitive to a portion 
of our galaxy. Sixteen events would be expected from a supernova at 10 kpc. A 
search for these events will probably require an accurate clock and a separate 
trigger. If simulations show that a supernova can be uniquely identified online 
with no background, a reactor detector experiment could join the SuperNova 
Early Warning System, of potentially invaluable use to astronomers. ^° 



^This might be also of some relevance for other reactor neutrino experiments without near 
detectors like KamLAND. A precise determination of the reactor neutrino spectrum could 
help to reduce the error implied by the flux uncertainty for such experiments. 
^"See http://cyclo.mit.edu/snnet/. 
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11 Tunneling 



11.1 Introduction 

This section provides guidance on the underground siting of the ^13 facihties. 
Requirements, design and construction issues are discussed and the basic phases 
in the planning and construction of a tunnel are outlined. The near-term need 
for site-specific geotechnical data to support project development is emphasized. 
Site-specific geologic and geotechnical data is needed to support the identifi- 
cation and study of candidate alignment(s) at the different sites. This same 
data will also serve to support the early selection of construction methods and 
means and the development-associated work schedule(s) and cost estimate(s). 
At present only a limited amount of geotechnical information has been gathered 
on individual sites so comments are necessarily general in nature. This section 
is limited to a discussion of the more general aspects of tunnel construction in 
rock. Shaft excavation is discussed in Section 111.51 and life safety issues are 
addressed in Section [12.21 of this White Paper. 

11.2 Factors Impacting Rock Tunnel Behavior 

A basic understanding of the factors that influence tunnel behavior is a prereq- 
uisite to the achievement of a practical and cost-effective underground design. A 
tunnel design developed with due regard to the constraints of the construction 
process will ensure realistic requirements setting, result in a cost effective design, 
and ultimately help provide for a more affordable and lower risk construction 
product. Conversely, if the design requirements do not pay due attention to the 
particular constraints of the underground construction process, they arc likely 
to place unreasonable demands on the contractor, mining methods and means 
and/or the rock mass. These unreasonable demands will translate to an increase 
in costs and risk that the project and/or funding agent may find prohibitive. A 
discussion of some construction constraints is provided below. 

11.2.1 Tunnel Size, Shape and Alignment Considerations 

The tunneled excavation, as designed, will satisfy the space demands of the end 
user while providing support to maintain an adequate safety margin against 
tunnel instability. The potential for instability around the tunneled excavation is 
impacted by both the size and shape of the excavation. Where tunnel instability 
is driven by the density and shear strength of natural fractures, the larger the 
opening excavated the greater the likelihood of more frequent, larger rock fall- 
outs occurring and the greater the density and size of the supporting structures 
needed to counter such fall-outs. Where in situ stresses are relatively high 
compared to the strength of the rock mass, the tunnel profile may be modified to 
an elliptical or circular shape. Such cross-sections are selected to avoid corners, 
at which high stress concentrations occur. If in situ stresses are relatively low 
compared to the rock strength, the selection of the cross-section shape may be 
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driven by the economic factors of the mining process itself (excavation, ground 
support and treatment). Water within the rock mass can have a deleterious 
impact on the stability of the excavation. Water pressure acting across planes 
of weakness or flushing through soil- like zones within the rock mass can instigate 
fall-out at the tunnel perimeter. Besides increasing the likelihood of fall-outs, 
the mere presence of water within the tunnel can have a significant detrimental 
impact on the efficiency of the tunneling operation. To reduce the degree of 
interference that water can have on the tunnel work, it is common practice to 
place tunnels on a slight gradient. This gradient provides for gravity drainage 
and collection of water at a sump located cither at the base of a shaft or portal 
and away from the tunnel heading. Slopes are commonly kept below about 3% 
where rail-mounted (steel wheels on steel track) operations are envisaged, but 
can be steeper if only rubber tired vehicle access is required. 

11.2.2 Tunneling Methods and Means 

In all but the weakest rock, three basic types of excavation methods are com- 
monly considered to be feasible for tunnel work. Two of these mining methods 
rely on mechanical breakage of the rock, namely the tunnel boring machine- 
system (TBM-systcm) , and the roadheader; the third method is "drill and blast" 
(D&B) which relies on the use of explosives to effc;c,t rock breakage. The TBM- 
system, which includes not only the cutterhead machine but also the ground 
support and muck evacuation systems, is often used to good effect in the ex- 
cavation of longer, smaller diameter, relatively straight tunnels with uniform 
cross-sectional requirements, mined under more uniform rock mass conditions. 
In such applications, the TBM-system will have the ability to mine tunnel faster 
and cheaper than either of the other two mining methods cited. However, for 
many tunnel projects the TBM is not an automatic choice. It has a relatively 
high capital cost, mines a fixed, circular cross-section, cannot mine tight bends 
or corners and requires an extended period for fabrication/refurbishment, mo- 
bilization and demobilization. When mining shorter, larger and/or more com- 
plex tunnel layouts that demand greater flexibility from the mining equipment, 
including tight turns, steep gradients, multiple cross-sections and variations in 
ground support and treatment en route, roadheader or drill and blast excavation 
methods may be preferred. The roadheader and drill and blast methods have 
similar degrees of mining flexibility, but the roadheader has a limited range of 
economic application. Roadheader viability is severely compromised in harder, 
more abrasive and massive rock masses conditions where mining rate is reduced 
and the abrasive wear on the cutting tools increased significantly. Drill and 
blast tunnel excavation offers the user the most flexible excavation system that 
can be used economically in even the hardest, most abrasive, rocks. However, if 
drill and blast methods are used, the rock mass surrounding the tunnel will be 
subject to additional fracturing due to blast damage. Where explosives are used 
as the means of excavating the tunnel an increment in support should be antic- 
ipated over that required for the mechanically-mined tunnel. Even if the drill 
and blast method is ultimately not selected to mine the tunnel, the excavation 
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of portals and shafts may require the use of explosives. 
11.2.3 Rock Support and Treatment 

In all but the most intrinsically stable rock mass, some ground support will 
be needed in the tunnel. Rock support will be installed to stabilize the tunnel 
periphery and ensure that the miners can work safely within the confines of the 
newly excavated tunnel. The support installed can serve either a temporary 
(during construction) or a permanent (for the life of the project) support role. 
Rock support installed at the heading will be adjusted locally in response to 
variations in the "as-encountered" ground conditions. Rock supports may be 
supplemented by ground treatment work performed around and/or ahead of 
the face where necessary. Treatment of the rock mass may be needed to im- 
prove the tunneling conditions and reduce the impact of construction on the 
surrounding area (for example, water table draw-down or surface settlement). 
Treatment (freezing, grouting etc.) may be undertaken to achieve a temporary 
or permanent increase in rock mass strength or a reduction or increment in rock 
mass permeability. Rock mass zones that may require significant amounts of 
such treatment should be identified early in the site investigation process in or- 
der that their presence, characteristics, extent and mining impacts can be fully 
evahiated during the siting process. Ground treatment, even if only required 
for a short stretch of tunnel along the alignment, can prove time consuming and 
costly. 

11.3 Underground Design Requirements 

At the conceptual stage of design, initial estimates of ground shielding, clear- 
ance envelopes (including tolerances), and general layout/environmental criteria 
required for the construction, installation, operation and maintenance phases of 
the Project should be defined in a drawing set that shows the tunnel(s) in plan 
and section (longitudinal and cross). The tunnel excavations should be laid-out 
to be compatible with selected construction methods and means. The tunnel(s) 
plan and section should show key geologic and hydrologic information. Some un- 
derground requirements that should be estimated during this conceptual period 
include: 

• External loading of floors, wall and crown anchorages including detector 
supports, transportation and lifting system 

• Electrical, electronic, communications networks (cables and cable trays) 

• Heating, cooling, ventilation and air conditioning (duct work, fans, door 
and louvers, drip ceilings and underground chilling/heating units etc.) 

• Groundwater collection and evacuation systems (excavations, drains, pumps 
and pipes) 

• Survey controls (including stations and lines of sight) 
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• Environmental requirements (spoil disposal, groundwater protection.) 

• Neighborhood issues (mitigation of construction/operation impacts on- 
and ofF-site). 

11.4 The Phases of a Tunnel Project 

The main design and construction phases of a tunnel project are outlined in the 
flowchart in Figure I^Hl The flowchart is modeled after the International Tun- 
neling Association guidelines for tunnel design. The flgure indicates a stepwise 
progression from site investigation through to construction and monitoring. In 
practice, the site investigation activity overlaps other planning activities to al- 
low for the detailed investigation of design and the mapping of the excavated 
geology. In addition to the activities listed in the flowchart, there will also be a 
need for estimating and scheduling work. As the project progresses it is likely 
that periodic reviews will be held to evaluate progress, improve management 
confidence in budget and time goals and enhance the practicality and economy 
of the tunnel work itself. 

The phases of the tunnel project are briefly discussed below. 

11.4.1 Site Investigation and Rock Mass Conditions 

To evaluate a sites suitability to be a "host" for an underground facility both 
regional and location-specific geologic information will need to be gathered, in- 
cluding information on rock units, structural folds and faults, groundwater and 
stress regimes. This basic geological information will need to be interpreted to 
characterize the rock mass along the alignment (s) and provide input for the req- 
uisite constructibility and engineering analyses. The interpretation will be used 
to support critical decisions on alignment and selection of methods and means. 
Early acquisition of site investigation data can help quickly identify difficult or 
showstopper situations along an alignment and expedite the short- listing of the 
more serious alignment candidates. Much, if not all, of the site investigation 
data necessary to support initial decisions on alignment and methods and means 
choices can be obtained from desk studies using published reports and papers 
of the regional geology, relevant construction case history data sets and field 
observations at the sites. As the design progresses from the conceptual to the 
alignment-specific stage, site-specific information will be needed to support the 
validation of methods and means for use in design and study rock mass issues 
that were noted as needing further investigation during the conceptual stage 
of the project. At this stage a modicum of alignment specific data will need 
to be acquired typically based on the use of trench and borehole investigation 
and laboratory testing. A general engineering description of the rock mass for 
tunneling purposes will typically include a geologic classification of the rock 
units (ideally with % minerals), an estimate of the intact rock strength, and 
a description of the natural block structure (condition, roughness, orientation, 
size and shape). The potential for the presence of atypical rock mass conditions 
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Figure 39: Flow Chart for Tunnel Design, after Reference |131| 
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also needs to be studied. Atypical conditions that merit investigation include 
soil-like zones within the rock mass, zones of faulting, shears, open fractures, 
solution zones, hydrothermal alteration, weathering and buried valleys. Inves- 
tigative efforts should also be made to evaluate the potential for encountering 
zones of high water inflow that may or may not be associated with soil-like zones 
of weakness within the rock mass. The potential for more pervasive rock-unit 
or regional adverse tunnel conditions, including the presence of relatively high 
in situ stresses, high ambient rock temperatures and more pervasive fluid/gas 
inflows should also be investigated. 

11.4.2 Excavation Methods and Means and Structural Design 

Once a preferred alignment (s) has been identified and basic rock engineering 
characteristics determined, the selection of an initial set of baseline methods and 
means can be made for layout purposes. Throughout the planning period, and 
most notably in support of the selection of methods and means, contractor input 
is highly desirable. Practicing contractors are best positioned to provide state- 
of-the-industry input for selection of safe, practical and cost-effective methods 
and means for tunnel construction. The flowchart in Figure 1391 identifies a 
discrete step for the structural design of a tunnel. This step may be eliminated if 
the rock mechanics or geotechnical engineer considers that a separate structural 
lining (reinforced cast-in-place concrete) is not required. The structural design 
phase may also be skipped or minimized if a similar design case history can be 
referenced. 

11.4.3 Tunnel Contracting and Construction 

Even the most thorough site investigation of the most uniform geologic condi- 
tions will not be able to completely define the scope of an underground con- 
struction contract. Some surprises from the natural material should always be 
anticipated along the way. Risk analyses should be conducted in order to prop- 
erly characterize the likelihood and severity of the impact of all such surprises 
on both the construction work and the project as a whole. Unacceptable risks 
should be mitigated by design, specification, contract provisions or insurance 
measures before the contract is let. The level of risk that tunnel construction 
brings to the overall project can be high, and is strongly influenced by factors, 
including: 

• the complexity of the geology, 

• the thoroughness of the site-specific investigation, 

• the amount and relevance of accumulated case history information, 

• the flexibility of the mining system and, perhaps most importantly of all, 

• by the skill-set of the owner's design and construction team that is assem- 
bled to plan and execute the work. 
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Figure 40: Responsibilities of Contractually-Bound Parties 



Within this context, there is again a strong argument to be made for more 
actively involving the contractors in the design process. Ironically, these key 
protagonists, who have the most relevant construction experience, are commonly 
excluded from the tunnel design. 

Figure 1^ illustrates how, under conventional general contracting practices, 
the responsibility for design and selection of methods and means falls under the 
control of the owner or his representative(s). The desire to involve the contrac- 
tor more actively in the planning of the tunnel project favors the adoption of 
a more integrated design and build approach, where responsibility for design 
of the equipment ("hybrid") or the equipment and tunnel (design and build) 
passes to the contractor. Cording |132| notes "The separation of design and 
specifications from the contractor's planning creates unnecessary impediments 
and adds unnecessary costs to the project." Ultimately, a more integrated de- 
sign strategy that involves the contractor can provide for a more integrated and 
innovative approach to tunneling \1'6'6\ . The use of a design and build contract 
may result in cost and time savings to the Project, but would require the owner 
to freeze requirements at an earlier stage in the development of the project. 
However, the form of contract is a secondary issue compared to the need for 
the owner to assemble and manage a core project team that has a thorough 
understanding of both the end-user needs and tunneling constraints. In the 
case where the owner does not initially have all the requisite skill-set in-house, 
he/she may need to supplement such a team with outside contractors. Adequate 
care should be paid to pre-qualification and selection of such outside contrac- 
tors (consultant-designers and/or builders). The pre-qualification and selection 
processes used by the project team should be project-specific. At a minimum, 
proposing outside contractors should be required to demonstrate a requisite level 
of individual and/or corporate technical expertise and provide work product for 
review that demonstrates the specific qualities that he/she/they can offer to 
the project. Proposals should include a description of recent past experience 
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working on similar jobs (similar requirements, geology, methods and means, 
etc.), and description of levels of project-specific responsibility on the relevant 
projects. References should be provided for follow-up |134) . During design and 
construction the owner's management team should take full responsibility for 
all aspects of project planning and be endowed with adequate responsibility and 
commensurate authority to be able to effectively administer all related design 
and construction activities jl35| . 

11.5 Shaft Hole 

For a reactor experiment, the baseline is already known to be around 1.5 km. 
The far detector should be placed at least several hundred m.w.e. underground 
to suppress cosmic ray muons. If a tunnel is used, the reactor site is to be close 
to steep mountains with a few hundred meters height, depending on the reactor 
power. This condition puts restrictions when selecting a reactor site. The 
situation becomes more problematic for the near detector. The near detector is 
to be located within a few hundred meters from the reactors, but still hundreds 
of mwe of overburden is necessary. If a shaft hole is used, it is possible to 
obtain necessary depth underground at a reactor site which is located on a flat 
landscape. Thus a shaft hole gives added flexibility in the choice of a reactor 
site. 

There are two techniques for digging shaft holes. One is blasting and ex- 
cavating. In this case debris is winched out and dumped outside. Industrial 
needs to dig wide and deep shaft holes are growing. For example, there has 
been intensive R&D to construct 6.5 m diameter and 300 m deep shaft holes as 
access tunnels to underground nuclear waste storage areas. If such a technique 
can be adopted, the shaft holes for our purposes are expected to be constructed 
rather easily and cheaply. 

Another method is drilling. The drilling method is suited for soft soil as the 
case of Kashiwazaki-Kariwa site which is located along a sea shore. The drilling 
is performed by filling the hole with water with higher head pressure than outside 
water level to prevent collapse of the shaft wall. The water helps to make the 
mud-rock softer and to make the drilling easier. The debris generated by the 
drilling are forms of mud and are pushed up to ground level by pressurized 
water, through the center pipe of the drilling machine. When the drilling is 
finished, the water is pumped out and steel and concrete rings are put in the 
hole. There exists a 6 m diameter drilling machine which can dig as deep as 
300 m. 

The capability and cost of the hole construction very much depends on the 
nature of the site and should be evaluated for each case. 

11.6 Summary 

Digging a hole is not as simple as it sounds. Cost and risk are potentially 
much higher than they are for equivalently sized conventional structures built 
on surface or as cut and cover structures using quality-controlled construction 
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materials. Success in design will be largely determined by the owner's ability to 
properly integrate the end-user needs of the facility with the construction needs 
of the tunneler. During construction, the need for good active management 
by the owner will continue in order to ensure that the contract provisions are 
met and, in particular, that the ground conditions as-encountered are properly 
recognized, responded to and documented in an appropriate and timely fashion. 
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12 Safety 



12.1 Safety Planning 

The purpose of this document is to give the present scientific case for a new 
neutrino oscillation experiment at a nuclear reactor. If and when such an ex- 
periment is built and operated, safety must be a paramount concern, and safety 
will be a paramount concern. This will be done by incorporating the following 
principles in the project at every stage: 

• Clear Management Responsibilities 

• Documentation 

• A Working Atmosphere of Safety 

• Integrated Safety Management 

The collaboration and its management will perform and document a haz- 
ard/risk analysis for each phase of the project in order to systematically identify 
the hazards that may be associated with it. This review will be intended to en- 
sure that matters of environmental protection and worker health and safety 
related to the project are identified and that they will be thoroughly addressed 
in the design, construction and operation of the project. Also, at the appropri- 
ate time, a collaboration will be established in accordance with all applicable 
ES&H regulations, standards and good practices. 

12.2 Civil Construction 

The major cost driver for this project is likely to be the construction of a tunnel 
or shaft in order to place the neutrino detector underground. This will require 
substantial civil construction. An outside contractor will be responsible for 
building the underground facility. 

Civil construction of an underground site offers special challenges. Some of 
these involve use of heavy construction equipment, while simultaneously dealing 
with fire protection, confined space issues and the potential for flooding. When 
a specific site is chosen, a set of procedures for instituting safety rules will have 
to be implemented and agreed to by a number of parties. These procedures will 
follow all applicable state laws, federal regulations, and industry practices for 
the selected site. 

In addition, there are a number of special considerations due to the fact that 

this facility will be built in conjunction with a nuclear reactor. Nuclear reactor 
safety in the United States is under the purview of the Nuclear Regulatory 
Commission. It seems clear that nothing involved in the neutrino experiment 
would have any significant consequence on the operation of the nuclear plant. 
This lack of impact may not be immediately obvious to everyone involved, so a 
careful analysis and documentation will be carried out. 



109 



Besides the normal safety requirements involved with the handling of ma- 
terials, there will be special security requirements involved in working near a 
nuclear reactor. A major issue with regard to siting of an experiment is the 
present boundary of the security fence around each facility. These issues are 
already being considered by the security experts who work for nuclear reactor 
power companies as we approach them to discuss possible experiments. In par- 
ticular, there will need to be careful regulation and procedures involving site 
access. These will be worked out in cooperation with all appropriate regulatory 
agencies. 

12.3 Safe Detector Construction and Operations 

As this project develops, there will be a number of technical reviews, cost reviews 
and safety reviews. Safety will be an integral part of documentation at all stages. 
At an early stage of the project, the collaboration will prepare a Preliminary 
Safety Assessment Document which will deal with the the following issues: 

General Construction Safety 

Fire Safety 

Flooding Hazards 

Mechanical Hazards During Installation 
Electrical Hazards During Installation 

Industrial Safety 

Environmental Protection During Construction 

Life Safety - Egress 

Fire Protection 

Electrical Hazards 

Radiological Hazards 

Mechanical Hazards 

Hazardous/Flammable Materials 

Cryogenics and Oxygen Deficiency Hazards 

Emergency Preparedness 

Emergency Communications 

Conduct of Operations 

Training 
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Qualification of Personnel 



Waste Handing; storage and disposal 

Calibration of the detector often involves the handUng of a number of ra- 
dioactive sources. Physicists who participate in experiments at national lab- 
oratories have a great deal of experience regarding the proper training and 
documentation for handling, operating and storing radioactive sources with a 
variety of compositions and strengths. 

Safety is an independent concern and priority for nuclear power generation 
in the United States, under the auspices of the Nuclear Regulatory Commission. 
In the event of any kind of incident at the power plant, experimenters would be 
trained to cooperate fully with plant personnel. 

12.4 Quality Assurance 

It will be the policy of any reactor ^13 project that all activities shall be per- 
formed at a level of quality appropriate to achieving the technical, cost and 
schedule objectives of the project and at the same time insuring that all related 
ES&H considerations are properly addressed. A Quality Assurance plan for the 
project will be developed in accordance with policies and recommendations from 
the Department of Energy and any other appropriate agencies. 
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13 Outreach and Education 

13.1 Goals of the Outreach &c Education Effort 

Outreach is a substantial part of any scientific enterprise. The goals of the ^13 
outreach efforts are: 

1. Providing meaningful education opportunities for the public. 

2. Giving students and teachers of K-12 schools an experience with modern 
science. 

3. Creating and maintaining a learning eommunity between high-school and 
college teaching professionals and the researchers and graduate students 
involved in the experiments. 

4. Developing young people's interests in pursuing science or engineering 
careers. 

5. Increasing participation of members of underrepresented groups in science 
activities. 

6. Emphasizing the impact of pure science on the technological capabilities 

of onr society. 

7. Providing undergraduate and graduate students and K-12 teachers with 
access to valuable research experiences. 

8. Involving local community colleges in the research effort. 

13.2 Outreach and Education Opportunities at the ^13 Re- 
actor Neutrino Experiment 

In order to grasp the outreach opportunities we need to look at our experiment 
from a more general perspective. 

The experiment measures neutrinos generated by a nuclear power plant to 
search for the disappearance of neutrinos as a function of the distance from the 
reactor. The neutrinos disappear as a consequence of oscillations. The flux is 
measured by detecting nuclear reactions initiated by the neutrinos in large vol- 
umes of scintillator oil. Due to the elusive nature of neutrinos the experiment 
has to be pushed to the limits of detection sensitivity. This makes the exper- 
iment also sensitive to background from cosmic rays and other astrophysical 
neutrino sources. 

This makes it clear that the 613 experiment is based on many different dis- 
ciplines, all of which will be a subject of our outreach activities. Our outreach 
efforts are thus primed for interdisciplinary science education which will happen 
within the framework of an international collaboration of scientists. The large 
and exciting list of topics that we can draw from reads as follows: 
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• History of the neutrino The neutrino was postulated by Pauli, its 
experimental discovery was worth a Nobel prize and even John Updike 
wrote a poem about it Particle and nuclear physicists discovered 
that neutrinos have mass and that they oscillate. 

• Particle Physics While much of particle physics is well established, but 
the neutrino sector is not. There are many open questions regarding the 
nature of neutrinos. Are they Dirac or Majorana particles? What are the 
parameters that define the mixing in the neutrino sector? Can we observe 
CP-violation with neutrinos? 

• Nuclear Physics Very rare neutrino-nucleon interactions make the 6*13 
experiment possible. The chain of events that turns the absorption of an 
elusive neutrino into a detectable signal is a fascinating story. 

• Nuclear Engineering Nuclear power generation and the relation be- 
tween neutrino flux and thermal power output broadens the outreach ef- 
forts signiflcantly. Why is it possible to draw so much power from a rela- 
tively small amount of material, measured in kilograms? This is Einstein's 
famous E = mc^ equation put into practice. 

• General Energy Management We see the opportunity to directly com- 
plement existing outreach efforts associated with nuclear energy genera- 
tion. Over one thousand nuclear reactors operate in the world, many of 
which are subject to public scrutiny. Showing that they contribute to the 
understanding of universal fundamental processes could be helpful in the 
public discourse. 

• Astrophysics The biggest producer of neutrinos in our neighborhood is 
of course the sun. Fusion processes in the center of the sun create the 
energy that ultimately fuels all activity on earth. In the fusion process 
neutrinos are generated with a well-known flux of Ve neutrinos. The fact 
that the amount of detected on earth is much smaller than expected 
gave scientists the first evidence that neutrinos oscillate. The oscillations 
were confirmed by studying neutrinos created in interaction of cosmic rays 
with the earth's higher atmosphere. The list of topics is even richer be- 
cause there are also neutrinos expected to come from exotic astrophysical 
objects. Supernovae for example are known to emit a huge burst of neu- 
trinos. 

The outreach efforts we envision will consist of solid but palatable expla- 
nations of each topic, descriptions of their common features, an explanation of 
technology transfer, and generally an invitation to participate, particularly to 
the young, in the science-based activities of the collaborating nations. 

13.3 Strategies for Outreach & Education 

The strategic plan for outreach and education has to be paced with the devel- 
opment of our experiment. We divide the ^13 experiment into three phases: 
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Pre-construction phase Design of the experiment proceeds and the inter- 
national cohaboration forms. First outreach activities are started, such as 
long-term cooperation between scientists and K-12 teachers to stimulate ad- 
ditional joint activities impacting K-12 education. Universities and national 
laboratories are likely places to start outreach efforts by binding 6'13-specific ac- 
tivities into ongoing outreach efforts at those institutions. Quarknet |137| and 
similar activities might serve as examples. To coordinate and further outreach 
efforts an office of outreach will be established within the 6*13 collaboration (see 
also section H 3. 4f) . The first coordinated effort will be to develop a strong online 
representation and to link it to several education databases (e.g. |138) '). The 
experiment will market itself, an effort that will increase throughout the initial 
phase. 

Construction phase With the K-12 teacher-researcher connections estab- 
lished in the first phase and a web presence in place, the second phase of our ex- 
periment will increase our outreach activities by targeting students and teachers 
to become actively involved. The construction of a modern physics experiment 
is a very exciting time and this excitement is likely to be contagious. The en- 
gagement can be achieved by means of research experiences (REU, Quarknet, 
Scientist in Residence) or workshops and seminars. This phase will also be char- 
acterized by a large cooperation with the American Physical Society (APS) and 
American Association of Physics Teachers ( AAPT) . 

Data-taking phase As the detector installation is realized and data taking 
commences, the outreach core content will no doubt be polished. Outreach 
efforts generally will increase. We foresee the possibility to develop an on- 
site welcome center that provides inquiry-driven experiences in content areas in 
which the materials developed for web-delivery are adapted for live engagements. 
The welcome center will be tied into tours of the power plant. We will actively 
work on attracting visitors groups from schools and other institutions. The 
site visit program will be complemented by a program designed to attract more 
minority groups. 

We also intend to focus our attention on aggressively pursuing media oppor- 
tunities in the towns and cities where the researchers/students/staff reside. 

International Activities The ^13 effort is an international one with scientists 
from more than a dozen nations being co-authors of this white paper. The 
outreach efforts have the great opportunity to be internationally oriented, and 
a cooperation of outreach groups from different regions (Europe |139j . Japan 
| 14U| ) will enhance the efforts. 

13.4 Office of Outreach and Education 

The experiment has a large potential to further the education and the interest 
of the general public in science. An office of outreach for the 613 experiment 
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would: 



• develop a rich palette of education materials and establish an strong online 
representation, 

• provide support for the development of the media used in public seminars 
and professional meetings. For example, 3D rendered detector images and 
photographs will be developed and cataloged. Delegating this responsibil- 
ity to the outreach office will minimize redundancy of effort and provide 
a central repository. 

• work with the media for the purpose of education and for marketing the 
experiment, 

• coordinate outreach efforts within the international collaboration 

• coordinate 613 activities with the various nuclear related agencies (e.g. 
Office of Nuclear Energy, Science, and Technology, or U.S. Nuclear Regu- 
latory Committee, DOE) 

• provide professional development opportunities for secondary and post- 
secondary science teachers, and public seminars with the ultimate goal of 
establishing an on-site welcome and outreach center. 

The office of outreach will also work to overcome any confusion that the 613 
experiment is specifically a nuclear physics project. A second challenge will be 
in the concern that any attention to an existing nuclear power plant is unwanted. 

It is important that the ^13 experiment receive regular attention in the press, 
particularly in newspapers with a national audience such as the Washington 
Post, and the New York Times, but also in the local newspapers near where 
research/students/staff reside. The outreach office will operate as a commu- 
nications office in order to "advertise" the program and to provide relevant 
materials for the articles. 

We suggest the outreach office should also "market" the experiment's ed- 
ucation and scientific goals by encouraging and supporting introductory talks 
at professional meetings. Here, we have an enhanced opportunity because of 
the multidisciplinary nature of the experiment. Talks will be coordinated at 
particle, solar, nuclear, astrophysical, and engineering conferences as well as at 
teacher related conferences such as those sponsored by AAPT, the American 
Astronomical Society (AAS) or the American Physical Society |141| . 

13.5 Ideas for Web-Based Outreach Efforts 

Providing rich content is the most important aspect to consider for reaching the 
public with news and descriptions of the 613 project. The initial low-cost/low- 
risk approach to developing content, while maintaining the highest potential 
for impact, is to develop teaching materials for delivery on the Internet. Thus, 
our initial strategy will include a small team of content developers who have an 
ability to integrate technology in both a meaningful and artful way. 
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An obvious choice is to house the development office at a university where 
similar activities take place among students and where a director can have 
ready access to the experiment and its researchers who participate. If selected 
carefuhy, software license costs, which can be significant, could be minimized if 
the host site were aheady participating in Ucense agreements. 

One interesting technology is Physlet IH2| . Physlets are Java-based platform- 
independent applets with physics contents that provide a real interactive expe- 
rience for a modest investment of resources and with the potential to be highly 
visible and accessible to the public. This technology has been prominently fea- 
tured in AAPT workshops and is perfect for outreach activities. 

13.6 Ideas for Working with K-12 institutions 

Physics principles are often regarded as out of reach for many. The art of 
knowing will begin with self-study but, to that end, attention should be given 
to high school teachers who are willing to deliver discourse on this subject in 
their classrooms. Sample lesson plans will need to be constructed. Who better 
to do this but high school teachers themselves? 

Teacher workshops could be organized at the summer AAPT meeting where 
the activities of the ^13 experiment will be the principal theme. However, be- 
cause of its affiliation with the neutrino industry, the four hour workshop will 
include an introduction to the neutrino, an overview of particle/nuclear/astro 
physics, and sample lesson plans that would be useful at both the high-school 
and college level. We expect this activity could be attended by 12-18 teachers 
on an annual basis with only limited costs. 

Another suitable place for representation would be at area planetariums 
where hands on teaching regularly takes place. Constructing a display which 
features the 613 experiment, the adjoining nuclear reactor, and a description of 
the background including astronomy-related backgrounds would draw attention 
to the experiment. 

Direct interaction of K-12 education institutions in the community surround- 
ing the site of the experiment is crucial for the success of the outreach effort. 
Programs like "Scientist in Residence" can help to achieve the goal. Another 
key to a successful outreach program would be to offer research experiences to 
interested K-12 teachers, supplemented by visits of the nuclear power plant and 
the welcome center. 

13.7 Ideas for Outreach to Community Colleges 

Local community colleges usually lack the funding to purchase research equip- 
ment and also lack travel funds to work in remote research facilities. Involving 
community colleges in the vicinity of the reactor site helps to overcome these 
problems and constitute a genuine opportunity for the community colleges to in- 
volve their students in research of international standard. By targeting minority 
colleges, an important goal of science outreach efforts can be met. 
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A Appendix - The Angra reactor in Brazil 



This appendix will focus on the features of the site at Angra dos Reis in Brazil. 
This site has many desirable attributes including good surrounding topography, 
a single powerful reactor, and good relations with the electrical and nuclear 
power companies. In addition, the local physics community has become very 
enthusiastic about this idea and is beginning the process of developing a realistic 
proposal. 

A.l The Reactor Site 

Angra dos Reis is located about 150 km south of Rio de Janeiro. The nuclear 
facility contains two operational reactors. The Angra-I reactor is an older low 
power (about 1.5 GWth) reactor that is not frequently operated. The Angra- 
II reactor, on the other hand, was brought on-line in 2000 and is consistently 
operated at about 4.1 GWth- The reactors are located on the coast and the 
reactor company controls a strip of land that stretches inland about 1-1.5 km 
and is approximately 4 or 5 km along the coast. All experimental constructions 
which will be considered here would be sited within the reactor companies site 
boundaries. 

Much of this terrain is mountainous granite with multiple peaks in the 200- 
600 m region. This allows good background reduction to be achieved for an 
experimental hall with relatively cheap civil construction by tunneling sideways 
into such a mountain. Also within the site boundaries there exists a town, 
Praia Brava, which houses most of the 2 or 3 thousand people which work at 
the reactor facility and also contains a hotel and stores. Such already existing 
infrastructure could make using this facility more attractive. 

A. 2 Communication with the Power Companies 

The company which runs Angra is state owned and operated. One of the unique 
features of attempting this experiment in Brazil is that the presidents of both 
the electric power company and its daughter nuclear power company are former 
particle physicists who used to do experiments at CERN. As a result, they 
are very receptive to communications from members of the Brazilian physics 
community and have been very helpful in providing resources and access to the 
facility. A one day site visit has already been performed to evaluate the viability 
of performing the experiment there. Significant assistance was provided by the 
director of operations from the reactor facility and significant time was spent 
with the director of civil construction on the site. With their help, possible 
experimental site locations were explored. The reactor company has agreed to 
supply full detailed cost estimates of any civil construction plan by using their 
knowledge of the site geology and known contractors. 
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Figure 41: A topographic map of the nuclear reactor site at Angra dos Reis. 
The concentric circles are at 500 meter radial intervals from the core of Angra- 
II. Proposed locations for the near and far detector experimental halls as well 
as the far detector access tunnel are shown. 

A. 3 The Experimental Design 

A topographic map of the site is shown in Figure 1411 The concentric circles 
are at 500 m radial intervals from the primary Angra-II reactor core. The near 
site location is 300-350 m from the core. It has the possibility to gain about 
15-20 m of rock overburden (30-50 meters of water equivalent). The far site 
location would exist under a 240 m hill at about 1350 m from the reactor core. 
Access to the far site would come from a 420 m tunnel which starts from the 
western edge of the hillside. This location is easily accessible from the town 
of Praia Brava and would be very near to the current location of their sewage 
treatment plant. 

It is envisaged to place identical 50ton fiducial detectors at each location. 
Exact detector designs have not yet been developed, but it is currently assumed 
that such detectors would build off of the developments from other groups. 
Most likely a 3 volume detector would be optimal: a central liquid scintillator 
volume that would be doped with gadolinium (target); a surrounding volume 
of liquid scintillator without gadolinium (gamma catcher); a non-scintillating 
buffer to shield the radioactivity of the photo-tubes which would be installed at 
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the outer edge of this volume. An active muon shield would then be required 
to surround this system. A spherical detector with an active target of 50 tons 
would have a total diameter of approximately 7.3 meters. The access tunnels 
and experimental halls are being designed to accommodate these dimensions. 

A. 4 Experimental Reach 

Preliminary estimations have been performed of the signal and background rates 
for the given detector configuration. The detector at the far location is expected 
to get about 120 signal events per day, while the near site would be expected 
to receive about 3000. Some very preliminary background estimates suggest 
that at the far detector, less than 10 Hz uncorrelated backgrounds, which would 
easily be vetoed by an active muon shield, would be expected and there would 
be about 1-2 correlated background events per day from muon induced radioac- 
tive isotopes. Similarly at the near detector there would be an uncorrelated 
background rate of about 830 Hz (yielding an active live time of 63% after 
muon vetoing) and a correlated background of approximately 150 events per 
day. Having a signal to noise rate of about 100 in the far detector and 20 in the 
near detector should allow reasonable background rejection while maintaining 
statistical sensitivity. Figure shows the expected statistical sensitivity as a 
function of time, for the best fit value and 90% allowed limits of Ato^ from 
Super-Kamiokande. As can be seen, a limit of sin'^{26i^) < 0.02 at 90% con- 
fidence level can be achieved within 3 years. Also in Figure 021 is shown the 
complete limit and 3a discovery potential for a 3 year run over all sm^(20i3) 
and Am^. 

A. 5 Brazilian Community and Support 

The Brazilian community has recently been having in-depth discussions about 
this possibility. There exists significant theoretical and phenomenological sup- 
port for neutrino oscillation work in Brazil already, primarily located in the 
Universidade Estadual de Campinas (UNICAMP) and the Univcrsidade do Sao 
Paulo. The director of the Centro Brasileiro de Pesquisas Fsicas (CBPF) in Rio 
de Janeiro has shown great enthusiasm for this project and has offered to make 
his institution the host of such an experiment. There exist strong experimental 
particle physics groups based mainly in Rio de Janeiro and Sao Paulo. They are 
already strongly involved with the Pierre Auger project as well as experiments 
based at Fermilab and CERN. They have expressed great interest in pursuing 
a project that is local to Brazil and several of the experimental members have 
already started working on producing a more realistic set of cost estimates and 
site plans. 

While it's not possible to put any real faith behind the numbers, some pre- 
liminary estimates have concluded that with the cheaper labor costs in Brazil, 
the civil construction could probably be achieved for 5-7 million dollars (US). 
The local Brazilian community is currently in consultation with the govern- 
ment to get a commitment to this project and for covering the complete costs 
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Figure 42: Expected sensitivity to sin^ (26*13) which could be achieved by an 
experiment at Angra dos Reis assuming a 1% systematic error. The plot on the 
left shows the sensitivity as a function of years of running for three different 
values of Am^. On the right, the full coverage of Am^ vs. sin^(20i3) is shown 
assuming three years of data taking. Curves for both the limit at 90% confidence 
level and the discovery at 3 cr are shown. The current limit at 90% confidence 
is sin2(26li3) < 0.2. 

of civil construction. For the detectors, a very conservative cost estimate con- 
cluded that they would require about 7 million dollars each. Thus the total 
project could probably be completed with an external contribution of less than 
15 million dollars (US). 
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B The Double Chooz Project 



Since the end of the year of 2002 a group of European (French, German, and 
Russian) physicists has been working on the possibility to measure the last 
undetermined neutrino mixing angle 613 with a two detector reactor neutrino 
experiment. During the year 2003, a site investigation has been performed, and 
three sites have been sclcctcid to be potentially suitable to performed such a 
project: the Penly site (Normandie, France), the Cruas site (Ardeche, France), 
and the Chooz site (Ardennes, France). Thanks to its existing infrastructure 
and to the support of the EDF (Electricite De France) power plant staff, Chooz 
is the most reliable location to perform an experiment and we'll only focus on 
it in the following. 

B.l The Double-CHOOZ concept 

The new Chooz project will run two identical detectors of medium size, con- 
taining 10 to 20 tons of liquid scintillator target. The (still existing) neutrino 

laboratory of the first Chooz experiment, located at 1.05 km from the two cores 
of the Chooz powerful nuclear plant (PWR, 8.4 thermal GW) will be used again. 
We label this detector the far detector or CHOOZ-far. The CHOOZ-far site is 
shielded by about 300 mwc of 2.8 g/cm"^ rocks. We plan to start to take data 
at CHOOZ-far within the first month of the year 2007. 

In order to cancel the systematic errors originating from the nuclear reactor 
(electron anti- neutrino flux and spectrum), as well as to reduce the set of sys- 
tematic errors related to the detector and to the event selection procedure, a 
second detector will be installed close to the nuclear cores, within a range from 
100 to 200 meters. We labc;! this detector the near detector or CHOOZ-near. 
Since no natural hills or underground cavity exists at this location, an artifi- 
cial overburden of a few tens of meter height has to be built. Prom the first 
discussions, this construction has been allowed by the power plant company au- 
thorities. A pre-study of the civil engineering work to be done for CHOOZ-near 
has been order and will be performed during the end of 2003. This pre-study 
which aims to determine the best combination location-overburden as well as 
to optimize the cost is financed by the French electricity power company EDF. 

From the simulation, a new experiment at CHOOZ would reach a sensi- 
tivity of sin^ (26*13) < 0.03, at 90% CL for Am^^^ = 2.0 x lO^^ eV^ (latest 
best fit value of the Super-Kamiokande experiment), after three years of oper- 
ation. A sensitivity of 0.05 would be reached within the first year of running 
with two detectors, improving the CHOOZ boimd by roughly a factor of flve. 
These estimates are based on the assumptions that the relative normalization 
error between the near and far detectors could be kept at 0.8%, and that the 
backgrounds at both sites amount to less than 1 % of the anti-neutrino signals. 

It is worth mentioning that the 1.05 km average baseline at CHOOZ is not 
optimal (a little too short) compared to the first oscillation maximum (for anti- 
neutrino energy of a few MeV) if the atmospheric mass splitting would remain 
at Ato^j^ = 2.0 X 10~^ eV^. A shorter baseline is compensated by higher 
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statistics for a fixed size detector, however. A value of Ato^^^ < 1.5x 10""^ el^^, 
close to the 90 % CL current lower bound, would restrict the absolute potential 
of the double-CHOOZ experiment. 

B.2 Detector design 

The detector design foreseen is an evolution of the detector of the first ex- 
periment, (CHOOZ-I). The first improvement with respect to the CHOOZ-I 
sensitivity comes from an increase of the exposure to get more than 40,000 anti- 
neutrino events at CHOOZ-far. This condition translates to a statistical error 
lower than 0.5%. This increase of signal is obtained by using a target cylinder 
of 120 cm radius and 280 cm height, providing a volume of lAm^, 2.5 higher 
than in CHOOZ-I. In addition, the operation period will be extended to at least 
three years, and the efficiency will be improved (the global load factor of the 
reactor is about 80 % while it was much lower for the CHOOZ-I experiment due 
to the power plant commissioning). In addition, the backgroimd level will be 
decreased in order to get a signal to noise ratio greater than 100 (about 30 in 
CHOOZ-I). 

The near and far detectors will be identical up to the PMTs surface. This 
will allow to have a relative normalization error less than 0.8%. However the 
outer shielding will not be identical since the cosmic background varies be- 
tween CHOOZ-near and CHOOZ-far, due to the different overburdens. The 
overburden of the near detector has been chosen in order to keep the signal to 
background ratio above 100; with this condition, even a poor measurement of 
the background keeps the associated systematics error well below 1 %. 

The detector design was intensively studied and tested by Monte-Carlo sim- 
ulation, using a code derived from the simulation of the CHOOZ-I experiment. 
The one meter low radioactive sand shielding of CHOOZ-I, used to reduce the 
external radioactivity contamination, will be replaced by a 15 cm stainless steel 
shielding. The size of the liquid active buffer can then be increased. Prom the 
center there will be: 

• a 120 cm radius 280 cm height acrylic target cylinder, filled with 0.1 % Gd 

loaded liquid scintillator 

• a 40 cm buffer of non-loaded liquid scintillator with the same optical prop- 
erties (light yield, attenuation length), to get the full positron energy as 
well as most of the neutron energy released after the neutron capture (this 
buffer is called the 7 catcher) 

• a 95 cm buffer of non scintillating liquid, to decrease the level of the 
accidental background 

• a 60 cm veto filled with liquid scintillator. 

The spatial reconstruction is not affected by the cylindrical shape. Each 
parameter of the detector is being studied by Monte-Carlo simulation, to define 
the tolerance on the differences between the two detectors. The inner volume 
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dimensions as well as the shape of the target vessels are still preliminary and 
could change prior to the completion of the proposal. The aim is to reduce 
the number of analysis cuts with respect to the CHOOZ-I experiment to push 
down the systematic errors related to the anti-neutrino event selection in both 
detectors independently. The non-scintillating buffer will reduce the singles 
rates in each detector by two orders of magnitude with respect to CHOOZ-I. 
This lowers the positron threshold down to 500 keV, well below the 1 MeV 
physical threshold of the inverse beta decay reaction. A very low threshold 
has three advantages: a) the systematic error due to this threshold is then 
suppressed, b) the background below the threshold can be measured, c) the 
beginning of the positron spectrum provides an additional calibration point 
between the near and far detectors. This reduction allows the relaxation or 
even suppression of the localization cuts, such as the distance of an event to the 
PMT surface and the distance between the positron and the neutron. These 
cuts, used in CHOOZ-I, are difficult to calibrate and have to be avoided in 
a precision experiment. The remaining event selection cuts will have to be 
calibrated between the two detectors with a very high precision. The main cut 
to calibrate is the energy selection of the delayed neutron after its capture on a 
Gd nucleus (with a mean energy release of 8 MeV gammas) . The requirement is 
50-100 keV on the precision of this cut between both detectors, which is feasible 
with standard techniques using radioactive sources at different positions and 
lasers. The main sensitivity of a reactor experiment of the Double-CHOOZ 
scale (400 GW t y) is the total number of events; the requirements on the 
positron energy scale are then less stringent since the weight of the spectrum 
distortion is low in the analysis. This is being studied by simulation. 

B.3 Backgrounds 

The correlated background was measured in CHOOZ-I, and found to be about 
one event per day. At CHOOZ-far the active buffer will be increased by a factor 
of two, with a solid angle for the out-coming background being unchanged. This 
together with a signal increased by a factor of 2.5 will fulfill the requirement of 
a signal to noise ratio greater than 100. 

At CHOOZ-near, due to the shallow depth, the cosmic background will 
be important. For a detector located at 150 meters from the cores, with an 
overburden of 55-65 mwe, the signal will be 5,000 events per day, while the 
muon rate is expected to be a factor of ten less. A dead time of 500 fxs will be 
applied after each muon, leading to a global dead time of about 25 %. About 20 
recoil proton events per day, mimicking the anti-neutrino signal, are expected 
while the estimate of the muon induced cosmogenic events {^Li and ^He) is 15 
per day with a large uncertainty (this last point is being studied). This fulfills 
the requirement of a signal to noise ratio greater than 100 at CHOOZ-near. 
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B.4 Conclusion and outlook 



A new reactor neutrino experiment at the Chooz nuclear power plant has the 
potential to reach a sensitivity of sin^(20i3) < 0.03, at 90% CL for Am^j^ = 
2.0 X 10~^ eV'^, after three years of operation. This potential is similar to what 
could be achieved by future long baseline neutrino experiments (JPARC and 
NuMI-Off axis) at the horizon 2013. Furthermore, both results are complemen- 
tary. 

The CHOOZ-far detector will be installed in the fall of 2006, in order to start 
data taking in the beginning of 2007. Due to civil engineering constraints, the 
CHOOZ-near detector will be installed one year after CHOOZ-far, and will start 
data taking beginning 2008. The Double-CHOOZ proposal will be written in 
the forthcoming months by a team of European physicists from three countries 
and five laboratories (about 20 physicists); the international collaboration is 
expected to grow during this period. 

We would like to thank the Electricite de France power company for its 
technical and financial support to the Double-CHOOZ project. 
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C Daya Bay 



The Daya Bay nuclear power plant is the largest nuclear power generating sta- 
tion in China, and the tenth largest in the world. It consists of two twin reactor 
cores, one is called Daya as shown in Figure03l and the other is called Lingao as 
shown in Figure^] each core can generate a thermal power of 2.9 GW, making 
a total of 11.6 GW. All reactors are the type of the pressurized water. A third 
twin-core with a total thermal power of 5.8 GW, is planed to be online at about 
2010. It will then be the second largest nuclear power plant in the world, only 
next to the Kashiwazaki in Japan. 

Daya Bay is located about 50 km from ShenZhen, Guangdong province, and 
it is only two hours drive from Hong Kong. The reactors are build on coast 
next to a mountain whose highest point is 700 m above the sea level. A group 
of Chinese physicist has inspected the site and it seems no problem to have a 
near experimental hall with an overburden of 400 mwe at a distance of about 
300 m to the core, and a far experimental hall with an overburden of 1200 mwe 
at a distance of about 1500-2000 m to the core. Initial discussions about the 
possible experiment with the authority of the power plant are positive and a 
conceptual design of the tunneling and experimental hall in collaboration with 
the power plant is underway. 

A group of physicists has been working on the this possibility, including ne- 
gotiation with the power plant, conceptual design of the experimental hall and 
tunneling, detector design and R&D, Monte Carlo simulation etc.,. Collabora- 
tion with all interested parts are welcome. 
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D The Diablo Canyon Power Plant 



A number of commercial nuclear power plants worldwide are currently under 
evaluation as sites for a next-generation reactor neutrino oscillation experiment. 
An overview of these activities is given in this paper and references herein. One 
of the candidate sites in the United States is the Diablo Canyon Power Plant 
(DCPP) on the Central Coast in California. The Diablo Canyon Power Plant 
is located in San Luis Obispo County north of Santa Barbara and provides 
electricity to northern and central California. DCPP is owned by Pacific Gas 
and Electric Company (PG&E) [HS- 

A powerful reactor and overburden in excess of 300 mwe to shield the an- 
tineutrino detectors from cosmic rays and associated backgrounds are principal 
features of a suitable site for a reactor neutrino experiment. An underground 
detector facility or artificial overburden are usually required to achieve sufficient 
shielding from cosmic rays. 

The Diablo Canyon Power Plant consists of two 1.1 MW reactor cores sep- 
arated by ^ 100 m with a total thermal energy of 6.2 GW th- Nearby coastal 
mountains provide the opportunity for good overburden and make the plant 
an almost ideal site for a reactor neutrino experiment |144| . The excavation of 
a horizontal tunnel in the coastal mountains can provide overburden of up to 
800 mwc with distances up to ^ 3 km. The general layout and topography of 
the site allows the construction of a 1-2 km-long tunnel for two or more movable 
detectors. A picture of the Diablo Canyon Power Plant is shown in Figure B51 
and Figure 1461 gives a topographic map of the area surrounding the DCPP site. 

Construction of a horizontal tunnel is required to place the neutrino detectors 
underground and shield them from cosmic rays and associated backgrounds. 
The local topography allows the construction of a tunnel at distances between 
0.9-2 km from the reactor. Longer tunnels up to a distance of 3 km are possible. 
A near detector within 500 m of the reactor cores requires the construction of 
a detector room and artificial overburden in excess of at least 50 mwe. 

The neutrino detectors can be placed in detector rooms located in side drifts 
off the main tunnel (see Figure HTjl . This allows the detectors to be surrounded 
by a hermetic muon veto and passive shielding reducing the cosmic-ray induced 
neutron and external 7-background from the rock. The side drifts can provide 
a clean-room environment for the sensitive neutrino detectors and air control 
for the reduction of radon. Each side drift will have enough space to accom- 
modate two neutrino detectors in order to allow for the relative calibration of 
the detectors' efficiencies at one particular distance from the reactor and under 
the same background conditions. Due to the natural topography of the site the 
overburden and hence the cosmic ray background will vary depending on the 
detector location and the distance from the reactor. 

A horizontal access tunnel will allow for the possibility to optimize the base- 
line of the detectors with respect to the reactor. The exact location of the side 
drifts can be chosen when the parameters that determine the neutrino oscilla- 
tion length (such as Am^^^) and hence the baseline of the experiment are better 
known. With multiple side drifts for the detectors one can imagine demonstrat- 
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Figure 45: Picture of the Diablo Canyon nuclear power plant in San Luis Obispo 
County, California, USA. The local topography at Diablo Canyon allows the 
construction of an underground tunnel between 0.9-2 km for the placement of 
two or more neutrino detectors. A longer tunnel up to a distance of 3 km is 
possible. It may be possible to place an additional near detector at 0.4 km under 
artificial overburden. 




Figure 46: Topographic map of the site of the Diablo Canyon Power Plant. 
The land boundary (black) as well as the power plant site boundary (red) are 
indicated. 
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Figure 47: Horizontal access tunnel and side drifts for the the detector rooms. 
The side drifts will accommodate two neutrino detectors. (Drawings provided 
by D. Oshatz, LBNL.) 
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ing conclusively the subdominant oscillation effect with measurements at various 
locations along the tunnel. 

Extensive geological studies of the site of the Diablo Canyon Power Plant 
have been done by PG&E. A preliminary evaluation of the geology for this 
project is in progress at Berkeley Lab. The geology of the site appears to be 
suitable for tunneling and the excavation of the underground detector rooms. 
A tunnel may even offer the opportunity for interesting geo-science studies. 
The stratigraphy of the Diablo Canyon region is dominated by three Pliocene- 
Miocene marine sedimentary units: the Pismo Formation, the Monterey For- 
mation and the Obispo Formation. The rocks present in the area around the 
prospective tunnel sites at Diablo Canyon are predominantly the older sedimen- 
tary and volcanic rocks of the Obispo Formation, consisting of (1) tuffaceous 
siltstones and claystones, and (2) basaltic flows, dikes, and sills. There are 
a variety of geologic hazards and issues associated with the construction and 
maintenance of a neutrino detector facility at Diablo Canyon: landslide haz- 
ards, rock quality issues, seismic hazards, and water quality. These issues are 
currently under investigation. The average bulk density of these rock formations 
varies between 2.2-2.8 g/cm^. This leads to a maximum effective overburden 
of ~ 800 mwe at various locations along the tunnel. Geologic cross-sections of 
four possible tunnels at the DCPP site are shown in Figure HSl These tunnel 
cross-sections start at an elevation of ^ 85 m at a location approximately 0.9 km 
East-North-East from the Diablo Canyon Power Plant. 

The physics potential of a next-generation reactor neutrino measurement is 
discussed in great detail in the previous sections of this paper. The ultimate 
sensitivity of such an experiment will strongly depend on the layout of the 
experiment including the distances from the reactor cores, the overburden, the 
size of the detectors, and our understanding of the relative detector response. 
The Diablo Canyon Power Plant is a site that offers the opportunity to design 
an experiment that meets these criteria and achieves an ultimate sensitivity of 
sin2(26'i3) = 0.01-0.02. 
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Figure 48: Geologic cross-sections of four possible tunnels at the DCPP site. 
The distance is given in meters from the tunnel portal. The tunnel portal is 
located at a distance of about 0.9 km from the reactor. A maximum overburden 
of up to 300 m of rock or 800 mwe can be obtained. (Figure provided by 
C. Onishi and P. Dobson, LBNL.) 
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Figure 49: Interaction of the Muon Veto Neutron Shield system with muon 
induced fast neutrons produced both inside and outside the shielding bunker. 




E An Illinois Reactor Experiment 

The Illinois proposal is to locate the experiment at one of the Exelon nuclear 
reactors: Braidwood, Byron, or La Salle. All three of these reactors are located 
in northern Illinois with 100 km of both Argonne National Lab and Fermilab. 
The typical elevation of northern Illinois is about 250 meters and the topology is 
flat. The baseline experiment design calls for three identical spherical detectors, 
each with a 25 ton gadolinium loaded scintillating target. There is one near 
detector located at a distance of about 200 meters from the reactor cores. The 
remaining two detectors are at a far location and share a baseline of 1400 to 
1800 meters. Both the near and the far detector halls are located below > 
100 meters of rock (~ 300 mwe). This shielding is achieved by digging shafts 
straight down to a depth of about 100 meters 

In addition to the standard source and flasher calibration systems, the rel- 
ative efficiency of the near and far detectors is measured head-to-head. To 
facilitate this measurement, the detectors are movable, and the near and far de- 
tector sites are connected by a tunnel (See Figure EH} ■ The relative efficiency is 
measured with the two detectors side by side in the intense neutrino flux at the 
near detector site. Calibration running for 10% of the experimental run results 
in an uncertainty on the relative efficiency that is smaller than the statistical 
error on the event rate in the far detector. 

Correlated coincident backgrounds, in which both parts of the event are 
generated by the same initial interaction, are the most difficult to account for. 
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Most correlated backgrounds are due to fast neutrons generated in cosmic ray 
^ spallation in the materials surrounding the detector. Inside the detector a 
fast neutron can strike a hydrogen nucleus giving it enough energy to mimic the 
positron. Alternatively, a neutron-carbon inelastic collision may result in the 
production of gamma rays. In both cases the fast neutron will thermalize and 
capture with nearly the same spatial and time correlation as the neutrons from 
inverse /3-decay. 

To deal with these correlated background events a Muon Veto Neutron Shield 
(MVNS) system is proposed. In the MVNS system (shown in Figure . the 
detector is housed in a bunker of dense material intended to range out neu- 
trons. The outside of this bunker is covered in an array of plastic scintillator. 
The scintillator array detects muons entering the bunker. Muons in the bunker 
material may kick out fast neutrons and cause a correlated background event in 
the detector. If a muon passes close to the bunker without passing through, neu- 
trons generated by it in the surrounding rock are ranged out by the bunker. In 
addition, matching the energy distribution for muon events to the neutrino can- 
didate events, outside the reactor energy range, the any surviving background 
events can be subtracted from the neutrino event sample. 

The MVNS bunker is fixed and does not move with the detector. 

The Illinois proposal is designed to control the dominant systematic errors, 
which are due to the relative efficiency of the near and far detectors and the 
uncertainty in the background rate. By keeping these systematic errors low, a 
sensitivity of sin^ 26*13 < 0.01 at 90% CL can be achieved in a three year run. 
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F The KASKA project 



In the KASKA experiment |145| . three identical detectors are to be buih in the 
site of Kashiwazaki-Kariwa NPP which has 7 reactors, producing total thermal 
energy of 24.3 GW. This is the most powerful NPP in the world. Using a large- 
power nuclear power plant is profitable for obtaining not only a high event rate 
but also a low background-to-signal ratio at a given depth underground. The 
relative locations of reactors and detectors are shown in the fig.-ISUI Although 
the far/near distance ratios between the reactors and detectors are not unique, 
the effective uncertainty introduced from the variations of the distances are 
estimated to be only 0.2%. 



Reactor and Detector Locations 
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Figure 50: Reactor (circles) and detector (squares) relative locations. Reactor 
#1 through #4 form a cluster and #5 through #7 form another cluster. The 
two clusters separate about 1.3 km apart. Two near detectors will be placed 
at around 400 m from each cluster. The far detector will be placed at around 
1.3 km from all the reactors. 



The detector is a CHOOZ like detector as shown in the Figure IFTl 



The Ve target is 8 tons of Gadolinium loaded liquid scintillator (GD-LS) con- 
tained in a acrylic vessel which is called region-I. The component of the Gd-LS 
is the Palo Verde type, which is proven to be stable in the acrylic container |75| . 

The Gd concentration is 0.1% (0.15% optional) and the neutron capture 
efficiency on Gd is 88%. The optional 0.15% Gd concentration is intended to 
increase the neutron capture efficiency and to reduce the systematic uncertain- 
ties associating with the inefficiency. The reactor Ve is detected by the inverse 
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Figure 51: Schematic view of the detector. 
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/3 decay reaction with proton. 

i7e+])^e+ + n (26) 

The kinetic energy of the positron is neutrino energy minus 1.8 MeV. The 
positron annihilates with electron after slowing down in the LS and then pro- 
duces two 0.511 MeV 7's. This process produces a prompt signal, whose energy 
is between 1 MeV and 8 MeV. 

The energy threshold for prompt signal is set to be below the minimum 
prompt signal energy of 1 MeV. In this way, no systematic ambiguities associated 
with threshold energy cut is introduced. 

The neutron produced in the inverse (3 decay reaction thermalizcs quickly 
and is absorbed by Gd, producing 7-rays whoso total energy amounts to 8 MeV. 
The neutron absorption signal occurs typically 30 /iS after the prompt signal. 

The neutrino signal is defined by the existence of correlated signals. That 
is, when 8 MeV of energy deposit is observed after associating prompt signal 
whose energy is between 1 MeV and 8 MeV, this event is considered to be 

event, regardless the positions of the signals. Becaiise there is no position 
cut, the measured neutrino rate is free from the position reconstruction error. 
Requirements for the timing correlation between the positron signal and the 
neutron signal suppress backgrounds severely 

The main source of systematic error comes from the relative difference of the 
LS mass in the three detectors. The total volume of the LS will be controlled 
not by the size of the container but by the liquid volume put in the container. 
The relative volume of the introduced LS will be calibrated before putting in 
each detector. The total volume of the LS in the oval part of the vessel can be 
measured precisely from the total volume of introduced LS and the liquid level 
in the thin calibration pipe even if there is a distortion of the vessel after the 
installation. The relative uncertainty of the LS mass between detectors will be 
less than 0.5%. 

If neutrino events occur near the boundary of region-I, 7 rays from the 
positron annihilation and the neutron absorption may escape from region-I. In 
such cases, the visible energies of the signals may become less than thresholds 
and this inefficiency introduces potential systematic error. This is also true for 
the case of positron detection. 

In order to catch such 7 rays, a 50 cm thick Gd unloaded LS (region-II) 
surrounds region-I acrylic container. The light output of region-II LS is adjusted 
to be same as region-I LS to obtain accurate original energies accurately. It is 
estimated that most of the 7 rays can be contained within regions I and II and 
the inefficiencies of the positron and neutron signal is less than a few percent and 
the systematic error associated with this inefficiency is a fraction of a percent. 

If the neutrino event occurs near the acyclic wall, the neutron may escape 
from region-I (spill-out effect). This kind of events produces a loss of efficiency. 
However, the inefficiency due to the spill-out effect is estimated to be only ~3% 
and the systematic error due to the spill-out effect is estimated to be on the 
order of 0.1%. On the other hand, there is the opposite case that the inverse 
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P decay event occurs in region-II and then the resulting neutron spills into 
region-I (spill-in effect) and absorbed by Gd. This event is also considered to 
be a neutrino event and it amounts to 20% of the total events. However, the 
spill-in efficiency is insensitive to the hydrogen concentration in region-II LS 
because the spill-in event rate is proportional to the neutrino event rate and 
the neutron absorption length and both parameters have opposite dependence 
on the hydrogen concentration. Because most of the events are contained in 
regions I and II, systematic error due to the spill-in effect is estimated to be less 
than 0.5%. 

The outer wall of region -II is also a acrylic container and the region-II is 
surrounded by a 90 cm thick scintillator with very slight light output (region- 
Ill). Region-Ill shields 7-rays from PMT glasses. Thanks to the region-Ill 
buffer, the singles rate is expected to become less than 10 Hz. 

The outer most layer (Region-IV) is muon anti-counter filled with the same 
scintillator as region-Ill. These layers work as a shield of gamma rays and as 
cosmic ray anti counters. The slight light output is to detect low energy muons 
whose velocity is below the Cerenkov threshold, and to obtain high muon tagging 
efficiency. There will be about 400 8-inch low background PMTs at the boarder 
of region-Ill and IV. The wall between region-Ill and IV is just for light shield 
and no hermeticity is necessary. 

There are cosmic-ray tracking devices on both the top and bottom of region- 
IV. The trackers are used to measure cosmic-ray tracks with accuracy of ~10 cm. 
The cosmic-ray spallation background rate, such as the ^Li rate, is estimated 
from the excess at small distances in the distribution of distance between the 
neutrino event candidates and the cosmic-ray tracks. In order to stop the 
hadronic and soft part of the direct cosmic rays, the thickness of the top part 
of region-IV will be made thick. Region-IV LS is contained in a stainless steel 
cylindrical hermetic container. Outside the stainless steel container, there is a 
thick iron shield to prevent 7-rays coming from the soil outside. 

The far detector will be placed at the bottom of a 200 m shaft hole with an 
inner diameter of 6 m. The cosmic-ray flux at the bottom is about 0.35 m~^s~^. 
Although there is open space above the detector, the cosmic-ray rate directly 
down the open space is only a fraction of that which reaches the detector pene- 
trating the soil. The near detectors will be placed at the bottom of a 70 m depth 
shaft hole. This depth is chosen to make muon/neutrino ratio approximately 
equal to or slightly larger than the far detector case. 

The major component of the background comes from fast neutrons produced 
in nuclear interactions caused by cosmic rays going through the rock near by. 
The visible energy distribution of the prompt signal in the fast neutron back- 
grounds was measured to be flat by CHOOZ group at the energy range below 
30 MeV |10j and this kind of background rate can be estimated by using the 
event rate within non-reactor-Pg energy range, such as below 1 MeV and above 
10 MeV. 

The absolute background rate is expected to be less than 2% and the error 
associated with estimation of the background will be less than 0.3%. 
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The systematic error in CHOOZ experiment was 1.7% (detector associated) 
® 2.1% (neutrino flux associated). The systematic error associated with the 
neutrino flux is reduced to 0.2% by comparing near and far detectors in the 
KASKA experiment. The systematic error associated with the detector is esti- 
mated to be between 0.5% and 1%. Because the neutrino osciflation at the near 
detectors is smafl, the consistency check of the systematic error estimation can 
be performed by comparing data of the two near detectors. 

In two years of operation, 40,000 neutrino events wiU be recorded in the 
far detector and ten times more in the near detectors and the corresponding 
statistical error wiU be 0.5%. The 90% CL sensitivity of this experiment is 
shown in the flg.-El At Am^ - 3 x 10'^ eV'^, the sensitivity of sin^ 26*13 ~ 0.02 
is expected. This is five times better limit than CHOOZ. 
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Figure 52: The expected 90%CL exclusion region of this experiment for the case 
of asys = i% and 0.5% obtained by rate only analysis. 



140 



G The Krasnoyarsk Reactor and KR2DET 



The KR2DET proposal is to place two identical liquid scintillation spectrome- 
ters are stationed at distances Lfar ~ 1000 m (far position) and Lnear ~ 115 
m from the underground Krasnoyarsk reactor. (Figure The overburden at 
Krasnoyarsk is ~ 600 m.w.e., which is twice as much as in the CHOOZ experi- 
ment. (At short distances form the reactor the one reactor - 2 detector approach 
was first probed at Rovno ^146; and later successfully used at Bugey |147| . 

Two types of analysis can be used. Analysis I is based on comparison of the 
shapes of positron spectra S{Ee)far and S{Ei,)near measured simultaneously in 
two detectors. In no oscillation case the ratio S{Ee) far/ S{Ei,)near is energy 
independent. Small deviations from the constant value of this ratio 

1 — sm 20 ■ sm I g — ^ J 

shave ~~ ^ ' 7 \~ (^"7) 

l-sin2 20.sin2(M7Ar^ ^ 

are searched for oscillation parameters. 

In the one reactor - two detector scheme 

• Results of the Analysis I do not depend on the exact knowledge of the 
reactor power, absolute v,, flux and energy spectrum, burn up effects, 
absolute values of hydrogen atom concentrations, detection efficiencies, 
target volumes and reactor - detector distances. 

• At Krasnoyarsk the detector backgrounds can be measured during reactor 
OFF periods, which periodically follow 50 day long reactor ON periods. 

Calculated ratios S{Ee) far/ S{Ee)near for a set of oscillation parameters are 
shown in Figure 15^ 

Analysis II is based on the ratio of the total number of neutrinos Nfar, Nfar 
detected at two distances: 

Xrateism' 29, Am'} = (^) ' . (^) ■ f ^) . (P^) (28) 

\^ near J \ ^far / \ ^far / near J 

Vfar, Vnear, ^far, ^near are the target volumcs and neutrino detection efficiencies. 
In no oscillation case Xrate — 1- 

Analysis II is also independent of the exact knowledge of the reactor neutrino 
flux and energy spectrum. The absolute values of detection efficiencies are 
practically canceled, only their small difference is to be considered here while 
the ratios {L far/ LnearY ^nd {V„ear/Vfar) should be known accurately. 

A miniature version of the KamLAND [HHI and BOREXINO [133 and a 
scaled up version of the CHOOZ three - concentric zone detector design is chosen 
for the construction of the spectrometers fFigure 1331) . KR2DET plans a 4.7 m 
diameter liquid scintillator target, enclosed in transparent spherical balloon. 
The target is viewed by ~800 8-inch EMI-9350 (9350 - 9356) photomultipliers 
trough ^90 cm layer of mineral oil of the zone-2 of the detector. The PMTs of 
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this type have successfully been used in the CHOOZ experiment and are used 
now in the BOREXINO and SNO detectors [123 . A 20% light collection and 
150 - 200 photoelectron signal is expected for 1 MeV positron energy deposition. 
The PMTs are mounted on the stainless steel screen, which separates external 
zone-3 from the central zones of the detector. The ~75 cm thick zone-3 is filled 
with mineral oil (or liquid scintillator) and serves as active (muon) and passive 
shielding from the external radioactivity. 

The ratio of measured positron spectra S[Ef.) far I S{Ee)near in Equation H27|l 
can be slightly distorted because of relative difference in response functions of 
the two "identical" spectrometers. 

The goal of calibration procedures is to measure this difference and introduce 
necessary corrections. This can be done by a combination of different methods. 
First there will be a periodic control of the energy scales in many points using 
7-sources shown by arrows in Figure E3 A useful continuous monitoring of the 
scales at 2.23 MeV can provide neutrons produced by trough going muons and 
captured by the target protons during veto time. 

The second method uses small spontaneous fission ^^^Cf or ^38^ sources 
periodically placed in the detectors. These sources generates continuous energy 
spectrum due to prompt fission gammas and neutron recoils (the dashed line 
in Figure IS7I ). Deviation from unity of the measured spectra can be used to 
calculate relevant corrections. 

The goal is that the systematic uncertainty due to detector spectrometric 
difference essential for Analysis I can be controlled down to 0.5%. 

In Analysis II the systematic uncertainty in the quantity {Lfar/LnearY ' 
(Vnear/Vfar) ' {cnear / £ f ar) in Equation H28|l Can hopcfuUy be kept within 0.8%. 

The choice of the scintillator has not been made so far. There should be 
progress in manufacturing Gd (~0.9 g/liter) loaded scintillators to improve the 
response to neutrons and suppress accidentals, which originate from U/Th gam- 
mas coming from surrounding rock. The Palo Verde Gd-scintillator showed bet- 
ter stability than the scintillator used in CHOOZ. The LENS project considers 
scintillators with rare earth contents as high as ~50 g/liter. One possibility is 
scintillator without Gd based on the mixture of isoparaffin or mineral oil and 
pseudocumene (~20%) with ~2g/liter PPO as primary flour. This scintillator 
has C/H ratio 1.85, density 0.85 kg/hter and 0.785 x 10^^ H atoms per ton. 

The neutrino events satisfy the following requirements: (i) a time window on 
the delay between and neutron signals 2—600 /is, (ii) energy window for the 
neutron candidate 1.7—3.1 MeV and for e"*" 1.2—8.0 MeV, (iii) distance between 
e+ and neutron less than 100 cm. At this stage no pulse shape analysis to reject 
proton recoils is planned. 

Under these assumptions neutrino detection efficiency of 75% was found and 
neutrino detection rate N{e~^,n) — 55/day calculated for the far detector. 

The time correlated background 0.1 per day per one target ton was found 
by extrapolation of the value 0.25/per day per target ton measured at CHOOZ: 

CHOOZ (300 mwe), 0.25/day • ton ^ Kr2DET (600 mwe), 0.1/dayton (29) 
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The accidental coincidences come from the internal radioactivity of detec- 
tor materials and U and Th contained in the surrounding rock. The internal 
component of the background was estimated to be less 0.3/day, which is an 
order of magnitude smaller than the rate of the correlated background (see hep- 
ph/0109277). In contrast to the KamLAND and Borexino experiments three 
orders higher concentrations of U, Th, K and Rn can be tolerated in the liquids 
used in the Kr2DET case. 

First estimations of accidentals coming from the radioactivity of the rock 
showed however that external passive shielding of the detector should be in- 
creased in case scintillator without Gd is used as the neutrino target. 

Calculated neutrino detection rates N{e~^,n) and backgrounds for scintilla- 
tor with no Gd are summarized in Table [TBI 



Table 16: 



Detector 


Distance, 


Target, 


N{e+,n), 


N{e+,n), 


Backgi 


., day ^ 




m 


mass, ton 


day-i 


year~^* 


correl. 


accid.** 


Far 


1000 


46 


55 


16.5- 10^ 


5 


-0.3 


Near 


115 


46 


4200 


12.5-10^ 


5 


-0.3 



* 300 days/year at full power. 



** due to internal radioactivity of the detector materials only. 

Expected 90% CL constraints on the oscillation parameters (Figure I56L 
curves K2Det) are obtained for 40000 detected 17^ in the far detector (750 days 
of full power). The systematic uncertainties (Jshape= 0.5% in the Analysis I 
("shape") and Urate— 0.8% in the Analysis II ("rate") have been assumed. The 
"shape" analysis is somewhat more sensitive and can shift (at Am^ = 2.5 x 10"'^ 
eV^ the SIT? 29 upper limit from 0.14 (CHOOZ) to 0.017. 

The one reactor - two detector approach fully eliminates uncertainties asso- 
ciated with the reactor neutrino source inherent to the absolute method used 
at CHOOZ. 

Small relative difference in conceptually identical detector properties can be 
minimized through calibration and monitoring procedures. 

The detector backgrounds can be measured during reactor OFF periods, 
which periodically follow 50 day long reactor ON periods. 

Good signal to background ratio can be achieved due to sufficiently deep 
underground position of the detectors. 

High statistics can be accumulated in reasonably short time period using 
detectors with —45 ton targets, which are relatively small if compared to modern 
neutrino detectors. 

Neutrino community has accumulated positive experience in building and 
running 3 concentric zone detectors similar to the Kr2DET detectors. 
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G.l Krasnoyarsk site details 

The reactor belongs to the Federal State-Owned Unitary Enterprise MINING & 
CHEMICAL COMBINE (MCC) 53, Lenin Str., Zhelezhnogorsk, Krasnoyarsk 
Territory, RUSSIA, 660972. 

The Krasnoyarsk neutrino laboratory is built in the MCC underground ter- 
ritory. 

There are two places to install the detectors. One of them at ~115 m from 

the reactor is 10 m high 15 x 15 m square room. The other is a 125 m long, 11.5 
high and 15 m wide corridor at ~1000 m from the reactor. More information 
on neutrino at Krasnoyarsk can soon be found at 

http: / / www.lngs.infn.it /site/oxppro/panagic/soction_indcxos / framc_particlcs.html 
(click "Laboratories and experiments" , then "Underground and underwater lab- 
oratories" and go to "Krasnoyarsk neutrino laboratory" ) 

Zhelezhnogorsk is located at about 70km from Krasnoyarsk on the bank 
of the Yenisei River. Zhelezhnogorsk is a very nice and clean town built in 
direct neighborhood to the Siberian taiga, rich of birds and animals. There is 
a beautiful large lake in the center of the town. Picturesque hills surround the 
town center. A musicale theater, hotel, rest home, restaurants, a lot of shops 
are in Zhelezhnogorsk.. The weather is comfortable; the number of sunny days 
is the same as in resort Sochi (at the Black Sea). Winter is cold but not so 
much compared with Moscow, air is dry. The summer and autumn are warmer 
and sunnier than in Moscow. 

Some information about tourism in Krasnoyarsk Territories is available at 
the site: http: / /tlcom.krs.ru/kalinka/indexe.htm, 
tours http : / /tlcom. krs.ru/kalinka/indexe. htm 

Every day there are flights from Moscow to Krasnoyarsk airport. Big com- 
fortable airbus IL86 in 4.5 hours time brings you from Moscow to Krasnoyarsk 
with good service of KrasAir company and a special minivan in 2 hours carries 
you from Krasnoyarsk airport Yemelianovo directly to the center of Zhelezhno- 
gorsk. 

MINING & CHEMICAL COMBINE has two own rest homes; one of them 

is in the town territory near the forest and another outside of the town not 
far from it on the bank of Yenisei River. Both of them have conference halls, 
comfortable living rooms and dining rooms. 
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Krasnoyarsk reactor underground site: £00 mwe 
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Figure 53: Scheme of the KR2DET experiment. 
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Figure 54: Calculated ratio of positron spectra S{Ee) far/S{Ef,)near for some 
oscillation parameters. Values of sin'^(20) are shown at the curves. 
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Figure 56: Reactor antineutrino oscillation plots. Curves "CHOOZ" , 
"KR2DET" (expected) "shape" and "rate" are 90% CL Vf. disappearance lim- 
its. The KR2DET limits are obtained assuming 40,000 detected antineutrinos 
in the far detector, 10:1 effect to background ratio and systematic uncertainties 
o'shape = 0.5% and arate = 0.8%. The shaded area represents the most probable 
atmospheric neutrino mass parameter region. 
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Visible energy, MeV 



Figure 57: Positron visible energy spectrum. 
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